Pericyclic reaction benchmarks: hierarchical computations targeting CCSDT(Q)/CBS and analysis of DFT performance

Explainable chemical artificial intelligence from accurate machine learning of real-space chemical descriptors

Machine-learned computational chemistry has led to a paradoxical situation in which molecular properties can be accurately predicted, but they are difficult to interpret. Explainable AI (XAI) tools can be used to analyze complex models, but they are highly dependent on the AI technique and the origin of the reference data. Alternatively, interpretable real-space tools can be employed directly, but they are often expensive to compute. To address this dilemma between explainability and accuracy, we developed SchNet4AIM, a SchNet-based architecture capable of dealing with local one-body (atomic) and two-body (interatomic) descriptors. The performance of SchNet4AIM is tested by predicting a wide collection of real-space quantities ranging from atomic charges and delocalization indices to pairwise interaction energies. The accuracy and speed of SchNet4AIM breaks the bottleneck that has prevented the use of real-space chemical descriptors in complex systems. We show that the group delocalization indices, arising from our physically rigorous atomistic predictions, provide reliable indicators of supramolecular binding events, thus contributing to the development of Explainable Chemical Artificial Intelligence (XCAI) models.

Quantum Tunneling Instability in Pericyclic Reactions

Several cycloreversion reactions of the retro-Diels-Alder type were computationally assessed to understand their quantum tunneling (QT) reactivity. N2, CO, and other leaving groups were considered based on their strong exothermicity, as it reduces their thermodynamic and kinetic stabilities. Our results indicate that for many of these reactions, it is essential to take into account their QT decomposition rate, which can massively weaken their molecular stability and shorten their half-lives even at deep cryogenic temperatures. In practical terms, this indicates that many supposedly stable molecules will actually be unsynthesizable or unisolable, and therefore trying to prepare or detect them would be a futile attempt. In addition, we discuss the importance of tunneling to correctly understand the enthalpy of activation and the collective atomic effect on the tunneling kinetic isotope effects to test if third-row atoms can tunnel in a chemical reaction. This project raises the question of the importance of in silico chemistry to guide in vitro chemistry, especially in cases where the latter cannot solve its own uncertainties.

Controlling Diastereoselectivity in Dearomatizing Diels-Alder Reactions of Nitroarenes with 2-Trimethylsilyloxycyclohexadiene

Dearomative Diels-Alder cycloadditions between nitroarenes and 2-trimethylsilyloxycyclohexadiene are carried out under high pressure at room temperature in the absence of any chemical promoter. Reactions are performed with different arenes, including the highly aromatic naphthalenes and quinolines. They lead to 3D-scaffolds with exquisite exo-diastereoselectivity. The exo approach is characterized by lower distortion of the substrates in a late TS and by more favorable orbital interactions presumably between the nitro group and the dienic part, explaining the stereoselectivity.

Enantiodetermining processes in the synthesis of alanine, serine, and isovaline

In this study, quantum chemical calculations were used to explore the synthesis of three chiral α-amino acids, specifically alanine, serine, and isovaline, from reactants found in interstellar space. Our focus is on the crucial step in the synthesis pathway that determines the chirality of the amino acids. The results indicate that in the case of alanine, the determination of enantiomer is primarily influenced by the direction of the collision of molecules or functional groups, which leads to the formation of a chirality center in a crucial intermediate. However, contrary to chemical expectations, the enantiodetermining/enantioselection step for serine and isovaline synthesis occurs prior to the creation of a chirality center.

(3+2)-Cycloadditions of Levoglucosenone (LGO) with Fluorinated Nitrile Imines Derived from Trifluoroacetonitrile: An Experimental and Computational Study

