Theoretical insights into enantioselective [2 + 1] cyclopropanation reactions of diazo compounds with electron-deficient olefins

• Published in: Journal of molecular modeling 2024-08, Vol.30 (8), p.274, Article 274
• Main Authors: Liu, Xudong, Zhu, Ruyu, Yang, Yongsheng, Xue, Ying, Xu, Dingguo
• Format: Article
• Language: English
• Subjects: Alkenes; Characterization and Evaluation of Materials; Chemical synthesis; Chemistry; Chemistry and Materials Science; Computer Appl. in Life Sciences; Computer Applications in Chemistry; Cycloaddition; Cyclopropane; Decomposition reactions; Enantiomers; Mathematical analysis; Molecular Medicine; Molecular structure; Original Paper; Potential energy; Software; Stereoselectivity; Steric hindrance; Theoretical and Computational Chemistry; Toluene
• ISSN: 1610-2940; 0948-5023
• DOI: 10.1007/s00894-024-06079-9

Context The cyclopropane skeleton plays a significant role in bioactive  molecules due to its distinctive structural properties. This has sparked keen interest and in-depth exploration in the field of stereoselective synthesis of cyclopropane derivatives. In the present study, the mechanism and the origin of stereoselectivity of diastereodivergent synthesis of cyclopropane derivatives via the catalyst-free [2 + 1]-cyclopropanation reactions of 3-diazo- N -methylindole ( R1 ) with two types of electron-deficient olefins ( R2 and R3 ) in both aqueous and toluene media have been studied using the DFT calculations. The findings indicate that these [2 + 1] cycloaddition reactions proceed in two stages, where the first step is not only the rate-determining step but also critically dictates the stereoselectivity of the product. The calculated diastereomeric ratios are in agreement with the experimental results. Furthermore, by utilizing non-covalent interaction (NCI) analysis and energy decomposition analysis based on molecular force fields (EDA-FF), we elucidated that the electrostatic interactions between reactant fragments in the transition state TS1s for the first step are the predominant factors determining the stereoselectivity, as opposed to the experimentally hypothesized steric hindrance and π-π stacking interactions. Methods The geometrical structures of all minima and transition states on the potential energy surface (PES) in solvents water and toluene were fully optimized using the DFT method at the M06-2X(D3)/SMD/6–31 + G(d,p) level of theory. Single-point energy calculations were carried out based on the optimized geometries in the solution at the M06-2X(D3)/6–311 + G(d,p) level. All the DFT calculations were performed using the Gaussian 09 software. The optimized molecular structures were visualized using CYLview software. NCI analysis was performed using the Multiwfn and VMD softwares. The Multiwfn program was also used for CDFT and EDA-FF analyses.