TS-tools: Rapid and automated localization of transition states based on a textual reaction SMILES input

Here, TS-tools is presented, a Python package facilitating the automated localization of transition states (TS) based on a textual reaction SMILES input. TS searches can either be performed at xTB or DFT level of theory, with the former yielding guesses at marginal computational cost, and the latter directly yielding accurate structures at greater expense. On a benchmarking dataset of mono- and bimolecular reactions, TS-tools reaches an excellent success rate of 95% already at xTB level of theory. For tri- and multimolecular reaction pathways - which are typically not benchmarked when developing new automated TS search approaches, yet are relevant for various types of reactivity, cf. solvent- and autocatalysis and enzymatic reactivity - TS-tools retains its ability to identify TS geometries, though a DFT treatment becomes essential in many cases. Throughout the presented applications, a particular emphasis is placed on solvation-induced mechanistic changes, another issue that received limited attention in the automated TS search literature so far.

A human-machine interface for automatic exploration of chemical reaction networks

Autonomous reaction network exploration algorithms offer a systematic approach to explore mechanisms of complex chemical processes. However, the resulting reaction networks are so vast that an exploration of all potentially accessible intermediates is computationally too demanding. This renders brute-force explorations unfeasible, while explorations with completely pre-defined intermediates or hard-wired chemical constraints, such as element-specific coordination numbers, are not flexible enough for complex chemical systems. Here, we introduce a STEERING WHEEL to guide an otherwise unbiased automated exploration. The STEERING WHEEL algorithm is intuitive, generally applicable, and enables one to focus on specific regions of an emerging network. It also allows for guiding automated data generation in the context of mechanism exploration, catalyst design, and other chemical optimization challenges. The algorithm is demonstrated for reaction mechanism elucidation of transition metal catalysts. We highlight how to explore catalytic cycles in a systematic and reproducible way. The exploration objectives are fully adjustable, allowing one to harness the STEERING WHEEL for both structure-specific (accurate) calculations as well as for broad high-throughput screening of possible reaction intermediates.

Toward De Novo Catalyst Discovery: Fast Identification of New Catalyst Candidates for Alcohol-Mediated Morita-Baylis-Hillman Reactions

Recently we have demonstrated how a genetic algorithm (GA) starting from random tertiary amines can be used to discover a new and efficient catalyst for the alcohol-mediated Morita-Baylis-Hillman (MBH) reaction. In particular, the discovered catalyst was shown experimentally to be eight times more active than DABCO, commonly used to catalyze the MBH reaction. This represents a breakthrough in using generative models for catalyst optimization. However, the GA procedure, and hence discovery, relied on two important pieces of information; 1) the knowledge that tertiary amines catalyze the reaction and 2) the mechanism and reaction profile for the catalyzed reaction, in particular the transition state structure of the rate-determining step. Thus, truly de novo catalyst discovery must include these steps. Here we present such a method for discovering catalyst candidates for a specific reaction while simultaneously proposing a mechanism for the catalyzed reaction. We show that tertiary amines and phosphines are potential catalysts for the MBH reaction by screening 11 molecular templates representing common functional groups. The method relies on an automated reaction discovery workflow using meta-dynamics calculations. Combining this method for catalyst candidate discovery with our GA-based catalyst optimization method results in an algorithm for truly de novo catalyst discovery.

Sequential On-Surface Cyclodehydrogenation in a Nonplanar Nanographene

On-surface synthesis has emerged as an attractive method for the atomically precise synthesis of new molecular nanostructures, being complementary to the widespread approach based on solution chemistry. It has been particularly successful in the synthesis of graphene nanoribbons and nanographenes. In both cases, the target compound is often generated through cyclodehydrogenation reactions, leading to planarization and the formation of hexagonal rings. To improve the flexibility and tunability of molecular units, however, the incorporation of other, nonbenzenoid, subunits is highly desirable. In this letter, we thoroughly analyze sequential cyclodehydrogenation reactions with a custom-designed molecular precursor. We demonstrate the step-by-step formation of hexagonal and pentagonal rings from the nonplanar precursor within fjord and cove regions, respectively. Computer models comprehensively support the experimental observations, revealing that both reactions imply an initial hydrogen abstraction and a final [1,2] hydrogen shift, but the formation of a pentagonal ring proceeds through a radical mechanism.

Benchmark of general-purpose machine learning-based quantum mechanical method AIQM1 on reaction barrier heights