The in situ-generated N-aryl nitrile imines derived from trifluoroacetonitrile smoothly undergo (3+2)-cycloadditions onto the enone fragment of the levoglucosenone molecule, yielding the corresponding, five-membered cycloadducts. In contrast to the 'classic' C(Ph),N(Ph) nitrile imine, reactions with fluorinated C(CF3),N(Ar) analogues lead to stable pyrazolines in a chemo- and stereoselective manner. Based on the result of X-ray single crystal diffraction analysis, their structures were established as exo-cycloadducts with the location of the N-Ar terminus of the 1,3-dipole at the α-position of the enone moiety. The DFT computation demonstrated that the observed reaction pathway results from the strong dominance of kinetic control over thermodynamic control.

Net Intermolecular Silyloxypyrone-Based (5+2) Cycloadditions Utilizing Amides as Enabling and Cleavable Tethers

Silyloxypyrone-based (5+2) cycloadditions were facilitated by amides that allowed for increased reactivity and a pathway for cleaving the tether to afford net intermolecular cycloadducts. Various amides underwent facile cycloaddition, and several experiments revealed steric and electronic factors that accelerate the reaction. tert-Butyl amides reacted faster than less hindered variants in multiple cases. In the case of dearomative oxidopyrylium-indole (5+2) cycloadditions, an amine-based tether was ineffective, whereas amides enabled this powerful transformation. Theoretical calculations evidenced a concerted asynchronous reaction in which the amide facilitates a conformational driving force enabling cycloaddition. Finally, a one-pot acylation/(5+2) cycloaddition/nucleophilic lactam opening and other examples of tosyl lactam opening of a modified cycloadduct were demonstrated.

Not antiaromaticity gain, but increased asynchronicity enhances the Diels-Alder reactivity of tropone

Tropone is an unreactive diene in normal electron demand Diels-Alder reactions, but it can be activated via carbonyl umpolung by using hydrazone ion analogs. Recently, the higher reactivity of hydrazone ion analogs was ascribed to a raised HOMO energy induced by antiaromaticity (L. J. Karas, A. T. Campbell, I. V. Alabugin and J. I. Wu, Org. Lett., 2020, 22, 7083). We show that this is incorrect, and that the activation barrier is lowered by increased asynchronicity.

Molecular Understanding and Practical In Silico Catalyst Design in Computational Organocatalysis and Phase Transfer Catalysis-Challenges and Opportunities

Through the lens of organocatalysis and phase transfer catalysis, we will examine the key components to calculate or predict catalysis-performance metrics, such as turnover frequency and measurement of stereoselectivity, via computational chemistry. The state-of-the-art tools available to calculate potential energy and, consequently, free energy, together with their caveats, will be discussed via examples from the literature. Through various examples from organocatalysis and phase transfer catalysis, we will highlight the challenges related to the mechanism, transition state theory, and solvation involved in translating calculated barriers to the turnover frequency or a metric of stereoselectivity. Examples in the literature that validated their theoretical models will be showcased. Lastly, the relevance and opportunity afforded by machine learning will be discussed.

Reaction profiles for quantum chemistry-computed [3 + 2] cycloaddition reactions

Bio-orthogonal click chemistry based on [3 + 2] dipolar cycloadditions has had a profound impact on the field of biochemistry and significant effort has been devoted to identify promising new candidate reactions for this purpose. To gauge whether a prospective reaction could be a suitable bio-orthogonal click reaction, information about both on- and off-target activation and reaction energies is highly valuable. Here, we use an automated workflow, based on the autodE program, to compute over 5000 reaction profiles for [3 + 2] cycloadditions involving both synthetic dipolarophiles and a set of biologically-inspired structural motifs. Based on a succinct benchmarking study, the B3LYP-D3(BJ)/def2-TZVP//B3LYP-D3(BJ)/def2-SVP level of theory was selected for the DFT calculations, and standard conditions and an (aqueous) SMD model were imposed to mimic physiological conditions. We believe that this data, as well as the presented workflow for high-throughput reaction profile computation, will be useful to screen for new bio-orthogonal reactions, as well as for the development of novel machine learning models for the prediction of chemical reactivity more broadly.

