Prediction of Stereochemistry using Q2MM

Enantiopure molecules form apparently racemic monolayers of chiral cyclic pentamers

Ultra-high vacuum scanning tunneling microscopy (UHV-STM) was used to investigate two related molecules pulse-deposited onto Au(111) surfaces: indoline-2-carboxylic acid and proline (pyrrolidine-2-carboxylic acid). Indoline-2-carboxylic acid and proline form both dimers and C5-symmetric "pinwheel" pentamers. Enantiomerically pure S-(-)-indoline-2-carboxylic acid and S-proline were used, and the pentamer structures observed for both were chiral. However, the presence of apparently equal numbers of 'right-' and 'left-handed' pinwheels is contrary to the general understanding that the chirality of the molecule dictates supramolecular chirality. A variety of computational methods were used to elucidate pentamer geometry for S-proline. Straightforward geometry optimization proved difficult, as the size of the cluster and the number of possible intermolecular interactions produced an interaction potential with multiple local minima. Instead, the Amber force field was used to exhaustively search all of phase space for chemically reasonable pentamer structures, producing a limited number of candidate structures that were then optimized as gas-phase clusters using density functional theory (DFT). The binding energies of the two lowest-energy pentamers on the Au(111) surface were then calculated by plane-wave DFT using the VASP software, and STM images predicted. These calculations indicate that the right- and left-handed pentamers are instead two different polymorphs.

Understanding the Nonlinear Hammett Relationship in Osmylation of Olefins with OsO4-Amine Ligands: Importance of Singlet-Diradical Character

Although the concerted [3 + 2] mechanism of osmium-catalyzed asymmetric dihydroxylation has been generally accepted, the unusual nonlinear Hammett relationship induced by amine-type ligands remains unexplained. To understand this, we carried out a density functional theory (DFT) study for the osmylation of substituted styrenes by the following: OsO4, OsO4-pyridine, OsO4-4-cyanopyridine, OsO4-4-pyrrolidinopyridine, and OsO4-quinuclidine. Calculations using the M06 functional successfully reproduce the experimentally observed nonlinear relationships. The transition states exhibit considerable singlet-diradical character, which causes the nonlinear Hammett relationship. Regardless of the presence or absence of an amine-type ligand, an electron donation from styrene to OsO4 is observed, indicating no mechanistic change. Calculations indicate that the electronic interaction between the amine-type ligand and styrene also influences the reaction rate.

Finding Relevant Retrosynthetic Disconnections for Stereocontrolled Reactions

Machine learning-driven computer-aided synthesis planning (CASP) tools have become important tools for idea generation in the design of complex molecule synthesis but do not adequately address the stereochemical features of the target compounds. A novel approach to automated extraction of templates used in CASP that includes stereochemical information included in the US Patent and Trademark Office (USPTO) and an internal AstraZeneca database containing reactions from Reaxys, Pistachio, and AstraZeneca electronic lab notebooks is implemented in the freely available AiZynthFinder software. Three hundred sixty-seven templates covering reagent- and substrate-controlled as well as stereospecific reactions were extracted from the USPTO, while 20,724 templates were from the AstraZeneca database. The performance of these templates in multistep CASP is evaluated for 936 targets from the ChEMBL database and an in-house selection of 791 AZ designs. The potential and limitations are discussed for four case studies from ChEMBL and examples of FDA-approved drugs.

Overcoming the Pitfalls of Computing Reaction Selectivity from Ensembles of Transition States

The prediction of reaction selectivity is a challenging task for computational chemistry, not only because many molecules adopt multiple conformations but also due to the exponential relationship between effective activation energies and rate constants. To account for molecular flexibility, an increasing number of methods exist that generate conformational ensembles of transition state (TS) structures. Typically, these TS ensembles are Boltzmann weighted and used to compute selectivity assuming Curtin-Hammett conditions. This strategy, however, can lead to erroneous predictions if the appropriate filtering of the conformer ensembles is not conducted. Here, we demonstrate how any possible selectivity can be obtained by processing the same sets of TS ensembles for a model reaction. To address the burdensome filtering task in a consistent and automated way, we introduce marc, a tool for the modular analysis of representative conformers that aids in avoiding human errors while minimizing the number of reoptimization computations needed to obtain correct reaction selectivity.

