Electrostatic origin of the catalytic effect of a supramolecular host catalyst

Computational Screening of Putative Catalyst Transition Metal Complexes as Guests in a Ga4L612- Nanocage

Metal-organic cages form well-defined microenvironments that can enhance the catalytic proficiency of encapsulated transition metal complexes (TMCs). We introduce a screening protocol to efficiently identify TMCs that are promising candidates for encapsulation in the Ga4L612- nanocage. We obtain TMCs from the Cambridge Structural Database with geometric and electronic characteristics amenable to encapsulation and mine the text of associated manuscripts to curate TMCs with documented catalytic functionality. By docking candidate TMCs inside the nanocage cavity and carrying out electronic structure calculations, we identify a subset of successfully optimized candidates (TMC-34) and observe that encapsulated guests occupy an average of 60% of the cavity volume, in line with previous observations. Notably, some guests occupy as much as 72% of the cavity as a result of linker rotation. Encapsulation has a universal effect on the electrostatic potential (ESP), systematically decreasing the ESP at the metal center of each TMC in the TMC-34 data set, while minimally altering TMC metal partial charges. Collectively these observations support geometry-based screening of potential guests and suggest that encapsulation in Ga4L612- cages could electrostatically stabilize diverse cationic or electropositive intermediates. We highlight candidate guests with associated known reactivity and solubility most amenable for encapsulation in experimental follow-up studies.

Origins of High-Activity Cage-Catalyzed Michael Addition

Cage catalysis continues to create significant interest, yet catalyst function remains poorly understood. Herein, we report mechanistic insights into coordination-cage-catalyzed Michael addition using kinetic and computational methods. The study has been enabled by the detection of identifiable catalyst intermediates, which allow the evolution of different cage species to be monitored and modeled alongside reactants and products. The investigations show that the overall acceleration results from two distinct effects. First, the cage reaction shows a thousand-fold increase in the rate constant for the turnover-limiting C-C bond-forming step compared to a reference state. Computational modeling and experimental analysis of activation parameters indicate that this stems from a significant reduction in entropy, suggesting substrate coencapsulation. Second, the cage markedly acidifies the bound pronucleophile, shifting this equilibrium by up to 6 orders of magnitude. The combination of these two factors results in accelerations up to 109 relative to bulk-phase reference reactions. We also show that the catalyst can fundamentally alter the reaction mechanism, leading to intermediates and products that are not observable outside of the cage. Collectively, the results show that cage catalysis can proceed with very high activity and unique selectivity by harnessing a series of individually weak noncovalent interactions.

Picking the lock of coordination cage catalysis

The design principles of metallo-organic assembly reactions have facilitated access to hundreds of coordination cages of varying size and shape. Many of these assemblies possess a well-defined cavity capable of hosting a guest, pictorially mimicking the action of a substrate binding to the active site of an enzyme. While there are now a growing collection of coordination cages that show highly proficient catalysis, exhibiting both excellent activity and efficient turnover, this number is still small compared to the vast library of metal-organic structures that are known. In this review, we will attempt to unpick and discuss the key features that make an effective coordination cage catalyst, linking structure to activity (and selectivity) using lessons learnt from both experimental and computational analysis of the most notable exemplars. We will also provide an outlook for this area, reasoning why coordination cages have the potential to become the gold-standard in (synthetic) non-covalent catalysis.

Manipulation of N-heterocyclic carbene reactivity with practical oriented electric fields

Electric fields can be used to tune the nucleophilicity and electrophilicity of N-heterocyclic carbenes and enhance their catalytic activity.

Catalysis by [Ga4 L6 ]12- Metallocage on the Nazarov Cyclization: The Basicity of Complexed Alcohol is Key

The Nazarov cyclization is investigated in solution and within K₁₂ [Ga₄ L₆ ] supramolecular organometallic cage by means of computational methods. The reaction needs acidic condition in solution but works at neutral pH in the presence of the metallocage. The reaction steps for the process are analogous in both media: (a) protonation of the alcohol group, (b) water loss and (c) cyclization. The relative Gibbs energies of all the steps are affected by changing the environment from solvent to the metallocage. The first step in the mechanism, the alcohol protonation, turns out to be the most critical one for the acceleration of the reaction inside the metallocage. In order to calculate the relative stability of protonated alcohol inside the cavity, we propose a computational scheme for the calculation of basicity for species inside cavities and can be of general use. These results are in excellent agreement with the experiments, identifying key steps of catalysis and providing an in-depth understanding of the impact of the metallocage on all the reaction steps.

