Rationalizing the Influence of the Binding Affinity on the Activity of Phosphoserine Phosphatases

The Sabatier principle states that catalytic activity can be maximized when the substrate binding affinity is neither too strong nor too weak. Recent studies have shown that the activity of several hydrolases is maximized at intermediate values of the binding affinity (Michaelis-Menten constant: Km ). However, it remains unclear whether this concept of artificial catalysis is applicable to enzymes in general, especially for those which have evolved under different reaction environments. Herein, we show that the activity of phosphoserine phosphatase is also enhanced at an intermediate Km value of approximately 0.5 mM. Within our dataset, the variation of Km by three orders of magnitude accounted for a roughly 18-fold variation in the activity. Owing to the high phylogenetic and physiological diversity of our dataset, our results support the importance of optimizing Km for enzymes in general. On the other hand, a 77-fold variation in the activity was attributed to other physicochemical parameters, such as the Arrhenius prefactor of kcat , and could not be explained by the Sabatier principle. Therefore, while tuning the binding affinity according to the Sabatier principle is an important consideration, the Km value is only one of many physicochemical parameters which must be optimized to maximize enzymatic activity.

Automated de Novo Design of Olefin Metathesis Catalysts: Computational and Experimental Analysis of a Simple Thermodynamic Design Criterion