Polarizability matters in enantio-selection

The prevalence of chirality, or, handedness in biological world is a fundamental phenomenon and a characteristic hallmark of life. Thus, understanding the origin of enantio-selection, i.e., the sense and magnitude of asymmetric induction, has been a long-pursued goal in asymmetric catalysis. Herein, we demonstrated a polarizability-derived electronic effect that was shown to be capable of rationalizing a broad range of stereochemical observations made in the field of asymmetric catalysis. This effect provided a consistent enantio-control model for the prediction of major enantiomers formed in a ruthenium-catalyzed asymmetric transfer hydrogenations of ketones. Direct and quantitative linear free energy relationships between substrates' local polarizabilities and observed enantio-selectivity were also revealed in three widely known asymmetric catalytic systems covering both reductions and oxidations. This broadly applicable polarizability-based electronic effect, in conjunction with conventional wisdom mainly leveraging on steric effect considerations, should aid rational design of enantio-selective processes for better production of chiral substances.

Conformational Sampling over Transition-Metal-Catalyzed Reaction Pathways: Toward Revealing Atroposelectivity

The Py-Conformational-Sampling (PyCoSa) technique is introduced as a systematic computational means to sample the configurational space of transition-metal-catalyzed stereoselective reactions. When applied to atroposelective Suzuki-Miyaura coupling to create axially chiral biaryl products, the results show a range of mechanistic possibilities that include multiple low-energy channels through which C-C bonds can be formed.

Differentiating the Yield of Chemical Reactions Using Parameters in First-Order Kinetic Equations to Identify Elementary Steps That Control the Reactivity from Complicated Reaction Path Networks

The yield of a chemical reaction is obtained by solving its rate equation. This study introduces an approach for differentiating yields by utilizing the parameters of the rate equation, which is expressed as a first-order linear differential equation. The yield derivative for a specific pair of reactants and products is derived by mathematically expressing the rate constant matrix contraction method, which is a simple kinetic analysis method. The parameters of the rate equation are the Gibbs energies of the intermediates and transition states in the reaction path network used to formulate the rate equation. Thus, our approach for differentiating the yield allows a numerical evaluation of the contribution of energy variation to the yield for each intermediate and transition state in the reaction path network. In other words, a comparison of these values automatically extracts the factors affecting the yield from a complicated reaction path network consisting of numerous reaction paths and intermediates. This study verifies the behavior of the proposed approach through numerical experiments on the reaction path networks of a model system and the Rh-catalyzed hydroformylation reaction. Moreover, the possibility of using this approach for designing ligands in organometallic catalysts is discussed.

On Accelerating Substrate Optimization Using Computational Gibbs Energy Barriers: A Numerical Consideration Utilizing a Computational Data Set

Substrate optimization is a time- and resource-consuming step in organic synthesis. Recent advances in chemo- and materials-informatics provide systematic and efficient procedures utilizing tools such as Bayesian optimization (BO). This study explores the possibility of reducing the required experiments further by utilizing computational Gibbs energy barriers. To thoroughly validate the impact of using computational Gibbs energy barriers in BO-assisted substrate optimization, this study employs a computational Gibbs energy barrier data set in the literature and performs an extensive numerical investigation virtually regarding the Gibbs energy barriers as virtual experimental results and those with systematic and random noises as virtual computational results. The present numerical investigation shows that even the computational reactivity affected by noises of as much as 20 kJ/mol helps reduce the number of required experiments.

Machine Learning in Unmanned Systems for Chemical Synthesis