Source of Rate Acceleration for Carbocation Cyclization in Biomimetic Supramolecular Cages

The results of quantum chemical and molecular dynamics calculations reveal that polyanionic gallium-based cages accelerate cyclization reactions of pentadienyl alcohols as a result of substrate cage interactions, preferential binding of reactive conformations of substrate/H₃O⁺ pairs, and increased substrate basicity. However, the increase in basicity dominates. Experimental structure-activity relationship studies in which the metal vertices and overall charge of the cage are varied confirm the model derived via calculations.

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.

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.

Origin of the Rate Acceleration in the C-C Reductive Elimination from Pt(IV)-complex in a [Ga4 L6 ]12- Supramolecular Metallocage

The reductive elimination on [(Me₃ P)₂ Pt(MeOH)(CH₃ )₃ ]⁺ , 2P, complex performed in MeOH solution and inside a [Ga₄ L₆ ]^12- metallocage are computationally analysed by mean of QM and MD simulations and compared with the mechanism of gold parent systems previously reported [Et₃ PAu(MeOH)(CH₃ )₂ ]⁺ , 2Au. The comparative analysis between the encapsulated Au(III) and Pt(IV)-counterparts shows that there are no additional solvent MeOH molecules inside the cavity of the metallocage for both systems. The Gibbs energy barriers for the 2P reductive elimination calculated at DFT level are in good agreement with the experimental values for both environments. The effect of microsolvation and encapsulation on the rate acceleration are evaluated and shows that the latter is far more relevant, conversely to 2Au. Energy decomposition analysis indicates that the encapsulation is the main responsible for most of the energy barrier reduction. Microsolvation and encapsulation effects are not equally contributing for both metal systems and consequently, the reasons of the rate acceleration are not the same for both metallic systems despite the similarity between them.

Modeling Kinetics and Thermodynamics of Guest Encapsulation into the [M4L6]12- Supramolecular Organometallic Cage

The encapsulation of molecular guests into supramolecular hosts is a complex molecular recognition process in which the guest displaces the solvent from the host cavity, while the host deforms to let the guest in. An atomistic description of the association would provide valuable insights on the physicochemical properties that guide it. This understanding may be used to design novel host assemblies with improved properties (i.e., affinities) toward a given class of guests. Molecular simulations may be conveniently used to model the association processes. It is thus of interest to establish efficient protocols to trace the encapsulation process and to predict the associated magnitudes ΔG_bind and ΔG_bind^⧧. Here, we report the calculation of the Gibbs energy barrier and Gibbs binding energy by means of explicit solvent molecular simulations for the [Ga₄L₆]^12– metallocage encapsulating a series of cationic molecules. The ΔG_bind^⧧ for encapsulation was estimated by means of umbrella sampling simulations. The steps involved were identified, including ion-pair formation and naphthalene rotation (from L ligands of the metallocage) during the guest’s entrance. The ΔG_bind values were computed using the attach–pull–release method. The results reveal the sensitivity of the estimates on the force field parameters, in particular on atomic charges, showing that higher accuracy is obtained when charges are derived from implicit solvent quantum chemical calculations. Correlation analysis identified some indicators for the binding affinity trends. All computed magnitudes are in very good agreement with experimental observations. This work provides, on one side, a benchmarked way to computationally model a highly charged metallocage encapsulation process. This includes a nonstandard parameterization and charge derivation procedure. On the other hand, it gives specific mechanistic information on the binding processes of [Ga₄L₆]^12– at the molecular level where key motions are depicted. Taken together, the study provides an interesting option for the future design of metal–organic cages.

Modeling the Reaction of Carboxylic Acids and Isonitriles in a Self-Assembled Capsule

