新型吡唑类配体对EGFR和HER2的靶标依赖性结合与动态稳定化：整合分子对接、分子动力学和DFT研究

《Letters in Drug Design & Discovery》：Target-dependent binding and dynamic stabilization of novel pyrazole ligands toward EGFR and HER2: An integrated docking, molecular dynamics, and DFT investigation

　　

背景：表皮生长因子受体（EGFR）和人表皮生长因子受体2（HER2）是癌症中的重要治疗靶点。理解配体结构与电子性质如何与蛋白质–配体相互作用相关，可能有助于发现具有前景的激酶靶向化合物。本研究考察两种吡唑基配体L1和L2对EGFR和HER2的靶标依赖性结合与动态稳定化行为。方法：采用整合计算方法，结合吸收、分布、代谢、排泄和毒性（ADMET）预测、分子对接、100 ns分子动力学（MD）模拟、分子力学/广义Born表面积（MM/GBSA）计算、密度泛函理论（DFT）、含时密度泛函理论（TD-DFT）和自然键轨道（NBO）分析。结果：分子对接显示，L1对EGFR和HER2的预测对接分数比L2更有利（EGFR：?8.5 对 ?8.3 kcal/mol；HER2：?8.3 对 ?7.3 kcal/mol）。然而，MD轨迹揭示了不同的受体依赖性行为。L1–EGFR表现出相对较低的结构波动和更持久的蛋白质–配体相互作用，而L2–HER2表现出更一致的构象轮廓和更高的分子间接触保守性。MM/GBSA计算进一步支持在所应用计算框架内存在有利相互作用，尽管某些复合物的能量排序与对接不同。ADMET预测表明两种配体均满足Lipinski五规则。DFT、TD-DFT和NBO分析揭示了电子分布、轨道离域和分子内电荷转移特征的差异，为解释其分子性质提供了补充信息。结论：综合计算结果表明，L1和L2对EGFR和HER2表现出不同的结合和动态稳定化行为。这些结果为进一步研究吡唑基激酶配体提供了基础。然而，这些发现具有轨迹特异性，需要实验结合和激酶抑制测定以及独立MD重复，以确立所观察差异的可重复性和生物学相关性。

论文解读：癌症仍是全球主要死亡原因之一。受体酪氨酸激酶（RTKs，receptor tyrosine kinases）调控细胞增殖、分化、迁移和存活，其中表皮生长因子受体（EGFR，epidermal growth factor receptor）和人表皮生长因子受体2（HER2，human epidermal growth factor receptor 2）在乳腺、肺、结直肠和胃癌等恶性肿瘤中常过表达或失调，是重要抗癌靶点。Erlotinib和Lapatinib等临床激酶抑制剂虽有治疗获益，但存在获得性耐药、毒性及对近缘激酶靶标区分有限等问题；EGFR与HER2结构高度相似，使设计具有差异化靶标结合性质的小分子面临挑战。吡唑衍生物因氮-rich芳香框架、可参与氢键、π–π堆积和静电相互作用，以及可模拟ATP结合位点识别所需药效团，成为激酶靶向研究热点。然而，静态分子对接难以充分捕捉长时间稳定性和受体依赖构象适应。为此，Younes Larbi、Oussama Khaoua、Souhila Laib对两个新合成并结构表征的吡唑配体L1和L2开展研究，重点考察其与EGFR和HER2的靶标依赖性结合与动态稳定化行为。得到结论：L1在EGFR复合物中表现出更持久稳定化，L2在HER2结合口袋中表现出增强的动态稳定化；电子结构分析提供补充解释，但受体特异结合与抑制未实验测定，结果应视为计算预测而非实验确认的选择性或抑制。该研究为吡唑类激酶配体进一步研究提供基础，并强调在激酶靶向配体计算研究中考虑受体动力学的重要性。

主要关键技术方法（不超过250字）：研究人员采用整合计算策略，未涉及人类或动物样本队列；研究对象为已合成并结构表征的L1、L2及Protein Data Bank中的EGFR（PDB ID = 4HJO）和HER2（PDB ID = 3PP0）。用SwissADME和Protox II预测吸收、分布、代谢、排泄和毒性（ADMET）与药物相似性；用AutoDock Vina对接，PyMOL验证；用Desmond/Schr?dinger Maestro进行100 ns MD，OPLS4力场、TIP4P水模型、NPT系综300 K和1 atm；用Prime模块进行分子力学/广义Born表面积（MM/GBSA）；用Gaussian 16在B3LYP/6–31++G(d,p)水平、DMSO中C-PCM模型优化，并进行密度泛函理论（DFT）、含时密度泛函理论（TD-DFT）、自然键轨道（NBO）、分子静电势（MEP）和全局反应描述符分析。

