近日，深圳湾实验室/北京大学吴云东/张欣豪课题组与德克萨斯大学圣安东尼奥分校Michael P Doyle课题组合作，通过高通量DFT计算实现了自由基诱导的卤甲烷原子转移反应的精准预测。相关研究成果以“Precise Introduction of the -CHnX3-n(X=F, Cl, Br, I) Moiety to Target Molecules by a Radical Strategy: A Theoretical and Experimental Study”为题发表在Journal of the American Chemical Society 期刊上。
在生物医药领域，通过改变化合物中卤素原子的种类和数量是调控生物活性分子物化性质和药理性能的有效方法，因此科研工作者们致力于发展将卤甲基基团引入目标分子的合成方法。
卤甲烷是结构最简单的卤代烃，即CH₄分子中的一个或多个H原子被卤素原子（F, Cl, Br, I）所取代，易于制备，是最常见的化工原料之一。碳自由基的形成来源于分子共价键的断裂，种类丰富，反应活性高。在化学反应中，碳自由基可以诱导卤甲烷发生卤素原子转移反应(XAT)或H原子转移反应(HAT)（图1），从而产生特定的卤甲基自由基并与烯烃偶联得到目标产物（图2）。这类合成策略原料丰富，绿色经济，有广阔的探索和发展空间。发展这类合成策略的瓶颈在于如何实现卤甲烷和碳自由基的的精准匹配，从而控制卤素原子转移反应或H原子转移反应的选择性。
图1
图2
随着大数据时代的来临和数据库技术的不断发展，数据驱动策略应运而生并为科学研究带来了深刻的范式转变，推动着科学家们挖掘复杂体系背后隐藏的更加有用的信息。借助新兴的数据驱动策略，研究人员采用高通量DFT方法计算了共912个碳自由基引发的卤甲烷的H原子转移和卤素原子转移反应。这些模型反应以全部种类的卤甲烷为底物，近40个带有不同性质官能团的碳自由基作为原子转移引发剂。高通量的DFT计算一共得到了2038个能量数据，包括每个模型反应的过渡态能垒（ΔG‡）和反应自由能（ΔG），以及全部卤甲烷和产物中关键C-H键、C-X键的键能（BDE）。此外，由于碳自由基电离势（IP）反映了碳自由基的亲核性质，研究人员计算了所有作为原子转移引发剂的碳自由基电离势（IP），以期探究H原子转移和卤素原子转移反应活性和选择性的影响因素（图3）。
图3
研究人员采用数学模型对计算得到的ΔG‡，ΔG，IP和BDE等参数进行了处理，得到了碳自由基引发卤甲烷原子转移反应的反应活性和选择性规律。原子转移反应的活性由原子转移反应的过渡态能垒（ΔG‡H和ΔG‡X）决定，线性关系分析表明过渡态能垒（ΔG‡）与反应自由能（ΔG）呈紧密正相关。原子转移反应的选择性由反应过渡态能垒之差（ΔG‡H-ΔG‡X）决定：ΔG‡H-ΔG‡X>0，发生卤素原子转移反应；ΔG‡H-ΔG‡X<0，发生H原子转移反应。多元回归分析表明（ΔG‡H-ΔG‡X）由反应能（ΔGH、ΔGX）和IP三参数共同描述最为准确（图4），意味着原子转移选择性由反应热力学和自由基的亲核性共同影响。值得一提的是，（ΔG‡H-ΔG‡X）与单参数的碳自由基电离势（IP）间线性关系良好，因此可以借助查阅实验手册中碳自由基电离势，快速判断反应能垒差值是否符合原子转移能量标准。
图4
研究人员依据理论上发生原子转移反应的能量判据，预测出能引发全部卤甲烷发生特定H原子转移或卤素原子转移的碳自由基。以CH₂Cl₂的Cl原子转移反应为例，能量判据需同时满足ΔG‡H-ΔG‡Cl>0和ΔGCl<0。因此，能引发Cl原子转移的碳自由基主要为α-氨基烷基自由基、α-烷氧基自由基（图5）。
图5
基于该预测，实验人员在α-氨基碳自由基诱导下引发了CH₂Br₂和CH₂I₂中Br和I原子转移，将-CH₂Br基团和-CH₂I基团引入目标分子中（图6）。
图6
该研究工作表明数据驱动策略在挖掘隐藏的化学的反应规律方面有潜在的应用，通过查阅已知的实验数据库可以精准匹配引发卤甲烷原子转移的碳自由基。
深圳湾实验室/北京大学吴云东教授，张欣豪副教授，德克萨斯大学圣安东尼奥分校Michael P Doyle教授为本文共同通讯作者，深圳湾实验室博士后（原北京大学化学学院博士毕业生）刘思琪为本文第一作者。
文章来源 | 刘思琪
编辑 | 鲍 啦