Chemical synthesis is state-of-the-art, and, therefore, it is generally based on chemical intuition or experience of researchers. The upgraded paradigm that incorporates automation technology and machine learning (ML) algorithms has recently been merged into almost every subdiscipline of chemical science, from material discovery to catalyst/reaction design to synthetic route planning, which often takes the form of unmanned systems. The ML algorithms and their application scenarios in unmanned systems for chemical synthesis were presented. The prospects for strengthening the connection between reaction pathway exploration and the existing automatic reaction platform and solutions for improving autonomation through information extraction, robots, computer vision, and intelligent scheduling were proposed.

A Quantum-Guided Molecular Mechanics Force Field for the Ferrocene Scaffold

Ferrocene derivatives have a wide range of applications, including as ligands in asymmetric catalysis, due to their chemical stability, rigid backbone, steric bulk, and ability to encode stereochemical information via planar chirality. Unfortunately, few of the available molecular mechanics force fields incorporate parameters for the accurate study of this important building block. Here, we present a MM3* force field for ferrocenyl ligands, which was generated using the quantum-guided molecular mechanics (Q2MM) method. Detailed validation by comparison to DFT calculations and crystal structures demonstrates the accuracy of the parameters and uncovers the physical origin of deviations through excess energy analysis. Combining the ferrocene force field with a force field for Pd-allyl complexes and comparing the crystal structures shows the compatibility with previously developed MM3* force fields. Finally, the ferrocene force field was combined with a previously published transition-state force field to predict the stereochemical outcomes of the aminations of Pd-allyl complexes with different amines and different chiral ferrocenyl ligands, with an R² of ∼0.91 over 10 examples.

Machine learning and semi-empirical calculations: a synergistic approach to rapid, accurate, and mechanism-based reaction barrier prediction

Modern QM modelling methods, such as DFT, have provided detailed mechanistic insights into countless reactions. However, their computational cost inhibits their ability to rapidly screen large numbers of substrates and catalysts in reaction discovery. For a C-C bond forming nitro-Michael addition, we introduce a synergistic semi-empirical quantum mechanical (SQM) and machine learning (ML) approach that allows the prediction of DFT-quality reaction barriers in minutes, even on a standard laptop using widely available modelling software. Mean absolute errors (MAEs) are obtained that are below the accepted chemical accuracy threshold of 1 kcal mol-1 and substantially better than SQM methods without ML correction (5.71 kcal mol-1). Predictive power is shown to hold when the ML models are applied to an unseen set of compounds from the toxicology literature. Mechanistic insight is also achieved via the generation of full SQM transition state (TS) structures which are found to be very good approximations for the DFT-level geometries, revealing important steric interactions in some TSs. This combination of speed, accuracy, and mechanistic insight is unprecedented; current ML barrier models compromise on at least one of these important criteria.

The (not so) simple prediction of enantioselectivity - a pipeline for high-fidelity computations

The computation of reaction selectivity represents an appealing complementary route to experimental studies and a powerful means to refine catalyst design strategies. Accurately establishing the selectivity of reactions facilitated by molecular catalysts, however, remains a challenging task for computational chemistry. The small free energy differences that lead to large variations in the enantiomeric ratio (er) represent particularly tricky quantities to predict with sufficient accuracy to be helpful for prioritizing experiments. Further complicating this problem is the fact that standard approaches typically consider only one or a handful of conformers identified through human intuition as pars pro toto of the conformational space. Obviously, this assumption can potentially lead to dramatic failures should key energetic low-lying structures be missed. Here, we introduce a multi-level computational pipeline leveraging the graph-based Molassembler library to construct an ensemble of molecular catalysts. The manipulation and interpretation of molecules as graphs provides a powerful and direct route to tailored functionalization and conformer generation that facilitates high-throughput mechanistic investigations of chemical reactions. The capabilities of this approach are validated by examining a Rh(iii) catalyzed asymmetric C-H activation reaction and assessing the limitations associated with the underlying ligand design model. Specifically, the presence of remarkably flexible chiral Cp ligands, which induce the experimentally observed high level of selectivity, present a rich configurational landscape where multiple unexpected conformations contribute to the reported enantiomeric ratios (er). Using Molassembler, we show that considering about 20 transition state conformations per catalysts, which are generated with little human intervention and are not tied to "back-of-the-envelope" models, accurately reproduces experimental er values with limited computational expense.

