Computational molecular refinement to enhance enantioselectivity by reinforcing hydrogen bonding interactions in major reaction pathway

Computational analyses have revealed that the distortion of a catalyst and the substrates and their interactions are key to determining the stability of the transition state. Hence, two strategies "distortion strategy" and "interaction strategy" can be proposed for improving enantiomeric excess in enantioselective reactions. The "distortion strategy" is used as a conventional approach that destabilizes the TS (transition state) of the minor pathway. On the other hand, the "interaction strategy" focuses on the stabilization of the TS of the major pathway in which an enhancement of the reaction rate is expected. To realize this strategy, we envisioned the TS stabilization of the major reaction pathway by reinforcing hydrogen bonding and adopted the chiral phosphoric acid-catalysed enantioselective Diels-Alder reaction of 2-vinylquinolines with dienylcarbamates. The intended "interaction strategy" led to remarkable improvements in the enantioselectivity and reaction rate.

Chiral Phosphoric Acid Catalyzed Conversion of Epoxides into Thiiranes: Mechanism, Stereochemical Model, and New Catalyst Design

Computations and experiments leading to new chiral phosphoric acids (CPAs) for epoxide thionations are reported. Density functional theory calculations reveal the mechanism and origin of the enantioselectivity of such CPA-catalyzed epoxide thionations. The calculated mechanistic information was used to design new efficient CPAs that were tested experimentally and found to be highly effective. Bulky ortho-substituents on the 3,3'-aryl groups of the CPA are important to restrict the position of the epoxide in the key transition states for the enantioselectivity-determining step. Larger para-substituents significantly improve the enantioselectivity of the reaction.

Nazarov Cyclizations Catalyzed by BINOL Phosphoric Acid Derivatives: Quantum Chemistry Struggles To Predict the Enantioselectivity

Quantum chemical calculations have successfully predicted the stereoselectivities of many BINOL phosphoric acid catalyzed reactions over the past 10-15 years. Herein we report a contrasting example: a reaction for which standard quantum chemistry techniques have proven unexpectedly ineffective at explaining the stereoselectivity. The Nazarov cyclizations of a divinyl ketone catalyzed by a BINOL phosphoric acid or H₈-BINOL dithiophosphoric acid were studied with a conventional contemporary quantum chemical approach, consisting of transition state optimizations with B3LYP-D3(BJ) and single-point calculations with several functionals in implicit solvent. Unexpectedly, different functionals gave widely different predictions of the level of enantioselectivity and were unable even to agree on which enantiomer of the product would predominate. Molecular dynamics simulations with the OPLS-AA force field provided evidence that the transition state geometries optimized with DFT in the gas phase or in implicit solvent are not good representations of the true transition states of these reactions in solution. One possible reason for this, which may also explain the failure of quantum chemical techniques to reliably predict the enantioselectivity, is the fact that the transition states contain ion pairs which are not highly organized and do not contain any strongly directional noncovalent interactions between the substrate and the catalyst.

Expanding the Range of Force Fields Available for ONIOM Calculations: The SICTWO Interface

The ONIOM scheme is one of the most popular QM/MM approaches, but its extended application has been so far hindered by the limited availability of force fields in most practical implementations. This paper describes a simple software code to overcome this limitation, and its application to three representative chemical problems. The "Shell Interface for Combining Tinker With ONIOM" (SICTWO) program gives access to all force fields available in the Tinker molecular mechanics program from a Gaussian09 or Gaussian16 calculation. The first application presented is the geometry optimization of dithienobicyclo-[4.4.1]-undeca-3,8-diene-11-one ethylene glycol ketal. A variety of force fields were tested for the MM part, and the molecular structure using the ONIOM(B3LYP:OPLS-AA) method shows good agreement with the experimental results. The second problem is related with enantioselectivity of the vanadium-catalyzed asymmetric oxidation of 1,2-bis( tert-butyl)-disulfide. ONIOM(B3LYP:MM3) results, obtained with SICTWO, are consistent with those of a previous IMOMM(B3LYP:MM3) study. In the third application, we have combined AMOEBA09 polarizable force field with ONIOM to study a water molecule binding on water ice. Calculated ONIOM(ωB97X-D:AMOEBA09) binding energies are consistent with the ωB97X-D results. These studies show SICTWO as an easy-to-use tool that can further expand the use of the multiscale ONIOM method within computational chemistry and computational biology.

Insights into the N-Heterocyclic Carbene (NHC)-Catalyzed Intramolecular Cyclization of Aldimines: General Mechanism and Role of Catalyst

