儿茶酚取代调控5-芳亚甲基-噻唑烷-2,4-二酮类化合物的抗氧化效力：合成、构效关系、密度泛函理论、分子对接与分子动力学研究

《Results in Chemistry》：Catechol substitution governs antioxidant potency in 5-arylidene-thiazolidine-2,4-diones: Synthesis, structure-activity relationships, DFT, docking, and molecular dynamics

　　

摘要翻译：
杂环骨架如噻唑烷-2,4-二酮(TZD)及其衍生物是满足多功能治疗药物需求的优良候选化合物。本研究通过以吗啉为催化剂的简单两步法合成了一系列TZD衍生物，并通过红外光谱(IR)、紫外-可见光谱(UV–vis)、¹H核磁共振(NMR)和¹³C NMR光谱进行了表征。在抗氧化测定中，5-(3,4-二羟基苯亚甲基)-噻唑烷-2,4-二酮(TZD-1)在DPPH(0.3283 ± 0.0046 ppm)和ABTS(0.3232 ± 0.0159 ppm)测定中均表现出最低的半数抑制浓度(IC₅₀)值，且铁离子还原抗氧化能力(FRAP)值最高(4.5424 ± 0.3558 Fe²⁺ mM eq/g)。总酚和总黄酮含量与抗氧化能力之间的皮尔逊相关性支持了构效关系。在初步MTT筛选中，TZD-1降低了HeLa细胞活力，IC₅₀为885 μg/mL(约3.7 mM)，而正常Vero细胞在高达4000 μg/mL浓度下基本不受影响(选择性指数>4.5)。密度泛函理论(DFT)计算证实了反应性与较低的最高占据分子轨道-最低未占分子轨道(HOMO-LUMO)能隙及化学软度之间的相关性。针对细胞周期蛋白依赖性激酶-2(CDK2)的分子对接(最高达?8.19 kcal/mol)揭示了通过氢键和π–π相互作用在催化相关残基处的有利结合。所有五种衍生物均满足利平斯基五规则(Lipinski's Rule of Five)，且TZD-1是SwissADME预测的唯一不抑制五种细胞色素P450(CYP450)亚型中任何一种的衍生物。MM/GBSA计算显示TZD-1–CDK2复合物的结合自由能为?34.54 ± 11.23 kcal/mol，分子动力学模拟(100 ns)也通过在整个平衡轨迹中一致的氢键相互作用确认了其稳定性。总体而言，TZD-1代表了用于未来优化的结构上有效的生物活性骨架。

论文解读

**一、研究背景与目的**

杂环化合物因其独特的电子性质和多样的生物活性，在药物化学与合成化学中占据核心地位。自然界中超过85%的生物活性化合物至少含有一个杂环。噻唑烷-2,4-二酮(TZD)作为一种重要的含硫杂环骨架，其结构特征在于2位和4位的羰基作为电子受体，而NH/S基团作为电子供体，这种结构使其能够与DNA、受体和酶等生物靶标有效相互作用。TZD衍生物已被证实具有广泛的药理活性，包括抗菌、抗疟、抗炎、抗糖尿病、抗癌、抗氧化和抗结核等。尽管已有多种催化剂被用于合成TZD衍生物，但许多现有方法依赖对水分敏感的催化剂、专用设备或多步纯化过程。

本研究旨在通过简单的两步合成法，以吗啉为有机催化剂制备一系列5-芳亚甲基-TZD衍生物，并从多个维度系统评估其生物活性潜力。研究的核心科学问题在于：不同取代基，特别是儿茶酚基团，如何影响TZD衍生物的抗氧化效力、细胞毒性和电子结构特征。研究创新性地将多测定抗氧化评估、总酚/总黄酮含量定量、皮尔逊相关性构效分析、癌细胞与正常细胞对比毒性评估，以及密度泛函理论(DFT)和分子对接等计算分析方法相结合，为这类化合物提供了前所未有的综合评价。

**二、主要技术方法**

研究人员合成了五种TZD衍生物(TZD-1至TZD-5)，所有化合物通过红外光谱(IR)、紫外-可见光谱(UV–vis)、¹H NMR和¹³C NMR光谱进行表征。抗氧化活性通过DPPH、ABTS和FRAP三种互补方法评估，同时测定了总酚含量(TPC)和总黄酮含量(TFC)。细胞毒性通过MTT法在HeLa(人宫颈癌)和Vero(非洲绿猴肾上皮)细胞系中评估。计算研究包括：在B3LYP/6-31+G(d,p)水平进行的DFT计算；针对CDK2(PDB: 7KJS)的AutoDock Vina分子对接；GROMACS 2023.2平台上100 ns分子动力学(MD)模拟及MM/GBSA结合自由能计算；以及SwissADME药代动力学预测。

**三、研究结果**

**3.1 化学合成与表征**

通过以吗啉为催化剂的Knoevenagel缩合反应，以57-75%的产率成功合成了五种TZD衍生物。构效关系(SAR)原则指导了衍生化策略：TZD-1中的儿茶酚基团因其已知的自由基清除药效团而被选择；给电子基团(二甲氨基、甲氧基)用于评估提高HOMO能量对生物活性的影响；邻位和对位甲基取代用于区分位阻效应与电子效应对共轭和结合亲和力的影响。

**3.2 DPPH和ABTS抗氧化活性**

TZD-1在ABTS测定(IC₅₀ = 0.3232 ± 0.0159 ppm)和DPPH测定(IC₅₀ = 0.3283 ± 0.0046 ppm)中均表现出最低的IC₅₀值，甚至低于标准抗氧化剂抗坏血酸(ABTS: 4.8734 ± 0.8436 ppm；DPPH: 10.3088 ± 0.3393 ppm)。三种测定的活性排序一致为TZD-1 > TZD-4 > TZD-3 > TZD > TZD-5 > TZD-2。TZD-2活性最低可能是因为邻位取代降低了分子平面性，削弱了自由基离域能力。TZD-5尽管具有给电子的二甲氨基，但因缺乏可交换质子，只能通过单电子转移(SET)机制发挥作用，动力学上不利，故活性较低。

**3.3 FRAP铁离子还原抗氧化能力**

TZD-1表现出最高的FRAP值(4.5424 ± 0.3558 Fe²⁺ mM eq/g)，与其他测定结果一致。FRAP测定反映的是热力学还原能力，这一结果表明TZD-1的抗氧化活性是热力学驱动的，而非纯动力学假象，在慢性氧化应激(如糖尿病肾病)条件下具有治疗意义。

**3.4-3.5 总酚和总黄酮含量**

TZD-1在TPC测定中表现出最高值(225.3795 ± 0.4307 mg/g GAE)，直接反映了其结构中活性儿茶酚单元的存在。在TFC测定中，TZD-1同样最高(33.6022 ± 0.4045 mg/g QE)，这与AlCl₃比色法检测羰基邻位OH基团和与黄酮骨架结构类似的共轭芳香体系的能力一致。

**3.6 细胞毒性评估**

MTT测定显示，TZD-1显著降低HeLa细胞活力(存活率20-30%)，而对正常Vero细胞基本无影响(存活率>95%)。剂量-反应曲线外推得到TZD-1对HeLa细胞的IC₅₀约为885 μg/mL(约3.7 mM)，对Vero细胞的IC₅₀超过4000 μg/mL，选择性指数(SI)>4.5，表明其对癌细胞具有选择性毒性。尽管与标准抗癌药物阿霉素(IC₅₀ = 0.983 μg/mL)相比活性较弱，但该结果将TZD-1定义为需要进一步验证的初步细胞毒性候选化合物。

**3.7 统计分析**

皮尔逊相关性分析显示，TPC与FRAP(r = 0.963, p < 0.01)和TFC与FRAP(r = 0.960, p < 0.01)之间存在高度显著的正相关，表明酚类和黄酮类含量越高，抗氧化能力越强。TPC/TFC与IC₅₀值之间的相关性较弱且不显著，这可能是因为酚类单元的结构环境(而非总数)控制着自由基清除动力学，TZD-1中的邻二羟基儿茶酚是结构优化的氢原子转移药效团。

**3.8 DFT和MESP分析**

DFT计算显示，TZD-5具有最高的HOMO能量(?5.39 eV)和最低的能隙(3.34 eV)，表明其具有较高的电子反应性和电荷转移能力。化学软度计算中TZD-5最高(0.599 eV^?1)，对应更高的分子极化率和分子间相互作用灵活性。分子静电势(MESP)分析揭示了沿分子骨架的电荷分离，氧原子为主要的电荷富集中心，而硫原子带高度正电荷(+0.55至+0.61 e)，这种极化电荷分布与CDK2对接中观察到的氢键相互作用一致。

**3.9-3.10 分子对接与CDK2结合**

重对接验证实验的均方根偏差(RMSD)为0.3578 ?，远低于2.0 ?的接受阈值，确认了对接方案的可靠性。对接结果显示，TZD-1的得分最有利(?8.188 kcal/mol)，明显优于TZD-2(?6.996 kcal/mol)和TZD-3至TZD-5(?5.518至?5.761 kcal/mol)。TZD-1与CDK2的ATP结合口袋中的保守残基Lys33和Glu51形成氢键，并与Phe80形成π–π堆积相互作用，这归因于其儿茶酚部分作为氢键供体以及富电子芳香环的堆叠能力。

**3.11 分子动力学模拟**

100 ns MD模拟显示，TZD-1–CDK2复合物在约50 ns后达到稳定平衡。在平衡的50-100 ns区间内，Cα RMSD稳定在0.251 ± 0.010 nm，配体RMSD为0.143 ± 0.030 nm，表明TZD-1保持了稳定的结合取向。回转半径(Rg)平均为2.584 ± 0.008 nm，表明复合物保持紧凑结构。溶剂可及表面积(SASA)为254.756 ± 3.272 nm2，表明溶剂暴露受控。氢键分析显示TZD-1在整个稳定阶段维持3-5个氢键。MM/GBSA计算得到总结合自由能为?34.54 ± 11.23 kcal/mol，主要由范德华力(?22.20 ± 6.88 kcal/mol)和静电相互作用(?43.86 ± 15.20 kcal/mol)驱动。每残基分解分析确定了ASP145、LEU83、LEU134等主要稳定残基。

**3.12 ADME分析**

所有五种化合物均满足Lipinski五规则，零违例，分子量219-265 Da，log P在0.95-2.07之间，生物利用度评分0.55。所有化合物均预测具有高胃肠道吸收且非P-糖蛋白底物。CYP抑制是主要区分因素：仅TZD-1被预测不抑制五种CYP450亚型中的任何一种，而TZD-2至TZD-5被标记为CYP1A2和CYP2C19的潜在抑制剂。综合来看，TZD-1因无CYP抑制、良好的溶解性和简单合成(合成可及性SA = 2.83)而成为该系列中ADME性质最优越的化合物。

**四、讨论与结论**

研究表明，儿茶酚取代基是TZD衍生物抗氧化活性的关键决定因素，TZD-1中的邻二羟基结构通过氢原子转移机制实现高效自由基清除，同时赋予其独特的CDK2结合模式和选择性细胞毒性。研究证实了电子结构参数(能隙、化学软度)与实验生物活性之间的相关性，为理性设计提供了理论基础。分子动力学模拟和MM/GBSA计算进一步验证了TZD-1–CDK2复合物的动态稳定性。

综合来看，TZD-1(5-(3,4-二羟基苯亚甲基)-噻唑烷-2,4-二酮)代表了具有良好成药性的生物活性骨架，值得进一步的修饰与开发。然而，该研究存在一些明确局限性：生物评价仅限于体外实验和计算机模拟，未进行体内研究；细胞毒性仅在单一癌细胞系(HeLa)中通过MTT法评估，缺乏凋亡分析、细胞周期分析等机制研究；MD模拟仅针对TZD-1–CDK2单一复合物进行了单次100 ns轨迹；化合物为已知结构且衍生物库规模有限。这些局限性为未来研究指明了方向，包括多细胞系评价、机制细胞生物学研究、体内实验以及更广泛的MD模拟。