Automated fitting of transition state force fields for biomolecular simulations

The generation of surrogate potential energy functions (PEF) that are orders of magnitude faster to compute but as accurate as the underlying training data from high-level electronic structure methods is one of the most promising applications of fitting procedures in chemistry. In previous work, we have shown that transition state force fields (TSFFs), fitted to the functional form of MM3* force fields using the quantum guided molecular mechanics (Q2MM) method, provide an accurate description of transition states that can be used for stereoselectivity predictions of small molecule reactions. Here, we demonstrate the applicability of the method for fit TSFFs to the well-established Amber force field, which could be used for molecular dynamics studies of enzyme reaction. As a case study, the fitting of a TSFF to the second hydride transfer in Pseudomonas mevalonii 3-hydroxy-3-methylglutaryl coenzyme A reductase (PmHMGR) is used. The differences and similarities to fitting of small molecule TSFFs are discussed.

Proofreading experimentally assigned stereochemistry through Q2MM predictions in Pd-catalyzed allylic aminations

The palladium-catalyzed enantioselective allylic substitution by carbon or nitrogen nucleophiles is a key transformation that is particularly useful for the synthesis of bioactive compounds. Unfortunately, the selection of a suitable ligand/substrate combination often requires significant screening effort. Here, we show that a transition state force field (TSFF) derived by the quantum-guided molecular mechanics (Q2MM) method can be used to rapidly screen ligand/substrate combinations. Testing of this method on 77 literature reactions revealed several cases where the computationally predicted major enantiomer differed from the one reported. Interestingly, experimental follow-up led to a reassignment of the experimentally observed configuration. This result demonstrates the power of mechanistically based methods to predict and, where necessary, correct the stereochemical outcome.

A Dataset of Computational Reaction Barriers for the Claisen Rearrangement: Chemical and Numerical Analysis

Theoretical reaction screening based on Gibbs energy barriers would be promising to accelerate chemical reactions mining. The number of quantum chemical calculations can be reduced by using an optimization algorithm such as genetic algorithm (GA) and Bayesian optimization (BO). The focus of this study is to generate a dataset of reaction barriers of size ∼100000. Such a dataset would be useful to quickly evaluate various implementations of an optimization algorithm such as GA and BO. The dataset includes Gibbs energy barriers of the Claisen rearrangement for ∼100000 molecules computed on the basis of a semiempirical theory PM7. After evaluating its chemical and numerical features, it is found that the dataset well reflects chemical trends of various substitutions and is useful in testing various implementations of an optimization algorithm. The dataset is available in the supplementary material of this paper.

Computational Discovery of Transition-metal Complexes: From High-throughput Screening to Machine Learning

Transition-metal complexes are attractive targets for the design of catalysts and functional materials. The behavior of the metal-organic bond, while very tunable for achieving target properties, is challenging to predict and necessitates searching a wide and complex space to identify needles in haystacks for target applications. This review will focus on the techniques that make high-throughput search of transition-metal chemical space feasible for the discovery of complexes with desirable properties. The review will cover the development, promise, and limitations of "traditional" computational chemistry (i.e., force field, semiempirical, and density functional theory methods) as it pertains to data generation for inorganic molecular discovery. The review will also discuss the opportunities and limitations in leveraging experimental data sources. We will focus on how advances in statistical modeling, artificial intelligence, multiobjective optimization, and automation accelerate discovery of lead compounds and design rules. The overall objective of this review is to showcase how bringing together advances from diverse areas of computational chemistry and computer science have enabled the rapid uncovering of structure-property relationships in transition-metal chemistry. We aim to highlight how unique considerations in motifs of metal-organic bonding (e.g., variable spin and oxidation state, and bonding strength/nature) set them and their discovery apart from more commonly considered organic molecules. We will also highlight how uncertainty and relative data scarcity in transition-metal chemistry motivate specific developments in machine learning representations, model training, and in computational chemistry. Finally, we will conclude with an outlook of areas of opportunity for the accelerated discovery of transition-metal complexes.

