构象预组织在顺式-1,2-二甲磺酰氧基-双(苄氧基)环辛烷反应性中的作用：活化应变视角

《ChemistryOpen》：The Role of Conformational Preorganization in the Reactivity of cis-1,2-Dimesylate-bis(benzyloxy)cyclooctane: An Activation Strain Perspective

　　

为阐明顺式取代环辛烷二甲磺酸酯对叠氮化钠的反应性控制因素，研究人员开展了一项计算研究。计算得到的能量剖面清晰表明，单取代中间体的形成在动力学上具有显著优先性，而第二步SN2取代则对应明显更高的活化能垒。活化应变模型（Activation Strain Model, ASM）分析指出，该能垒升高主要源于更大的畸变能以及更不利的相互作用能。与这些结果一致，非共价相互作用（Non–Covalent Interaction, NCI）分析显示，在较高能量过渡态中，稳定化相互作用减弱且位阻排斥增强。对环辛烷构象的进一步分析表明，两条反应路径在扭转畸变模式上存在差异。较低能量过渡态经由一种更具预组织特征的几何构型进行，而较高能量路径则需要额外的构象重组，才能使离去基达到可反应排布。综合来看，这些结果确立了构象预组织是控制该柔性环体系反应性的关键因素。

该文发表于《ChemistryOpen》，围绕一个实验上表现出强烈单取代选择性的中等环体系亲核取代反应展开，核心问题是：在对顺式-1,2-二甲磺酸酯取代的环辛烷底物进行叠氮化钠亲核取代时，为什么第一步取代能够发生，而第二步在同样条件下却几乎完全被抑制。氨基环醇类（aminocyclitols，多羟基氨基环烷烃）及其前体叠氮环醇在药物化学与生物活性分子合成中具有重要价值，尤其C₈-氨基环醇与C₈-叠氮环醇是相关合成研究中的关键骨架。因此，阐明该类底物的选择性取代规律，不仅有助于解释既有实验现象，也有助于理解柔性中环体系中构象因素对SN2反应的控制机制。

现有问题在于，前期实验已观察到二甲磺酸酯1在DMF中与叠氮化钠反应4天后，仅分离得到单叠氮产物2，收率为83%，而二叠氮产物3未被检测到。虽然这一结果清楚显示第一步取代具有显著动力学优势，但单取代终止的根本原因并未得到解释。特别是，构象效应、位阻限制与电子因素各自对活化屏障的贡献尚不清楚。基于这一背景，研究人员采用密度泛函理论（DFT）对该逐步SN2反应进行系统分析，旨在从分子层面解释该显著动力学选择性的来源。

研究人员的主要工作是：建立二甲磺酸酯1与N₃^?逐步发生两次SN2取代的反应路径，比较第一步与第二步反应的自由能垒、过渡态几何、畸变能与相互作用能差异，并结合非共价相互作用分析与动力学建模，解析反应停留在单叠氮阶段的原因。研究得到的结论十分明确：第一步取代经由TS1进行，活化自由能为31.7 kcal mol^?1，可以在实验温度下发生；第二步取代经由TS2进行，活化自由能升高至43.8 kcal mol^?1，尽管热力学上仍有利，但在动力学上几乎不可及。其本质原因并非外部溶剂或温度条件不足，而是底物内部构象预组织程度不同所导致的活化应变差异。该研究的重要意义在于，将单取代选择性的来源归结为柔性环体系中离去基预组织不足与额外构象重组代价，为理解中等环体系的反应性提供了清晰的结构—能量关系图景。

研究人员主要采用了以下关键技术方法：首先利用CREST程序在GFN2-xTB水平进行构象采样，筛选低能构象；随后在SMD/DMF溶剂模型下，采用M06-2X/6-311++G(d, p)水平进行几何优化与频率分析，并以IRC计算确认过渡态连接的反应物与产物极小值；热化学校正借助GoodVibes完成，并引入准谐振修正。之后，研究人员应用活化应变模型（ASM）分解活化能来源，利用Multiwfn与VMD进行非共价相互作用（NCI）/约化密度梯度（RDG）分析，并根据计算速率常数建立时间依赖的动力学模型。本文未涉及外部样本队列。

在研究结果部分，论文首先通过“逐步叠氮取代的相对自由能剖面”说明两步反应存在显著动力学分化。第一步中，N₃^?从环辛烷骨架C2位发生亲核进攻，经TS1形成单叠氮衍生物2，反应自由能与活化自由能分别为?14.0和31.7 kcal mol^?1。第二步则在形成单叠氮化合物2之后发生，但其最低能构象中剩余甲磺酸酯基（OMs）并不适于背面进攻，因此在成键前必须先进行构象重组，使离去基转变为适合SN2的拟轴向取向。正是这一额外重组代价，使TS2只能从更高能构象进入，最终导致43.8 kcal mol^?1的高活化能垒。计算结果与实验完全一致，即反应在单叠氮阶段终止，而不继续生成二取代产物。