How Bases Catalyze Diels-Alder Reactions

We have quantum chemically studied the base-catalyzed Diels-Alder (DA) reaction between 3-hydroxy-2-pyrone and N-methylmaleimide using dispersion-corrected density functional theory. The uncatalyzed reaction is slow and is preceded by the extrusion of CO₂ via a retro-DA reaction. Base catalysis, for example, by triethylamine, lowers the reaction barrier up to 10 kcal mol^-1 , causing the reaction to proceed smoothly at low temperature, which quenches the expulsion of CO₂ , yielding efficient access to polyoxygenated natural compounds. Our activation strain analyses reveal that the base accelerates the DA reaction via two distinct electronic mechanisms: i) by the HOMO-raising effect, which enhances the normal electron demand orbital interaction; and ii) by donating charge into 3-hydroxy-2-pyrone which accumulates in its reactive region and promotes strongly stabilizing secondary electrostatic interactions with N-methylmaleimide.

Connections between the accuracy of rotational constants and equilibrium molecular structures

Rotational spectroscopy is the technique of choice for investigating molecular structures in the gas phase. Indeed, rotational constants are strongly connected to the geometry of the molecular system under consideration. Therefore, they are powerful tools for assessing the accuracy that quantum chemical approaches can reach in structural determinations. In this review article, it is shown how it is possible to measure the accuracy of a computed equilibrium geometry based on the comparison of rotational constants. But, it is also addressed what accuracy is required by computations for providing molecular structures and thus rotational constants that are useful to experiment. Quantum chemical methodologies for obtaining the "0.1% accuracy" for rotational constants are reviewed for systems ranging in size from small molecules to small polycyclic aromatic hydrocarbons. This accuracy for systems containing two dozen or so atoms opens the way towards future applications such as the accurate characterization of non-covalent interactions, which play a key role in several biological and technological processes.

Automated Variable Electric-Field DFT Application for Evaluation of Optimally Oriented Electric Fields on Chemical Reactivity

Recent theoretical work and experiments at molecular junctions have provided a strong conceptualization for the effects of oriented electric fields (OEFs) on organic reactions. Depending on the axis of application, OEFs can increase (or decrease) the reaction rate or distinguish between isomeric pathways. Despite the conceptual elegance of OEFs, which may be applied externally or induced locally, as tools for catalyzing organic reactions, implementation in synthetically relevant systems has been hampered by inefficiencies in evaluating reaction sensitivity to field effects. Herein, we describe the development of the Automated Variable Electric-Field DFT Application (A.V.E.D.A.) for streamlined evaluation of a reaction's susceptibility to OEFs. This open-source software was designed to be accessible for nonexpert users of computational and programming tools. Following initiation by a single command (and with no subsequent intervention) the Linux workflow manages a series of density functional theory calculations and mathematical manipulations to optimize local-minimum and transition-state structures in oriented electric fields of increasing magnitude. The resulting molecular and reaction dipole moments, field-perturbed geometries, and net effective activation energies are compiled for user interpretation. Ten representative pericyclic reactions that showcase the development and evaluation of A.V.E.D.A. are described.

Mechanism and the Origins of Periselectivity in Cycloaddition Reactions of Benzyne with Dienes

Density functional theory calculations have been used to explore the reaction mechanism of (4 + 2) and (2 + 2) cycloadditions of benzyne with classical dienes. The results indicate the following: (1) (4 + 2) products arise via concerted pathways, (2) (2 + 2) products arise via stepwise pathways with diradical intermediates, and (3) these diradical intermediates are formed via isomerization of carbene intermediates. The origins of periselectivity in these reactions are analyzed using distortion/interaction analysis for the key steps, and they indicate that the tiny distortion in the very early [4 + 2] transition structure, coupled with an entropic favorability, controls selective (4 + 2) cycloaddition.