Artificial intelligence-enhanced quantum mechanical method 1 (AIQM1) is a general-purpose method that was shown to achieve high accuracy for many applications with a speed close to its baseline semiempirical quantum mechanical (SQM) method ODM2*. Here, we evaluate the hitherto unknown performance of out-of-the-box AIQM1 without any refitting for reaction barrier heights on eight datasets, including a total of ∼24 thousand reactions. This evaluation shows that AIQM1's accuracy strongly depends on the type of transition state and ranges from excellent for rotation barriers to poor for, e.g., pericyclic reactions. AIQM1 clearly outperforms its baseline ODM2* method and, even more so, a popular universal potential, ANI-1ccx. Overall, however, AIQM1 accuracy largely remains similar to SQM methods (and B3LYP/6-31G* for most reaction types) suggesting that it is desirable to focus on improving AIQM1 performance for barrier heights in the future. We also show that the built-in uncertainty quantification helps in identifying confident predictions. The accuracy of confident AIQM1 predictions is approaching the level of popular density functional theory methods for most reaction types. Encouragingly, AIQM1 is rather robust for transition state optimizations, even for the type of reactions it struggles with the most. Single-point calculations with high-level methods on AIQM1-optimized geometries can be used to significantly improve barrier heights, which cannot be said for its baseline ODM2* method.

Guided discovery of chemical reaction pathways with imposed activation

Computational power and quantum chemical methods have improved immensely since computers were first applied to the study of reactivity, but the de novo prediction of chemical reactions has remained challenging. We show that complex reaction pathways can be efficiently predicted in a guided manner using chemical activation imposed by geometrical constraints of specific reactive modes, which we term imposed activation (IACTA). Our approach is demonstrated on realistic and challenging chemistry, such as a triple cyclization cascade involved in the total synthesis of a natural product, a water-mediated Michael addition, and several oxidative addition reactions of complex drug-like molecules. Notably and in contrast with traditional hand-guided computational chemistry calculations, our method requires minimal human involvement and no prior knowledge of the products or the associated mechanisms. We believe that IACTA will be a transformational tool to screen for chemical reactivity and to study both by-product formation and decomposition pathways in a guided way.

Best-Practice DFT Protocols for Basic Molecular Computational Chemistry

Nowadays, many chemical investigations are supported by routine calculations of molecular structures, reaction energies, barrier heights, and spectroscopic properties. The lion's share of these quantum-chemical calculations applies density functional theory (DFT) evaluated in atomic-orbital basis sets. This work provides best-practice guidance on the numerous methodological and technical aspects of DFT calculations in three parts: Firstly, we set the stage and introduce a step-by-step decision tree to choose a computational protocol that models the experiment as closely as possible. Secondly, we present a recommendation matrix to guide the choice of functional and basis set depending on the task at hand. A particular focus is on achieving an optimal balance between accuracy, robustness, and efficiency through multi-level approaches. Finally, we discuss selected representative examples to illustrate the recommended protocols and the effect of methodological choices.

Theoretical calculations of formation and reactivity of o-quinomethide derivatives of resorcin[4]arene with reference to empirical data

This paper describes theoretical reaction pathways of alkoxybenzyl derivatives of resorcin[4]arene leading to the formation of o-quinomethide derivatives of resorcin[4]arene (o-QMR[4]A). For each case, the activation energies for the formation of one o-QMR[4]A unit and the activation energies for the backward reaction were calculated. Based on the calculated reaction pathways, the reaction mechanism of o-QMR[4]A formation was proposed. Using the example of o-QMR[4]A generated from a methoxy derivative of resorcin[4]arene, the activation energies with selected nucleophiles were calculated and the reaction mechanisms discussed. Reaction path calculations were performed using the nudged elastic band method and semiempirical extended tight-binding method (GFN2-xTB). Using hydroxybenzyl derivatives of resorcin[4]arene as an example, a comparison of calculated activation energies by selected density-functional theory methods with GFN2-xTB and B97-3c geometries was performed. B97-3c and wB97XD methods were used to calculate the energies of the reactants (R), transition states (TS) and products (P) of the analysed reactions. Theoretical reaction mechanisms were discussed with respect to the orbital-weighted Fukui dual descriptor (Δfw ) and experimental data.

Chemoton 2.0: Autonomous Exploration of Chemical Reaction Networks

Fueled by advances in hardware and algorithm design, large-scale automated explorations of chemical reaction space have become possible. Here, we present our approach to an open-source, extensible framework for explorations of chemical reaction mechanisms based on the first-principles of quantum mechanics. It is intended to facilitate reaction network explorations for diverse chemical problems with a wide range of goals such as mechanism elucidation, reaction path optimization, retrosynthetic path validation, reagent design, and microkinetic modeling. The stringent first-principles basis of all algorithms in our framework is key for the general applicability that avoids any restrictions to specific chemical systems. Such an agile framework requires multiple specialized software components of which we present three modules in this work. The key module, Chemoton, drives the exploration of reaction networks. For the exploration itself, we introduce two new algorithms for elementary-step searches that are based on Newton trajectories. The performance of these algorithms is assessed for a variety of reactions characterized by a broad chemical diversity in terms of bonding patterns and chemical elements. Chemoton successfully recovers the vast majority of these. We provide the resulting data, including large numbers of reactions that were not included in our reference set, to be used as a starting point for further explorations and for future reference.

