Liu Xudong, Zhu Ruyu, Yang Yongsheng, Xue Ying, Xu Dingguo
College of Chemistry, Key Lab of Green Chemistry and Technology in Ministry of Education, Sichuan University, Chengdu, 610064, People's Republic of China.
School of Pharmacy and Institute of Pharmacy, North Sichuan Medical College, Nanchong, 637100, People's Republic of China.
J Mol Model. 2024 Jul 18;30(8):274. doi: 10.1007/s00894-024-06079-9.
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.
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.
环丙烷骨架因其独特的结构特性,在生物活性分子中发挥着重要作用。这引发了环丙烷衍生物立体选择性合成领域的浓厚兴趣和深入探索。在本研究中,使用密度泛函理论(DFT)计算,研究了在水相和甲苯介质中,3-重氮基-N-甲基吲哚(R1)与两种缺电子烯烃(R2和R3)通过无催化剂的[2 + 1]环丙烷化反应,非对映发散合成环丙烷衍生物的立体选择性的机理和起源。研究结果表明,这些[2 + 1]环加成反应分两个阶段进行,其中第一步不仅是速率决定步骤,而且对产物的立体选择性起着关键作用。计算得到的非对映体比例与实验结果一致。此外,通过利用基于分子力场的非共价相互作用(NCI)分析和能量分解分析(EDA-FF),我们阐明了第一步过渡态TS1s中反应物片段之间的静电相互作用是决定立体选择性的主要因素,这与实验假设的空间位阻和π-π堆积相互作用相反。
在M06-2X(D3)/SMD/6-31 + G(d,p)理论水平下,使用DFT方法对溶剂水和甲苯中势能面(PES)上的所有极小值和过渡态的几何结构进行了全优化。基于溶液中优化后的几何结构,在M06-2X(D3)/6-311 + G(d,p)水平下进行单点能量计算。所有DFT计算均使用Gaussian 09软件进行。使用CYLview软件对优化后的分子结构进行可视化。使用Multiwfn和VMD软件进行NCI分析。Multiwfn程序还用于CDFT和EDA-FF分析。