Accelerating the optimization of enzyme-catalyzed synthesis conditions via machine learning and reactivity descriptors

Enzyme-catalyzed synthesis reactions are of crucial importance for a wide range of applications. An accurate and rapid selection of optimal synthesis conditions is crucial and challenging for both human knowledge and computer predictions. In this work, a new scenario, which combines a data-driven machine learning (ML) model with reactivity descriptors, is developed to predict the optimal enzyme-catalyzed synthesis conditions and the reaction yield. Fourteen reactivity descriptors in total are constructed to describe 125 reactions (classified into five categories) included in different reaction mechanisms. Nineteen ML models are developed to train the dataset and the Quadratic support vector machine (SVM) model is found to exhibit the best performance. The Quadratic SVM model is then used to predict the optimal reaction conditions, which are subsequently used to obtain the highest yield among 109 200 reaction conditions with different molar ratios of substrates, solvents, water contents, enzyme concentrations and temperatures for each reaction. The proposed protocol should be generally applicable to a diverse range of chemical reactions and provides a black-box evaluation for optimizing the reaction conditions of organic synthesis reactions.

Stereoselectivity Predictions for the Pd-Catalyzed 1,4-Conjugate Addition Using Quantum-Guided Molecular Mechanics

The conjugate addition of aryl boronic acids to enones is a powerful synthetic tool to introduce quaternary chiral centers, but the experimentally observed stereoselectivities vary widely, and the identification of suitable substrate-ligand combinations requires significant effort. We describe the development and application of a transition-state force field (TSFF) by the quantum-guided molecular mechanics (Q2MM) method that is validated using an automated screen of 9 ligands, 38 aryl boronic acids, and 22 enones, leading to a MUE of 1.8 kJ/mol and a R² value of 0.877 over 82 examples. A detailed error analysis identified the structural origin for the deviations in the small group of outliers. The TSFF was then used to predict the stereoselectivity for 27 ligands and 59 enones. The vast majority of the virtual screening results are in line with the expected results. Selected results for 6-substituted pyrox ligands, which were not part of the training set, were followed up by density functional theory and experimental studies.

Reaction-based machine learning representations for predicting the enantioselectivity of organocatalysts

Hundreds of catalytic methods are developed each year to meet the demand for high-purity chiral compounds. The computational design of enantioselective organocatalysts remains a significant challenge, as catalysts are typically discovered through experimental screening. Recent advances in combining quantum chemical computations and machine learning (ML) hold great potential to propel the next leap forward in asymmetric catalysis. Within the context of quantum chemical machine learning (QML, or atomistic ML), the ML representations used to encode the three-dimensional structure of molecules and evaluate their similarity cannot easily capture the subtle energy differences that govern enantioselectivity. Here, we present a general strategy for improving molecular representations within an atomistic machine learning model to predict the DFT-computed enantiomeric excess of asymmetric propargylation organocatalysts solely from the structure of catalytic cycle intermediates. Mean absolute errors as low as 0.25 kcal mol-1 were achieved in predictions of the activation energy with respect to DFT computations. By virtue of its design, this strategy is generalisable to other ML models, to experimental data and to any catalytic asymmetric reaction, enabling the rapid screening of structurally diverse organocatalysts from available structural information.

Microsecond timescale MD simulations at the transition state of PmHMGR predict remote allosteric residues