接着，在“活化参数与反应性差异来源”部分，研究人员通过比较两步反应的活化自由能、畸变能、相互作用能、速率常数与理论半衰期，定量揭示第二步为何明显更慢。第一步到第二步的活化自由能增加了12.1 kcal mol^?1。ASM分析表明，这一变化主要来自更高的畸变能与更差的相互作用能。第二步过渡态需要骨架承受更大程度的几何变形，说明环辛烷框架必须为实现反应构型而付出更高结构代价。

在“NCI分析揭示过渡态中的稳定化与排斥效应”部分，论文指出TS1中的两个苯环采取包覆分子核心的空间布局，由此形成更明显的色散稳定作用，对应较有利的相互作用能?16.7 kcal mol^?1。而TS2中，由于第二步取代要求发生从准赤道向准轴向的重组，苯环被迫远离核心区域，导致稳定性色散作用减弱。与此同时，TS2还出现新引入叠氮基与相邻OBn取代基之间不利的gauche相互作用，以及两个叠氮基之间增强的静电排斥。分子紧凑性的降低还加剧了环辛烷骨架内部跨环氢原子相互作用。这些因素共同导致畸变能增加8.0 kcal mol^?1，从而显著抬高第二步活化能垒。相应地，计算半衰期显示第一步约为2.2天，与实验4天的反应时间尺度相符；第二步则高达60 900年，足以解释为何二取代产物在实验条件下无法形成。

在“过渡态几何分析”部分，研究人员进一步比较了两步反应中亲核体的进攻轨迹。TS1的N-C-O角为162.1°，C-N-N角为102.8°，更接近理想四面体sp³几何所允许的顺畅背面进攻，因此进入反应中心的路径更平缓、位阻更小。相比之下，TS2的N-C-O角降至154.9°，C-N-N角增至125.8°，显示叠氮离子采取更“正面”的进攻方向，明显偏离优选几何。这样的角度畸变迫使亲核体进入更拥挤的空间区域，是第二步动力学熄灭的重要原因之一。该部分结果从几何角度支持了ASM与NCI给出的能量学解释。

在“动力学模型与时间依赖浓度变化”部分，论文指出反向速率常数比正向速率常数低若干数量级，因此可将体系视为在研究条件下有效不可逆。基于此，研究人员建立了逐步反应动力学方程，并利用计算得到的速率常数求解浓度随时间变化的解析表达式。模拟结果显示，二甲磺酸酯1迅速消耗，单叠氮中间体2同步生成，并在模拟时间内接近完全转化；而二叠氮产物3在整个过程中始终维持可忽略水平。该结果进一步说明，反应体系并非处于“缓慢向双取代推进”的状态，而是由于第二步本征活化障碍过高，被实质性锁定在单取代阶段。

在“温度与溶剂效应的解释”部分，研究人员强调，尽管反应在378.15 K高温下进行，且DMF作为极性非质子溶剂通常有利于提升N₃^?的亲核性并促进SN2过程，但第二步仍然无法发生。这表明反应结果并不受限于外部反应介质，而主要由底物自身结构约束决定。尤其是环辛烷骨架的构象刚性与位阻限制，对第二次取代施加了显著能量惩罚。

讨论部分的核心在于：双重亲核取代中的反应性分化，本质上不是由“第二次成键本身更困难”单独造成，而是由第二步在成键前必须完成的离去基重新定位所主导。第一步反应中，被进攻的OMs离去基原本就位于有利的反应角位，能够以较小构象调整直接发生SN2取代；第二步则必须将OMs从head区域重新转移至corner位置，才能达到可反应排布，因此引入额外扭转畸变与更广泛的构象重组。扭转分析进一步支持这一认识：低能路径中的二面角调整更加局部且协同，而高能路径中的畸变则更分散、更不均匀。由此，论文将“构象预组织”确立为该柔性中环体系反应性的决定性因素，并指出这种由head-to-corner重组惩罚导致的反应性抑制，可能是中等环体系中的一般行为特征。

研究结论部分可译为：总之，本研究为顺式-1,2-二甲磺酸酯环辛烷1在双重亲核取代中表现出的反应性分化提供了详细的机理解释。计算能量剖面明确表明，第一步取代经TS1进行，活化能垒为31.7 kcal mol^?1，而第二步取代经TS2进行，其能垒显著升高至43.8 kcal mol^?1，因此构成整个转化的速率控制步骤。畸变—相互作用分析显示，这一增加主要由畸变能显著上升所驱动，说明环辛烷骨架为了到达第二过渡态必须经历更大的几何变形；同时，相互作用能也变得不利，这与NCI分析所揭示的位阻排斥增加和稳定化相互作用减弱相一致。与低能路径相比，TS2中的更高能垒可归因于OMs离去基缺乏预组织。在TS1路径中，亲核进攻发生于已位于反应角位的OMs基团，可在最小构象调整下直接发生SN2离去；而TS2则要求OMs在成键前从head区域迁移至corner位置，从而引入额外扭转畸变与构象重组。扭转分析进一步支持了这一解释：低能过渡态仅涉及更局部、更协同的二面角调整，而高能结构则表现为更分散且不均一的畸变模式。总体而言，这些结果表明，在该柔性环体系中，决定反应性的关键因素并非单纯的成键过程，而是构象预组织。这一行为对于存在head-to-corner重组惩罚的中等环体系可能具有普遍意义。