Quantum chemical calculations were used to study the reaction of carboxylic acids with isonitriles inside a resorcinarene-based self-assembled capsule. Experimentally, it has been shown that the reactions between p-tolylacetic acid and n-butyl isonitrile or isopropyl isonitrile behave differently in the presence of the capsule compared both with each other and also with their solution counterparts. Herein, the reasons for these divergent behaviors are addressed by comparing the detailed energy profiles for the reactions of the two isonitriles inside and outside the capsule. An energy decomposition analysis was conducted to quantify the different factors affecting the reactivity. The calculations reproduce the experimental findings very well. Thus, encapsulation leads to lowering of the energy barrier for the first step of the reaction, the concerted α-addition and proton transfer, which in solution is rate-determining, and this explains the rate acceleration observed in the presence of the capsule. The barrier for the final step of the reaction, the 1,3 O→N acyl transfer, is calculated to be higher with the isopropyl substituent inside the capsule compared with n-butyl. With the isopropyl substituent, the transition state and the product of this step are significantly shorter than the preceding intermediate, and this results in energetically unfavorable empty spaces inside the capsule, which cause a higher barrier. With the n-butyl substituent, on the other hand, the carbon chain can untwine and hence uphold an appropriate guest length.

Gold-Catalyzed Cross-Coupling Reactions: An Overview of Design Strategies, Mechanistic Studies, and Applications

Transition-metal-catalyzed cross-coupling reactions are central to many organic synthesis methodologies. Traditionally, Pd, Ni, Cu, and Fe catalysts are used to promote these reactions. Recently, many studies have showed that both homogeneous and heterogeneous Au catalysts can be used for activating selective cross-coupling reactions. Here, an overview of the past studies, current trends, and future directions in the field of gold-catalyzed coupling reactions is presented. Design strategies to accomplish selective homocoupling and cross-coupling reactions under both homogeneous and heterogeneous conditions, computational and experimental mechanistic studies, and their applications in diverse fields are critically reviewed. Specific topics covered are: oxidant-assisted and oxidant-free reactions; strain-assisted reactions; dual Au and photoredox catalysis; bimetallic synergistic reactions; mechanisms of reductive elimination processes; enzyme-mimicking Au chemistry; cluster and surface reactions; and plasmonic catalysis. In the relevant sections, theoretical and computational studies of Au^I /Au^III chemistry are discussed and the predictions from the calculations are compared with the experimental observations to derive useful design strategies.

Microsolvation and Encapsulation Effects on Supramolecular Catalysis: C-C Reductive Elimination inside [Ga4L6]12- Metallocage

The host effect of the supramolecular [Ga₄L₆]^12- tetrahedral metallocage on reductive elimination of substrate by encapsulated Au(III) complexes is investigated by means of computational methods. The behavior of the reactants in solution and within the metallocage is initially evaluated by means of classical molecular dynamics simulations. These results guided the selection of proper computational models to describe the reaction in solution and inside the metallocage at the DFT level. The calculated Gibbs energy barriers are in very good agreement with experiment both in solution and inside the metallocage. The analysis in solution revealed that microsolvation around the Au(III) complex increases the Gibbs energy barrier. The analysis within the metallocage shows that its encapsulation favors the reaction. The process can be formally described as removing explicit microsolvation around the gold complex and encapsulating the metal complex inside the metallocage. Both processes are important for the reaction, but the removal of the solvent molecules surrounding the Au(III) metal complex is fundamental for the reduction of the reaction barrier. The energy decomposition analysis of the barrier among strain, interaction, and thermal terms shows that strain term is very low whereas the contribution of thermal (entropic) effects is moderate. Interestingly, the key term responsible for reducing the Gibbs energy barrier is the interaction. This term can be mainly associated with electrostatic interactions in agreement with previous examples in the literature.

Electrostatics Generated by a Supramolecular Capsule Stabilizes the Transition State for Carbon-Carbon Reductive Elimination from Gold(III) Complex

Tetrahedral assemblies of stoichiometry M₄L₆ have been proven to catalyze a range of chemical reactions including the carbon-carbon reductive elimination reaction from transition metals such as gold. Here, we perform quantum chemical calculations of Gold(III) transition metal complexes in vacuum, and encapsulated in Ga₄L₆^12- or Si₄L₆^8- assemblies within both a reaction field continuum solvent and in an aqueous molecular environment with counterions, to rationalize the rate enhancements observed experimentally for the reductive elimination reaction. We find that the Ga₄L₆^12- assembly lowers the energy barrier of the reaction compared to Si₄L₆^8-, which is consistent with kinetic trends observed experimentally. We have determined that the primary factor for catalytic rate acceleration stems from the electrostatic environment emanating from the Ga₄L₆^12- capsule as opposed to the water or counterions.

Deconvoluting the Role of Charge in a Supramolecular Catalyst