Understanding the mechanisms of enzymatic catalysis requires a detailed understanding of the complex interplay of structure and dynamics of large systems that is a challenge for both experimental and computational approaches. More importantly, the computational demands of QM/MM simulations mean that the dynamics of the reaction can only be considered on a timescale of nanoseconds even though the conformational changes needed to reach the catalytically active state happen on a much slower timescale. Here we demonstrate an alternative approach that uses transition state force fields (TSFFs) derived by the quantum-guided molecular mechanics (Q2MM) method that provides a consistent treatment of the entire system at the classical molecular mechanics level and allows simulations at the microsecond timescale. Application of this approach to the second hydride transfer transition state of HMG-CoA reductase from Pseudomonas mevalonii (PmHMGR) identified three remote residues, R396, E399 and L407, (15-27 Å away from the active site) that have a remote dynamic effect on enzyme activity. The predictions were subsequently validated experimentally via site-directed mutagenesis. These results show that microsecond timescale MD simulations of transition states are possible and can predict rather than just rationalize remote allosteric residues.

Organic reactivity from mechanism to machine learning

As more data are introduced in the building of models of chemical reactivity, the mechanistic component can be reduced until 'big data' applications are reached. These methods no longer depend on underlying mechanistic hypotheses, potentially learning them implicitly through extensive data training. Reactivity models often focus on reaction barriers, but can also be trained to directly predict lab-relevant properties, such as yields or conditions. Calculations with a quantum-mechanical component are still preferred for quantitative predictions of reactivity. Although big data applications tend to be more qualitative, they have the advantage to be broadly applied to different kinds of reactions. There is a continuum of methods in between these extremes, such as methods that use quantum-derived data or descriptors in machine learning models. Here, we present an overview of the recent machine learning applications in the field of chemical reactivity from a mechanistic perspective. Starting with a summary of how reactivity questions are addressed by quantum-mechanical methods, we discuss methods that augment or replace quantum-based modelling with faster alternatives relying on machine learning.

Building a Toolbox for the Analysis and Prediction of Ligand and Catalyst Effects in Organometallic Catalysis

Computers have become closely involved with most aspects of modern life, and these developments are tracked in the chemical sciences. Recent years have seen the integration of computing across chemical research, made possible by investment in equipment, software development, improved networking between researchers, and rapid growth in the application of predictive approaches to chemistry, but also a change of attitude rooted in the successes of computational chemistry-it is now entirely possible to complete research projects where computation and synthesis are cooperative and integrated, and work in synergy to achieve better insights and improved results. It remains our ambition to put computational prediction before experiment, and we have been working toward developing the key ingredients and workflows to achieve this.The ability to precisely tune selectivity along with high catalyst activity make organometallic catalysts using transition metal (TM) centers ideal for high-value-added transformations, and this can make them appealing for industrial applications. However, mechanistic variations of TM-catalyzed reactions across the vast chemical space of different catalysts and substrates are not fully explored, and such an exploration is not feasible with current resources. This can lead to complete synthetic failures when new substrates are used, but more commonly we see outcomes that require further optimization, such as incomplete conversion, insufficient selectivity, or the appearance of unwanted side products. These processes consume time and resources, but the insights and data generated are usually not tied to a broader predictive workflow where experiments test hypotheses quantitatively, reducing their impact.These failures suggest at least a partial deviation of the reaction pathway from that hypothesized, hinting at quite complex mechanistic manifolds for organometallic catalysts that are affected by the combination of input variables. Mechanistic deviation is most likely when challenging multifunctional substrates are being used, and the quest for so-called privileged catalysts is quickly replaced by a need to screen catalyst libraries until a new "best" match between the catalyst and substrate can be identified and the reaction conditions can be optimized. As a community we remain confined to broad interpretations of the substrate scope of new catalysts and focus on small changes based on idealized catalytic cycles rather than working toward a "big data" view of organometallic homogeneous catalysis with routine use of predictive models and transparent data sharing.Databases of DFT-calculated steric and electronic descriptors can be built for such catalysts, and we summarize here how these can be used in the mapping, interpretation, and prediction of catalyst properties and reactivities. Our motivation is to make these databases useful as tools for synthetic chemists so that they challenge and validate quantitative computational approaches. In this Account, we demonstrate their application to different aspects of catalyst design and discovery and their integration with computational mechanistic studies and thus describe the progress of our journey toward truly predictive models in homogeneous organometallic catalysis.