Bio-inspired chemical space exploration of terpenoids

Many computational methods are devoted to rapidly generating pseudo-natural products to expand the open-ended border of chemical spaces for natural products. However, the accessibility and chemical interpretation were often ignored or underestimated in conventional library/fragment-based or rule-based strategies, thus hampering experimental synthesis. Herein, a bio-inspired strategy (named TeroGen) is developed to mimic the two key biosynthetic stages (cyclization and decoration) of terpenoid natural products, by utilizing physically based simulations and deep learning models, respectively. The precision and efficiency are validated for different categories of terpenoids, and in practice, more than 30 000 sesterterpenoids (10 times as many as the known sesterterpenoids) are predicted to be linked in a reaction network, and their synthetic accessibility and chemical interpretation are estimated by thermodynamics and kinetics. Since it could not only greatly expand the chemical space of terpenoids but also numerate plausible biosynthetic routes, TeroGen is promising for accelerating heterologous biosynthesis, bio-mimic and chemical synthesis of complicated terpenoids and derivatives.

A Theoretical Study on the Structural, Electronic, and Magnetic Properties of Bimetallic Pt13−nNin (N = 0, 3, 6, 9, 13) Nanoclusters to Unveil the Catalytic Mechanisms for the Water-Gas Shift Reaction

In this work, first-principles calculations by using density functional theory at the GFN-xTB level, are performed to investigate the relative stability and structural, electronic, and magnetic properties of bimetallic Pt13−nNin (n = 0, 3, 6, 9, 13) nanoclusters by using corrected Hammer and Nørskov model. In addition, by employing the reaction path and the energetic span models, the energy profile and the turnover frequency are calculated to disclose the corresponding reaction mechanism of the water-gas shift reaction catalyzed by these nanoclusters. Our findings render that Ni causes an overall shrinking of the nanocluster’s size and misalignment of the spin channels, increasing the magnetic nature of the nanoclusters. Pt7Ni6 nanocluster is the most stable as a result of the better coupling between the Pt and Ni d-states. Pt4Ni9 maintains its structure over the reaction cycle, with a larger turnover frequency value than Pt7Ni6. On the other hand, despite Pt10Ni3 presenting the highest value of turnover frequency, it suffers a strong structural deformation over the completion of a reaction cycle, indicating that the catalytic activity can be altered.

Fast and automated identification of reactions with low barriers using meta-MD simulations

We test our meta-molecular dynamics (MD) based approach for finding low-barrier (<30 kcal/mol) reactions on uni- and bimolecular reactions extracted from the barrier dataset developed by Grambow, Pattanaik & Green (2020). For unimolecular reactions the meta-MD simulations identify 25 of the 26 products found by Grambow, Pattanaik & Green (2020), while the subsequent semiempirical screening eliminates an additional four reactions due to an overestimation of the reaction energies or estimated barrier heights relative to DFT. In addition, our approach identifies 36 reactions not found by Grambow, Pattanaik & Green (2020), 10 of which are <30 kcal/mol. For bimolecular reactions the meta-MD simulations identify 19 of the 20 reactions found by Grambow, Pattanaik & Green (2020), while the subsequent semiempirical screening eliminates an additional reaction. In addition, we find 34 new low-barrier reactions. For bimolecular reactions we found that it is necessary to “encourage” the reactants to go to previously undiscovered products, by including products found by other MD simulations when computing the biasing potential as well as decreasing the size of the molecular cavity in which the MD occurs, until a reaction is observed. We also show that our methodology can find the correct products for two reactions that are more representative of those encountered in synthetic organic chemistry. The meta-MD hyperparameters used in this study thus appear to be generally applicable to finding low-barrier reactions.

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

A synergistic approach that combines machine learning with semi-empirical methods enables the fast and accurate prediction of DFT-quality reaction barriers, with mechanistic insights available from semi-empirical transition state geometries.

Chemical reactivity prediction: current methods and different application areas

The ability to predict chemical reactivity of a molecule is highly desirable in drug discovery, both ex vivo (synthetic route planning, formulation, stability) and in vivo: metabolic reactions determine pharmacodynamics, pharmacokinetics and potential toxic effects, and early assessment of liabilities is vital to reduce attrition rates in later stages of development. Quantum mechanics offer a precise description of the interactions between electrons and orbitals in the breaking and forming of new bonds. Modern algorithms and faster computers have allowed the study of more complex systems in a punctual and accurate fashion, and answers for chemical questions around stability and reactivity can now be provided. Through machine learning, predictive models can be built out of descriptors derived from quantum mechanics and cheminformatics, even in the absence of experimental data to train on. In this article, current progress on computational reactivity prediction is reviewed: applications to problems in drug design, such as modelling of metabolism and covalent inhibition, are highlighted and unmet challenges are posed.