研究结果：
ADMET性质：通过SwissADME和Protox II预测，L1和L2分子量均低于500 g/mol，Lipinski五规则零违例，多种药物相似性过滤器通过；L1溶解度更优，L2亲脂性更高；两者高胃肠吸收、非P-糖蛋白（P-gp）底物、非CYP3A4抑制剂；毒性分类均为4类，L2预测半数致死量（LD₅₀）较高；两者预测非肝毒、非免疫毒、非致突变、非细胞毒，但致癌性预测需谨慎，需实验验证。
Docking simulation：通过重新对接验证，Erlotinib–EGFR和03Q–HER2的均方根偏差（RMSD，root mean square deviation）分别为0.631 ?和0.824 ?。对接显示Erlotinib对EGFR分数?11.3 kcal/mol，L1?8.5，L2?8.3；Lapatinib对HER2?10.6 kcal/mol，L1?8.3，L2?7.3。L1对两靶标预测分数均优于L2。L1–EGFR形成Lys721、Thr830、Asp831氢键及疏水接触；L2–EGFR形成Asn818、Asp831氢键及π–cation、π–sulfur作用；L1–HER2形成Arg784、Thr862、Asp863氢键及疏水接触；L2–HER2形成Leu726、Cys805氢键和Asp808静电作用。对接仅表明接触类型和分布差异，氢键数不能单独解释结合强弱。
MD simulation：100 ns MD显示，L1–EGFR蛋白骨架RMSD约1.5–2.0 ?波动，配体RMSD较一致，残基均方根涨落（RMSF，root mean square fluctuation）低，Asp831等氢键和Val702、Ala719、Leu764疏水接触持久，配体回转半径（rGyr，radius of gyration）、溶剂可及表面积（SASA，solvent-accessible surface area）、极性表面积（PSA，polar surface area）变化有限；L2–EGFR偏差更大、波动尖峰、140–170区残基活动更高、水介导作用更多，配体构象与表面性质变化更大。EGFR–L1与Erlotinib–EGFR行为更接近。HER2体系中，L1–HER2配体RMSD波动明显、结合方向调整较大；L2–HER2在约20–30 ns后达到较一致RMSD并保持较窄波动，残基涨落较低，疏水、反复氢键及Asp残基持久水介导接触更一致，配体性质描述符较稳定；L2–HER2与Lapatinib–HER2行为较接近。首末帧比较显示，EGFR–L1保留Asp831、Thr830、Val702、Ala719等作用，EGFR–L2接触重排更明显；HER2–L1结合模式重排更明显，HER2–L2相互作用网络更保守。结论：单条100 ns轨迹中，L1在EGFR中更一致，L2在HER2中更一致，但轨迹特异性，需独立重复。
Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) Analysis：MM/GBSA显示所有复合物ΔG_bind为负，范围?34.87至?64.63 kcal mol^?1；范德华贡献主导，静电提供额外稳定，ΔG_GB正贡献不利，ΔG_SA负贡献有利。EGFR–L1 ΔG_bind?44.01，EGFR–L2?34.87，EGFR–Erlotinib?39.66；HER2–L1?38.67，HER2–L2?40.01，HER2–Lapatinib?64.63 kcal mol^?1。L1对EGFR计算结合能更有利，L2对HER2略更有利；MM/GBSA与对接排序不完全一致，因方法框架不同，应视为互补。
Electronic origin of dynamic kinase stabilization：优化几何显示L1较紧凑，给体中心主要位于吡唑氮和酚氧，电子分布较局域，可能有助于解释EGFR–L1持久相互作用；L2更延展，含两个吡唑环、中心叔胺和酚基，杂原子更多、共轭框架更大，轨道离域、极化率和构象适应性更大，可能有助于解释HER2–L2动态稳定，但不建立因果机制。
Frontier Molecular Orbital (FMO) Analysis：前沿分子轨道（FMO，frontier molecular orbital）分析显示L1的最高占据分子轨道（HOMO）主要位于吡唑和苯环，最低未占分子轨道（LUMO）主要位于酚基；L2的HOMO广泛离域于两个吡唑环和中心胺桥。L1与L2的HOMO–LUMO能隙相近（5.30867 eV对5.34513 eV），不能仅凭能隙区分反应性；FMO分布差异提供补充电子信息。
Global Reactivity Descriptors Analysis：全局反应描述符表明L1的E_HOMO较高、能隙略小、硬度较低、软度较高，具有相对更高的电子给予倾向和方向性响应；L2的E_LUMO较低、电负性和亲电指数较高、化学势更负，显示较强电子接受倾向和电子重分布能力。
Molecular electrostatic potential (MESP) analysis：MEP显示L1负电势主要集中在吡唑环，酚氧附近较不延展，整体静电分布较局域；L2负电势更广泛分布于两个吡唑氮和酚氧，正电势区域更明显，静电环境更异质、极化更强，可能支持更多瞬态接触，但不直接证明稳定机制。
theoretical (TDDFT) electronic spectral analysis：TD-DFT显示两者紫外吸收差异。L1最强跃迁236.88 nm（f=0.2679），S1为269.36 nm（4.603 eV，f=0.0982，HOMO→LUMO），S3为248.93 nm（4.98 eV，f=0.012，HOMO–2→LUMO）。L2主跃迁268.99 nm（4.61 eV，f=0.147，HOMO→LUMO），248.05 nm（5.00 eV，f=0.140，HOMO→LUMO+2），221.82 nm（5.59 eV，f=0.030）。这些结果仅描述孤立配体激发性质，不用于推断受体选择性或抑制。
Natural Bond Orbital (NBO) Analysis：NBO显示L1具有较平衡、互联的π共轭和电荷转移网络，包括N₂–C₅→N₄–C₈（25.21 kcal mol^?1）、C₆–C₇→N₄–C₈（30.39）、N₂–C₅→C₆–C₇（40.26）及O₁→C₁₂–C₁₄（24.99）、N₃→C₉–C₁₃（28.49）等LP→π*作用；L2显示更强单个LP→π*作用，如N₄→N₆–C₁₅（26.94）、N₄→C₁₁–C₁₃（43.21）、N₃→C₁₀–C₁₂（42.37），表明更强的分子内电荷离域和电子极化。NBO稳定化能描述孤立体分子内给体–受体稳定，不等于蛋白–配体结合强度。
讨论部分总结：原文讨论指出，已报道合成化合物的实验生物学评价显示抗增殖活性且不同癌细胞系响应有差异，但细胞抗增殖/活力实验不能直接证明游离配体对EGFR或HER2的受体特异性结合或抑制，可能涉及其他靶标或非特异机制。计算上，L1在EGFR结合位点动态稳定性更高，L2在HER2结合位点动态稳定化更强，电子特征提供补充信息但不能建立直接因果。需要直接生化EGFR和HER2结合或激酶抑制实验、相关激酶比较测试、独立MD重复、以及细胞/生理相关模型（如organoid、无标记细胞反应监测、三维生物打印癌症模型、正交实验）验证。
结论翻译：在本研究中，两种吡唑基配体L1和L2通过整合ADMET预测、分子对接、MD模拟、DFT、TD-DFT和NBO分析的综合计算策略，被综合评价为EGFR和HER2候选配体。结果显示，尽管两种配体均对EGFR和HER2表现出有利的预测相互作用，但其计算观察到的动态行为因受体环境不同而不同。MD模拟表明，L1与EGFR形成更稳定和持久的复合物，表现为较低结构波动、稳定相互作用网络和模拟期间持续相互作用。相反，尽管L2对接亲和力较低，但在HER2结合口袋中表现出增强的动态稳定化，表明受体依赖构象适应和相互作用持久性不能仅从对接分数完全推断。这些发现突出长时间尺度MD模拟在表征模拟条件下靶标依赖结合和稳定化行为方面的价值。电子结构分析为两种配体的内在电子特征提供额外信息。L1表现出相对局域电荷分布、有效分子内电荷转移、方向性电子响应和平衡π共轭模式；L2表现出更大轨道离域、更高极化率、更强单个LP→π*稳定化作用和更大电子灵活性。这些差异提供补充电子信息，可能有助于解释两种配体不同的计算观察到的相互作用和稳定化行为，但未建立电子描述符与受体结合的直接机制关系。总体而言，这些发现表明这些吡唑衍生物内的结构和电子差异与不同的计算观察到的受体依赖结合和动态稳定化行为相关。综合MD和电子结构分析为解释L1和L2观察到的不同EGFR和HER2相关计算行为提供可能框架，并强调在激酶靶向配体计算研究中考虑受体动力学的重要性。然而，由于未实验测定受体特异性结合对接分数和抑制活性，所观察差异应解释为计算预测的靶标依赖结合和稳定化行为，而非确认的分子选择性或EGFR/HER2抑制。此外，先前报道的细胞抗增殖活性不能视为游离配体EGFR或HER2特异性活性的直接证据。因此，需要使用游离配体进行直接生化结合和功能性激酶抑制研究，以确定预测相互作用是否对应实验可测的受体活性，以及计算观察到的差异是否转化为受体特异性分子活性。