Computational methods for training set selection and error assessment applied to catalyst design: guidelines for deciding which reactions to run first and which to run next

Different subset selection methods are examined to guide catalyst selection in optimization campaigns. Error assessment methods are used to quantitatively inform selection of new catalyst candidates from in silico libraries of catalyst structures.

Predicting glycosylation stereoselectivity using machine learning

Predicting the stereochemical outcome of chemical reactions is challenging in mechanistically ambiguous transformations. The stereoselectivity of glycosylation reactions is influenced by at least eleven factors across four chemical participants and temperature. A random forest algorithm was trained using a highly reproducible, concise dataset to accurately predict the stereoselective outcome of glycosylations. The steric and electronic contributions of all chemical reagents and solvents were quantified by quantum mechanical calculations. The trained model accurately predicts stereoselectivities for unseen nucleophiles, electrophiles, acid catalyst, and solvents across a wide temperature range (overall root mean square error 6.8%). All predictions were validated experimentally on a standardized microreactor platform. The model helped to identify novel ways to control glycosylation stereoselectivity and accurately predicts previously unknown means of stereocontrol. By quantifying the degree of influence of each variable, we begin to gain a better general understanding of the transformation, for example that environmental factors influence the stereoselectivity of glycosylations more than the coupling partners in this area of chemical space.

Development of a Computer-Guided Workflow for Catalyst Optimization. Descriptor Validation, Subset Selection, and Training Set Analysis

Modern, enantioselective catalyst development is driven largely by empiricism. Although this approach has fostered the introduction of most of the existing synthetic methods, it is inherently limited by the skill, creativity, and chemical intuition of the practitioner. Herein, we present a complementary approach to catalyst optimization in which statistical methods are used at each stage to streamline development. To construct the optimization informatics workflow, a number of critical components had to be subjected to rigorous validation. First, the critically important molecular descriptors were validated in two case studies to establish the importance of conformation-dependent molecular representations. Next, with a large data set available, it was possible to investigate the amount of data necessary to make predictive models with different modeling methods. Given the commercial availability of many catalyst structures, it was possible to compare models generated with algorithmically selected training sets and commercially available training sets. Finally, the augmentation of limited data sets is demonstrated in a method informed by unsupervised learning to restore the accuracy of the generated models.

Transition State Force Field for the Asymmetric Redox-Relay Heck Reaction

A transition state force field (TSFF) was developed using the quantum-guided molecular mechanics (Q2MM) method to describe the stereodetermining migratory insertion step of the enantioselective redox-relay Heck reaction for a range of multisubstituted alkenes. We show that the TSFF is highly predictive through an external validation of the TSFF against 151 experimentally determined stereoselectivities resulting in an R² of 0.89 and MUE of 1.8 kJ/mol. In addition, limitations in the underlying force field were identified by comparison of the TSFF results to DFT level calculations. A novel application of the TSFF was demonstrated for 31 cases where the enantiomer predicted by the TSFF differed from the originally published values. Experimental determination of the absolute configuration demonstrated that the computational predictions were accurate, suggesting that TSFFs can be used for the rapid prediction of the absolute stereochemistry for a class of reactions. Finally, a virtual ligand screen was conducted utilizing both the TSFF and a simple molecular correlation method. Both methods were similarly predictive, but the TSFF was able to show greater utility through transferability, speed, and interpretability.

Unusual KIE and dynamics effects in the Fe-catalyzed hetero-Diels-Alder reaction of unactivated aldehydes and dienes

