一、 研究作者、机构与发表信息
本研究由美国普林斯顿大学(Princeton University)的Eve Y. Xu, Casey B. Roos, Andrew J. Bendelsmith和Robert R. Knowles,以及犹他大学(University of Utah)的Jacob Werth和Matthew S. Sigman共同完成。通讯作者为普林斯顿大学的Robert R. Knowles教授和犹他大学的Matthew S. Sigman教授。该研究以题为“Noncovalent Stabilization of Radical Intermediates in the Enantioselective Hydroamination of Alkenes with Sulfonamides”的论文形式,发表于化学领域顶级期刊《Journal of the American Chemical Society》(简称J. Am. Chem. Soc.)2022年第144卷(第18948–18958页)。
二、 学术背景与研究目的
本研究属于不对称催化与自由基化学的交叉领域。不对称催化是现代合成化学的核心,旨在高选择性地合成手性分子(如药物分子)。其中,利用非共价相互作用来引导底物与手性催化剂的结合,是近年来在极性反应(涉及带电荷或极性中间体的反应)中实现高手性选择性的重要策略。然而,对于反应性极高、通常为电中性的自由基中间体,其与手性催化剂之间的非共价相互作用本质尚不清晰,相关研究非常有限。这制约了自由基不对称催化反应的发展,也限制了对生物酶催化、材料科学和超分子化学中类似自由基-分子相互作用的理解。
本研究旨在深入探究一个具体反应体系中的“自由基非共价相互作用”。该体系为研究团队在2020年开发的一种高对映选择性烯烃氢胺化反应(图1a),其通过质子耦合电子转移机制活化磺酰胺,产生中性的磺酰胺自由基中间体,随后与分子内的烯烃发生环化。初步研究表明,对映选择性的来源可能在于该中性自由基与手性磷酸催化剂(CPA)之间独特的非共价相互作用网络。特别是,催化剂分子中的三唑基团对于获得高选择性至关重要。因此,本研究的具体目标是通过详尽的实验、计算和数据分析,阐明该反应的对映选择性决定步骤,揭示磺酰胺自由基与手性磷酸之间稳定过渡态的具体非共价相互作用类型、结构和能量贡献,从而为自由基非共价相互作用的普适性理解提供分子层面的深刻见解。
三、 详细研究流程与方法
本研究是一个多方法联用的综合性机理研究,主要包含以下三个相互印证的研究阶段:
第一阶段:实验动力学研究——确定对映选择性决定步骤
研究对象: 以模型底物1(含有烯基和磺酰胺基团的特定分子)为核心,在标准反应条件下进行动力学分析。研究体系还包括手性磷酸催化剂P1、光催化剂PC和氢原子转移(HAT)试剂三苯甲基硫醇。
研究过程与实验方法: 1. 整体反应动力学监测: 采用¹H光核磁共振(¹H PhotonMR)光谱技术,在-20°C下实时监测反应进程(图2)。这是一种将LED光源集成到NMR谱仪中的先进技术,允许在光照引发反应的同时,原位、实时地通过NMR信号追踪反应物消耗和产物生成,特别适用于光催化反应的动力学研究。 2. 反应级数确定: 应用初速率法,通过改变各组分(底物1、硫醇、光催化剂PC、手性磷酸P1)的初始浓度,分别测量其对反应初始速率的影响(图3)。这通过分析低转化率(<10%)下产物浓度随时间变化的线性斜率实现。 3. 催化剂结合亲和力测定: 通过³¹P NMR滴定实验,研究P1与PC或底物1的结合。通过监测P1的³¹P化学位移随滴定物浓度的变化,利用Benesi-Hildebrand方程处理数据,计算结合常数K,以评估P1与各组分结合的相对强弱。 4. 非线性效应(NLE)研究: 使用不同对映体纯度的P1进行反应,观察产物对映选择性(ee值)与催化剂对映体纯度(ee值)之间的关系,以判断对映决定步骤中涉及的手性催化剂分子数。
第二阶段:计算化学研究——建立过渡态模型与相互作用分析
研究对象: 对映选择性决定步骤(C–N键形成)的微观过程。计算模型包括:磺酰胺自由基1a、模拟P1共轭酸的手性磷酸PA1(去除了用于增加溶解度的烷基侧链以降低计算成本)、以及环化后生成的碳中心自由基1c。
研究过程与计算方法: 1. 反应路径构建与能量计算: 使用密度泛函理论(DFT)分别计算有PA1存在和无PA1存在时的C–N键形成反应坐标(图4)。首先使用B3LYP/6-31G(d)理论水平进行几何结构优化和构象搜索,然后使用M06/6-311+G(d,p)/CPCM(甲苯溶剂)理论水平进行单点能校正。这种组合旨在更准确地描述非共价相互作用,特别是色散作用。 2. 过渡态(TS)结构分析与非共价相互作用识别: 对导致主产物(major)和次产物(minor)的两个非对映异构过渡态(1b–PA1maj 和 1b–PA1min)进行详细几何结构和电子结构分析(图5)。使用NCIplot程序可视化非共价相互作用区域。通过测量关键原子间距离、观察自旋密度分布、分析分子轨道和静电势,识别并归类稳定过渡态的具体相互作用类型(如氢键、C–H···π相互作用、孤对电子-π相互作用等)。 3. 模型普适性检验: 将建立的过渡态相互作用模型应用于另外两个结构不同的底物:高选择性的噻唑底物3和选择性较差的偕二甲基底物4(图6)。计算它们的过渡态结构和能量差,与实验选择性对比,验证模型的解释和预测能力。
第三阶段:统计建模研究——量化结构与选择性关系
研究对象: 一个包含7种不同磺酰胺底物和10种不同手性磷酸催化剂的组合反应数据集,共计40个反应,其选择性范围从12%到96% ee。
研究过程与数据分析方法: 1. 描述符库构建: 对底物和催化剂进行一系列计算,提取多种描述符。包括:自然键轨道(NBO)电荷、Sterimol立体参数、全局参数(如HOMO/LUMO能级、极化率)。尤为重要的是,构建了“探针分子”来计算关键的“非共价相互作用”参数(图7b)。例如,模拟噻唑底物3与催化剂三唑环之间孤对电子-π相互作用,使用对称性匹配微扰理论(SAPT)在SSAPT0/jun-cc-pvdz水平上计算该相互作用的平衡距离和总相互作用能,以此作为描述符。 2. 多元线性回归分析: 以实验测得的对映选择性(ee值)为因变量,以上述计算得到的描述符为自变量,进行多元线性回归建模(图7a)。采用前向逐步算法筛选出与选择性最相关的描述符组合。 3. 模型验证: 使用留一法交叉验证(LOO)和k折交叉验证(k=4)等统计方法评估模型的稳健性和预测能力(得到预测R² = 0.83)。
四、 主要研究结果及其逻辑关系
1. 动力学结果确立了C–N键形成是对映选择性和速率决定步骤: * 反应对底物1呈一级动力学关系,但对硫醇浓度呈零级关系(图3a)。这表明氢原子转移(HAT)步骤不是决速步,从而排除了HAT步骤决定对映选择性的可能性。 * 反应速率对手性磷酸P1的浓度呈现饱和动力学特征(图3c),并与P1–PC复合物的浓度相关。NMR滴定显示P1与PC的结合常数远高于与底物1的结合,说明反应中P1–PC复合物可能是质子耦合电子转移(PCET)过程的活性物种。 * 非线性效应研究表明,产物ee值与催化剂ee值呈线性关系,表明在对映决定步骤中仅涉及一个手性磷酸分子。 * 综合以上,结合PCET步骤可逆且快速的特点,研究得出结论:C–N键形成步骤是不可逆的、速率决定的,同时也是对映选择性决定的步骤。 这为后续聚焦于该步骤的计算研究提供了坚实的实验基础。
2. 计算研究揭示了强氢键和次级非共价相互作用网络是选择性的来源: * 强氢键络合阻止背景反应: 计算发现,PCET产生的磺酰胺自由基1a与手性磷酸PA1能形成极强的氢键络合物1a–PA1,结合自由能达-12.4 kcal/mol。这个强相互作用使得自由基在后续反应中始终与手性催化剂结合,有效抑制了无催化剂参与的非手性背景反应路径(其能垒仅1.94 kcal/mol)。 * 过渡态能量差与实验吻合: 计算得到导致主、次产物的两个过渡态之间的自由能差ΔΔG‡为1.7 kcal/mol(对应94% ee),与实验测得的1.6 kcal/mol(92% ee)高度吻合,验证了计算模型的可靠性。 * 关键非共价相互作用的可视化与比较(图5): * 主产物过渡态(TS-maj): 包含一个双重氢键网络:一是CPA的酸性质子与磺酰胺氮原子之间的强主氢键(~1.54 Å);二是CPA的磷酸氧原子与新生碳中心自由基邻近的两个脂肪族C–H键之间的较弱次级氢键。此外,存在关键的次级非共价相互作用:底物芳环的两个C–H键与CPA三唑环的π面形成C–H···π作用;底物的异丙基C–H键与CPA一个苯环的π面也形成C–H···π作用。这些相互作用共同稳定了主产物过渡态。 * 次产物过渡态(TS-min): 同样存在双重氢键网络。然而,其CPA的三唑环远离底物,缺乏有效的C–H···π稳定作用。相反,CPA的三唑和苯环上的C–H键转而与底物的磺酰氧原子形成一系列弱氢键。 * 模型普适性验证(图6): * 对于高选择性的噻唑底物3,在主产物过渡态中,CPA的三唑环与底物噻唑的硫原子之间形成了孤对电子-π相互作用,取代了模型底物的C–H···π作用,计算结果与实验选择性一致。 * 对于低选择性的底物4,在其次产物过渡态中,一个甲基的C–H键意外地与CPA的三唑环产生了较强的C–H···π相互作用,从而稳定了次产物过渡态,降低了选择性。计算成功预测了这一趋势。 * 逻辑递进: 从确定关键步骤,到建立该步骤的精确过渡态模型,再到用不同底物验证模型的解释力,计算研究逐步构建并证实了一个以强氢键锚定、以三唑环参与的次级非共价相互作用为关键区分因素的立体控制分子模型。
3. 统计建模量化了催化剂和底物结构对选择性的影响: * 回归分析得到了一个具有良好预测能力的三参数模型(R² = 0.85)。三个关键描述符及其意义如下: 1. SAPT相互作用能(催化剂-探针): 负系数表明,催化剂与模拟关键相互作用的探针(如二甲基硫醚)之间的总相互作用能越低(即结合越强),反应选择性越高。这直接量化并强调了催化剂参与特定非共价相互作用能力的重要性。 2. CPA三唑芳环上C5原子的NBO电荷: 正系数表明,吸电子取代基(使该点更正电性)有利于高选择性。这反映了3,5-二取代催化剂通过调整催化口袋的电子性质和形状来优化相互作用。 3. 底物产物的2D拓扑描述符(PW5): 该参数成功捕捉到底物4因偕二甲基取代导致的立体位阻负面影响,从而选择性降低。 * 逻辑关系: 统计建模独立地从大量结构-活性数据中提取出了关键因子,其结果与基于单个过渡态深度分析的计算结论相互印证。特别是,SAPT相互作用能描述符直接证实了三唑环参与的相互作用在决定全局选择性趋势中的核心作用。
五、 研究结论与价值
本研究通过多学科交叉的方法,首次对一个涉及中性氮中心自由基的高对映选择性反应进行了深入、定量的机理剖析,并得出以下核心结论:
科学价值: 本研究超越了简单提出“非共价相互作用可能重要”的假设,而是通过精确的实验设计和先进的计算工具,提供了关于自由基-分子非共价相互作用的强度、结构和能量贡献的定量分子图景。它深化了学术界对高反应性开壳层中间体如何通过与手性环境进行弱相互作用来实现高度立体控制的理解,为理性设计更多基于非共价相互作用的自由基不对称反应奠定了理论基础。
应用价值: 研究所揭示的“强氢键锚定+次级芳环相互作用微调”的催化模式,可以作为未来设计新型手性催化剂(特别是用于自由基反应)的指导原则。对三唑基团关键作用的深刻理解,将启发化学家在催化剂骨脑中合理引入特定芳杂环,以构建高效的立体控制口袋。
六、 研究亮点
七、 其他有价值内容
本研究还展示了现代物理有机化学研究中一些精妙的实验和计算技术细节的价值: * ¹H PhotonMR技术的应用: 为光催化反应的原位、实时动力学研究提供了强大工具,避免了传统淬灭取样法可能带来的误差,特别适合研究快反应和涉及不稳定中间体的反应。 * 非线性效应(NLE)实验的巧妙运用: 用一个简单而经典的实验,高效地确定了对映决定步骤中催化剂的聚合状态(单体作用),从而极大地简化了后续复杂的计算模型构建,体现了实验设计的重要性。 * “探针分子+SAPT计算”的描述符构建策略: 在统计建模中,为了量化难以直接计算的复杂非共价相互作用,研究团队创造性地采用构建简化模型探针并用SAPT理论计算其相互作用能的方法,成功地将关键的微观相互作用信息转化为可用于宏观回归分析的定量描述符,这是一个非常具有启发性的方法学创新。 2. 这篇文献解决了什么问题? 该研究旨在探究自由基参与的烯烃与磺酰胺的不对称氢胺化反应中,对映选择性产生的根本原因。其核心观点是:在该反应中,高水平的手性控制源自于反应中间体——中性的磺酰胺基氮自由基——与手性磷酸催化剂之间形成的一组吸引性非共价相互作用,特别是N···H强氢键以及涉及芳基(如三唑环)的次级非共价相互作用网络。这些相互作用稳定了生成优势对映体的过渡态。
实验动力学研究:使用¹H光NMR光谱进行原位监测,确定了C-N键形成是反应的决速步和对映选择性决定步骤。 密度泛函理论计算:对关键C-N键形成步骤进行了详细分析,构建并比较了生成主、副产物对映体的过渡态结构,首次通过计算揭示了中性磺酰胺基氮自由基作为强氢键受体的能力,并可视化识别了稳定优势过渡态的特定非共价相互作用(如C-H···π相互作用)。 数据科学的统计建模:通过多变量线性回归分析,将实验测得的大量对映选择性数据与从催化剂和底物计算的多种描述符(如SAPT计算的相互作用能、NBO电荷、拓扑描述符)相关联,从统计学上确认了催化剂三唑部分参与的相互作用对选择性具有关键贡献,从而补充并验证了过渡态模型。 4. 作者承认的局限性是什么? 文章在计算部分明确指出:
B3LYP泛函在描述色散作用和结合能方面存在不足(会低估结合能)。 尽管通过使用M06类泛函进行单点能校正来弥补这一缺陷,但作者也尝试了包含经验性色散校正的其他泛函,发现它们与实验值存在显著差异,可能是因为过度稳定了计算中的复合物结构。 这暗示了对涉及自由基的非共价相互作用的精确理论计算仍存在挑战。更广泛地说,文章引言部分指出,尽管非共价相互作用在极性反应的不对称催化中已被广泛研究,但人们对涉及自由基中间体的非共价相互作用的结构、性质及其在不对称催化中的作用机制知之甚少,而这正是本研究和该领域需要克服的障碍。