We have demonstrated that the microenvironment of a highly anionic supramolecular catalyst can mimic the active sites of enzymes and impart rate accelerations of a million-fold or more. However, these microenvironments can be challenging to study, especially in the context of understanding which specific features of the catalyst are responsible for its high performance. We report here the development of an experimental mechanistic probe consisting of two isostructural catalysts. When examined in parallel transformations, the behavior of these catalysts provides insight relevant to the importance of anionic host charge on reactivity. These two catalysts exhibit similar host-substrate interactions, but feature a significant difference in overall anionic charge (12^- and 8^-). Within these systems, we compare the effect of constrictive binding in a net neutral aza-Cope rearrangement. We then demonstrate how the magnitude of anionic host charge has an exceptional influence on the reaction rates for a Nazarov cyclization, evidenced by an impressive 680-fold change in reaction rate as a consequence of a 33% reduction in catalyst charge.

Quantum Chemical Modeling of Cycloaddition Reaction in a Self-Assembled Capsule

Dispersion-corrected density functional theory is used to study the cycloaddition reaction between phenyl acetylene and phenyl azide inside a synthetic, self-assembled capsule. The capsule is first characterized computationally and a previously unrecognized structure is identified as being the most stable. Next, an examination of the free energies of host-guest complexes is conducted, considering all possible reagent, solvent, and solvent impurity combinations as guests. The experimentally observed relative stabilities of host-guest complexes are quite well reproduced, when the experimental concentrations are taken into account. Experimentally, the presence of the host capsule has been shown to accelerate the cycloaddition reaction and to yield exclusively the 1,4-regioisomer product. Both these observations are reproduced by the calculations. A detailed energy decomposition analysis shows that reduction of the entropic cost of bringing together the reactants along with a geometric destabilization of the reactant supercomplex are the major contributors to the rate acceleration compared to the background reaction. Finally, a sensitivity analysis is conducted to assess the stability of the results with respect to the choice of methodology.

Quantum mechanical calculation of electric fields and vibrational Stark shifts at active site of human aldose reductase

Recent advance in biophysics has made it possible to directly measure site-specific electric field at internal sites of proteins using molecular probes with C = O or C≡N groups in the context of vibrational Stark effect. These measurements directly probe changes of electric field at specific protein sites due to, e.g., mutation and are very useful in protein design. Computational simulation of the Stark effect based on force fields such as AMBER and OPLS, while providing good insight, shows large errors in comparison to experimental measurement due to inherent difficulties associated with point charge based representation of force fields. In this study, quantum mechanical calculation of protein's internal electrostatic properties and vibrational Stark shifts was carried out by using electrostatically embedded generalized molecular fractionation with conjugate caps method. Quantum calculated change of mutation-induced electric field and vibrational Stark shift is reported at the internal probing site of enzyme human aldose reductase. The quantum result is in much better agreement with experimental data than those predicted by force fields, underscoring the deficiency of traditional point charge models describing intra-protein electrostatic properties.

Computer aided enzyme design and catalytic concepts

Gaining a deeper understanding of enzyme catalysis is of great practical and fundamental importance. Over the years it has become clear that despite advances made in experimental mutational studies, a quantitative understanding of enzyme catalysis will not be possible without the use of computer modeling approaches. While we believe that electrostatic preorganization is by far the most important catalytic factor, convincing the wider scientific community of this may require the demonstration of effective rational enzyme design. Here we make the point that the main current advances in enzyme design are basically advances in directed evolution and that computer aided enzyme design must involve approaches that can reproduce catalysis in well-defined test cases. Such an approach is provided by the empirical valence bond method.

Supramolecular catalysis. Part 2: artificial enzyme mimics

The design of artificial catalysts able to compete with the catalytic proficiency of enzymes is an intense subject of research. Non-covalent interactions are thought to be involved in several properties of enzymatic catalysis, notably (i) the confinement of the substrates and the active site within a catalytic pocket, (ii) the creation of a hydrophobic pocket in water, (iii) self-replication properties and (iv) allosteric properties. The origins of the enhanced rates and high catalytic selectivities associated with these properties are still a matter of debate. Stabilisation of the transition state and favourable conformations of the active site and the product(s) are probably part of the answer. We present here artificial catalysts and biomacromolecule hybrid catalysts which constitute good models towards the development of truly competitive artificial enzymes.