Hetero-Diels-Alder (HDA) reaction is an important synthetic method for many natural products. An iron(III) catalyst was developed to catalyze the challenging HDA reaction of unactivated aldehydes and dienes with high selectivity. Here we report extensive density-functional theory (DFT) calculations and molecular dynamics simulations that show effects of iron (including its coordinate mode and/or spin state) on the dynamics of this reaction: considerably enhancing dynamically stepwise process, broadening entrance channel and narrowing exit channel from concerted asynchronous transition states. Also, our combined computational and experimental secondary KIE studies reveal unexpectedly large KIE values for the five-coordinate pathway even with considerable C-C bond forming, due to equilibrium isotope effect from the change in the metal coordination. Moreover, steric and electronic effects are computationally shown to dictate the C=O chemoselectivity for an α,β-unsaturated aldehyde, which is verified experimentally. Our mechanistic study may help design homogeneous, heterogeneous and biological catalysts for this challenging reaction.

Quantitative Structure-Selectivity Relationships in Enantioselective Catalysis: Past, Present, and Future

The dawn of the 21st century has brought with it a surge of research related to computer-guided approaches to catalyst design. In the past two decades, chemoinformatics, the application of informatics to solve problems in chemistry, has increasingly influenced prediction of activity and mechanistic investigations of organic reactions. The advent of advanced statistical and machine learning methods, as well as dramatic increases in computational speed and memory, has contributed to this emerging field of study. This review summarizes strategies to employ quantitative structure-selectivity relationships (QSSR) in asymmetric catalytic reactions. The coverage is structured by initially introducing the basic features of these methods. Subsequent topics are discussed according to increasing complexity of molecular representations. As the most applied subfield of QSSR in enantioselective catalysis, the application of local parametrization approaches and linear free energy relationships (LFERs) along with multivariate modeling techniques is described first. This section is followed by a description of global parametrization methods, the first of which is continuous chirality measures (CCM) because it is a single parameter derived from the global structure of a molecule. Chirality codes, global, multivariate descriptors, are then introduced followed by molecular interaction fields (MIFs), a global descriptor class that typically has the highest dimensionality. To highlight the current reach of QSSR in enantioselective transformations, a comprehensive collection of examples is presented. When combined with traditional experimental approaches, chemoinformatics holds great promise to predict new catalyst structures, rationalize mechanistic behavior, and profoundly change the way chemists discover and optimize reactions.

Holistic prediction of enantioselectivity in asymmetric catalysis

When faced with unfamiliar reaction space, synthetic chemists typically apply reported conditions (reagents, catalyst, solvent, additives) from closely-related reactions to new substrate types. Unfortunately, this approach often fails due to subtle, albeit important, differences in reaction requirements. Consequently, a significant goal in synthetic chemistry is the ability to transfer chemical observations from one reaction to another, quantitatively. Here, we present such a platform by developing a holistic, data-driven workflow for deriving statistical models for one set of reactions that can be applied to predict out-of-sample examples. As a validating case study, published enantioselectivity data sets that employ BINOL-derived chiral phosphoric acids for a range of nucleophilic addition reactions to imines were combined and statistical models developed. These models reveal the general interactions imparting asymmetric induction and allow the quantitative transfer of this information to new reaction components. The disclosed techniques create opportunities for translating comprehensive reaction analysis to diverse chemical space, streamlining both catalyst and reaction development.

After the database of the reactions is constructed, the experimental output, enantiomeric ratios, were mathematically modelled through linear regression techniques to reveal which of the proposed parameters allow for the prediction of new outcomes. The detailed acquisition of parameters can be found in the Supplementary Information and the descriptor tables attached as an accompanying spreadsheet. The models produced were evaluated for their goodness of fit, R², and their robustness is demonstrated by external validations’ goodness of fit, _predR². The nearer the R² and slope values are to one (indicating a tight, one-to-one correlation between predicted and measured outcomes) and the nearer the intercept is to zero (indicating minimal systematic error), the more robust the model. Potential models were refined through number of parameters, because this allows for a mechanistically informative interrogation and cross-validation scores. Leave one reaction out (LORO) analysis was performed to probe general mechanistic principles, which provides the basis for mechanistic transfer of experimental observations and tested further by predicting out-of-sample.