One of the most challenging questions in the Lewis base organocatalyst field is how to predict the most electrophilic carbon for the complexation of N-heterocyclic carbene (NHC) and reactant. This study provides a valuable case for this issue. Multiple mechanisms (A, B, C, D, and E) for the intramolecular cyclization of aldimine catalyzed by NHC were investigated by using density functional theory (DFT). The computed results reveal that the NHC energetically prefers attacking the iminyl carbon (AIC mode, which is associated with mechanisms A and C) rather than attacking the olefin carbon (AOC mode, which is associated with mechanisms B and D) or attacking the carbonyl carbon (ACC mode, which is associated with mechanism E) of aldimine. The calculated results based on the different reaction models indicate that mechanism A (AIC mode), which is associated with the formation of the aza-Breslow intermediate, is the most favorable pathway. For mechanism A, there are five steps: (1) nucleophilic addition of NHC to the iminyl carbon of aldimine; (2) [1,2]-proton transfer to form an aza-Breslow intermediate; (3) intramolecular cyclization; (4) the other [1,2]-proton transfer; and (5) regeneration of NHC. The analyses of reactivity indexes have been applied to explain the chemoselectivity, and the general principles regarding the possible mechanisms would be useful for the rational design of NHC-catalyzed chemoselective reactions.

Computational insight into the cooperative role of non-covalent interactions in the aza-Henry reaction catalyzed by quinine derivatives: mechanism and enantioselectivity

Density functional theory (DFT) calculations were performed to elucidate the mechanism and the origin of the high enantioselectivity of the aza-Henry reaction of isatin-derived N-Boc ketimine catalyzed by a quinine-derived catalyst (QN). The C-C bond formation step is found to be both the rate-determining and the stereo-controlled step. The results revealed the important role of the phenolic OH group in pre-organizing the complex of nitromethane and QN and stabilizing the in situ-generated nitronate and protonated QN. Three possible activation modes for C-C bond formation involving different coordination patterns of catalyst and substrates were studied, and it was found that both the ion pair-hydrogen bonding mode and the Brønsted acid-hydrogen bonding mode are viable, with the latter slightly preferred for the real catalytic system. The calculated enantiomeric excess (ee) favouring the S enantiomer is in good agreement with the experimental result. The high reactivity and enantioselectivity can be ascribed to the cooperative role of the multiple non-covalent interactions, including classical and non-classical H bonding as well as anionπ interactions. These results also highlight the importance of the inclusion of dispersion correction for achieving a reasonable agreement between theory and experiment for the current reaction.

Insights into Ag(i)-catalyzed addition reactions of amino alcohols to electron-deficient olefins: competing mechanisms, role of catalyst, and origin of chemoselectivity

The competing mechanisms of Ag(i)-catalyzed chemoselective addition reactions of amino alcohols and electron-deficient olefins leading to the O-adduct or N-adduct products were systematically studied with density functional theory methods. Calculations indicate that the AgHMDS/dppe versus AgOAc/dppe catalytic systems can play different roles and thereby generate two different products. The AgHMDS/dppe system works as a Brønsted base to deprotonate the amino alcohol OH to form the Ag–O bond, which leads to formation of the O-adduct. In contrast, the AgOAc/dppe system mainly acts as a Lewis acid to coordinate with O and N atoms of the amino alcohol, but it cannot act as the Brønsted base to further activate the OH group because of its weaker basicity. Therefore, the AgOAc/dppe catalyzed reaction has a mechanism that is similar to the non-catalyzed reaction, and generates the same N-adduct. The obtained insights will be important for rational design of the various kinds of cooperatively catalyzed chemoselective addition reactions, including the use of the less nucleophilic hydroxyl groups of unprotected amino alcohols.

The origin of the chemoselectivities of Ag(i)-catalyzed addition reactions of amino alcohols to olefin has been predicted for the first time.

Computational characterization of the mechanism for the light-driven catalytic trichloromethylation of acylpyridines

DFT and DFT/MM calculations are applied to a photocatalytic enantioselective reaction and shown to be able to characterize the mechanism.

Origins of Selectivity and General Model for Chiral Phosphoric Acid-Catalyzed Oxetane Desymmetrizations

The origins of the high enantioselectivity of chiral phosphoric acid-catalyzed oxetane desymmetrizations were investigated by density functional theory (DFT) calculations. Distortion of the catalyst structure, caused by steric crowding in the catalyst pocket of one enantiomeric transition state, is the main cause for stereochemical preference. A general model was developed to assist in the rational design of new catalysts for related transformations.

QM/MM study on the enantioselectivity of spiroacetalization catalysed by an imidodiphosphoric acid catalyst: how confinement works

An explanation of why confined imidodiphosphoric acid catalyst succeeds where other chiral phosphoric acid catalysts fail.