Methods for computational de novo design of inorganic molecules have paved the way for automated design of homogeneous catalysts. Such studies have so far relied on correlation-based prediction models as fitness functions (figures of merit), but the soundness of these approaches has yet to be tested by experimental verification of de novo-designed catalysts. Here, a previously developed criterion for the optimization of dative ligands L in ruthenium-based olefin metathesis catalysts RuCl2(L)(L')(═CHAr), where Ar is an aryl group and L' is a phosphine ligand dissociating to activate the catalyst, was used in de novo design experiments. These experiments predicted catalysts bearing an N-heterocyclic carbene (L = 9) substituted by two N-bound mesityls and two tert-butyl groups at the imidazolidin-2-ylidene backbone to be promising. Whereas the phosphine-stabilized precursor assumed by the prediction model could not be made, a pyridine-stabilized ruthenium alkylidene complex (17) bearing carbene 9 was less active than a known leading pyridine-stabilized Grubbs-type catalyst (18, L = H2IMes). A density functional theory-based analysis showed that the unsubstituted metallacyclobutane (MCB) intermediate generated in the presence of ethylene is the likely resting state of both 17 and 18. Whereas the design criterion via its correlation between the stability of the MCB and the rate-determining barrier indeed seeks to stabilize the MCB, it relies on RuCl2(L)(L')(═CH2) adducts as resting states. The change in resting state explains the discrepancy between the prediction and the actual performance of catalyst 17. To avoid such discrepancies and better address the multifaceted challenges of predicting catalytic performance, future de novo catalyst design studies should explore and test design criteria incorporating information from more than a single relative energy or intermediate.

Constructing and interpreting volcano plots and activity maps to navigate homogeneous catalyst landscapes

Volcano plots and activity maps are powerful tools for studying homogeneous catalysis. Once constructed, they can be used to estimate and predict the performance of a catalyst from one or more descriptor variables. The relevance and utility of these tools has been demonstrated in several areas of catalysis, with recent applications to homogeneous catalysts having been pioneered by our research group. Both volcano plots and activity maps are built from linear free energy scaling relationships that connect the value of a descriptor variable(s) with the relative energies of other catalytic cycle intermediates/transition states. These relationships must be both constructed and postprocessed appropriately to obtain the resulting plots/maps; this process requires careful execution to obtain meaningful results. In this protocol, we provide a step-by-step guide to building volcano plots and activity maps using curated reaction profile data. The reaction profile data are obtained using density functional theory computations to model the catalytic cycle. In addition, we provide volcanic, a Python code that automates the steps of the process following data acquisition. Unlike the computation of individual reaction energy profiles, our tools lead to a holistic view of homogeneous catalyst performance that can be broadly applied for both explanatory and screening purposes.

Uncovering the Activity of Alkaline Earth Metal Hydrogenation Catalysis Through Molecular Volcano Plots

Recent advances in alkaline earth (Ae) metal hydrogenation catalysis have broadened the spectrum of potential catalysts to include candidates from the main group, providing a sustainable alternative to the commonly used transition metals. Although Ae-amides have already been demonstrated to catalyze hydrogenation of imines and alkenes, a lucid understanding of how different metal/ligand combinations influence the catalytic activity is yet to be established. In this article, we use linear scaling relationships and molecular volcano plots to assess the potential of the Ae metal-based catalysts for the hydrogenation of alkenes. By analyzing combinations of eight metals (mono-, bi-, tri-, and tetravalent) and seven ligands, we delineate the impact of metal-ligand interplay on the hydrogenation activity. Our findings highlight that the catalytic activity is majorly determined by the charge and the size of the metal ions. While bivalent Ae metal cations delicately regulate the binding and the release of the reactants and the products, respectively, providing the right balance for this reaction, ligands play only a minor role in determining their catalytic activity. We show how volcano plots can be utilized for the rapid screening of prospective Ae catalysts to establish a guideline to achieve maximum activity in facilitating the hydrogenation process.

SUPPLEMENTARY INFORMATION

The online version of this article at 10.1007/s11244-021-01480-7.

Trends in computational molecular catalyst design

Computational methods have emerged as a powerful tool to augment traditional experimental molecular catalyst design by providing useful predictions of catalyst performance and decreasing the time needed for catalyst screening. In this perspective, we discuss three approaches for computational molecular catalyst design: (i) the reaction mechanism-based approach that calculates all relevant elementary steps, finds the rate and selectivity determining steps, and ultimately makes predictions on catalyst performance based on kinetic analysis, (ii) the descriptor-based approach where physical/chemical considerations are used to find molecular properties as predictors of catalyst performance, and (iii) the data-driven approach where statistical analysis as well as machine learning (ML) methods are used to obtain relationships between available data/features and catalyst performance. Following an introduction to these approaches, we cover their strengths and weaknesses and highlight some recent key applications. Furthermore, we present an outlook on how the currently applied approaches may evolve in the near future by addressing how recent developments in building automated computational workflows and implementing advanced ML models hold promise for reducing human workload, eliminating human bias, and speeding up computational catalyst design at the same time. Finally, we provide our viewpoint on how some of the challenges associated with the up-and-coming approaches driven by automation and ML may be resolved.

The Genesis of Molecular Volcano Plots

ConspectusFor the past two decades, linear free energy scaling relationships and volcano plots have seen frequent use as computational tools that aid in understanding and predicting the catalytic behavior of heterogeneous and electrocatalysts. Based on Sabatier's principle, which states that a catalyst should bind a substrate neither too strongly nor too weakly, volcano plots provide an estimate of catalytic performance (e.g., overpotential, catalytic cycle thermodynamics/kinetics, etc.) through knowledge of a descriptor variable. By the use of linear free energy scaling relationships, the value of this descriptor is employed to estimate the relative energies of other catalytic cycle intermediates/transition states. Postprocessing of these relationships leads to a volcano curve that reveals the anticipated performance of each catalyst, with the best species appearing on or near the peak or plateau. While the origin of volcanoes is undoubtedly rooted in examining heterogeneously catalyzed reactions, only recently has this concept been transferred to the realm of homogeneous catalysis. This Account summarizes the work done by our group in implementing and refining "molecular volcano plots" for use in analyzing and predicting the behavior of homogeneous catalysts.We begin by taking the reader through the initial proof-of-principle study that transferred the model from heterogeneous to homogeneous catalysis by examining thermodynamic aspects of a Suzuki-Miyaura cross-coupling reaction. By establishing linear free energy scaling relationships and reproducing the volcano shape, we definitively showed that volcano plots are also valid for homogeneous systems. On the basis of this key finding, we further illustrate how unified pictures of C-C cross-coupling thermodynamics were created using three-dimensional molecular volcanoes.The second section highlights an important transformation from "thermodynamic" to "kinetic" volcanoes by using the descriptor variable to directly estimate transition state barriers. Taking this idea further, we demonstrate how volcanoes can be used to directly predict an experimental observable, the turnover frequency. Discussion is also provided on how different flavors of molecular volcanoes can be used to analyze aspects of homogeneous catalysis of interest to experimentalists, such as determining the product selectivity and probing the substrate scope.The third section focuses on incorporating machine learning approaches into molecular volcanoes and invoking big-data-type approaches in the analysis of catalytic behavior. Specifically, we illustrate how machine learning can be used to predict the value of the descriptor variable, which facilitates nearly instantaneous screening of thousands of catalysts. With the large amount of data created from the machine learning/volcano plot tandem, we show how the resulting database can be mined to garner an enhanced understanding of catalytic processes. Emphasis is also placed on the latest generation of augmented volcano plots, which differ fundamentally from earlier volcanoes by elimination of the use of linear free energy scaling relationships and by assessment of the similarity of the complete catalytic cycle energy profile to that for an ideal reference species that is used to discriminate catalytic performance.We conclude by examining a handful of applications of molecular volcano plots to interesting problems in homogeneous catalysis and offering thoughts on the future prospects and uses of this new set of tools.

Scaling Relationships and Volcano Plots in Homogeneous Catalysis

Scaling relations and volcano plots are widely used in heterogeneous catalysis. In this Perspective, we discuss the prospects and challenges associated with the application of similar concepts in homogeneous catalysis using examples from the literature that have appeared recently.

Data-powered augmented volcano plots for homogeneous catalysis

Given the computational resources available today, data-driven approaches can propel the next leap forward in catalyst design. Using a data-driven inspired workflow consisting of data generation, statistical analysis, and dimensionality reduction algorithms we explore trends surrounding the thermodynamics of a model hydroformylation reaction catalyzed by group 9 metals bearing phosphine ligands. Specifically, we introduce "augmented volcano plots" as a means to easily visualize the similarity of each catalyst's complete catalytic cycle energy profile to that of a hypothetical ideal reference profile without relying upon linear scaling relationships. In addition to quickly identifying catalysts that most closely match the ideal thermodynamic catalytic cycle energy profile, these maps also enable a more refined comparison of closely lying species in standard volcano plots. For the reaction studied here, they inherently uncover the presence of multiple sets of scaling relationships differentiated by metal type, where iridium catalysts follow distinct relationships from cobalt/rhodium catalysts and have profiles that more closely match the ideal thermodynamic profile. Reconstituted molecular volcano plots confirm the findings of the augmented volcanoes by showing that hydroformylation thermodynamics are governed by two distinct volcano shapes, one for iridium catalysts and a second for cobalt/rhodium species.

Scaling relationships and volcano plots of homogeneous transition metal catalysis

This Frontier article highlights the recent applications of linear scaling relationships and volcano plots in homogeneous transition metal catalysis.

Machine learning meets volcano plots: computational discovery of cross-coupling catalysts

The application of modern machine learning to challenges in atomistic simulation is gaining attraction. We present new machine learning models that can predict the energy of the oxidative addition process between a transition metal complex and a substrate for C-C cross-coupling reactions. In turn, this quantity can be used as a descriptor to estimate the activity of homogeneous catalysts using molecular volcano plots. The versatility of this approach is illustrated for vast libraries of organometallic catalysts based on Pt, Pd, Ni, Cu, Ag, and Au combined with 91 ligands. Out-of-sample machine learning predictions were made on a total of 18 062 compounds leading to 557 catalyst candidates falling into the ideal thermodynamic window. This number was further refined by searching for candidates with an estimated price lower than 10 US$ per mmol. The 37 catalyst finalists are dominated by palladium phosphine ligand combinations but also include the earth abundant transition metal (Cu) with less common ligands. Our results indicate that modern statistical learning techniques can be applied to the computational discovery of readily available and promising catalyst candidates.