《Results in Chemistry》：Computational insights into Baicalein as a potential inhibitor of EGFR T790M: Comparative analysis with Osimertinib
编辑推荐：
非小细胞肺癌（Non-small cell lung cancer, NSCLC）仍然是癌症相关死亡的主要原因之一，其中表皮生长因子受体（Epidermal Growth Factor Receptor, EGFR T790M）突变是对第一代和第二代酪氨酸激酶
非小细胞肺癌（Non-small cell lung cancer, NSCLC）仍然是癌症相关死亡的主要原因之一，其中表皮生长因子受体（Epidermal Growth Factor Receptor, EGFR T790M）突变是对第一代和第二代酪氨酸激酶抑制剂（Tyrosine Kinase Inhibitors, TKIs）产生耐药的主要机制。在本研究中，研究人员采用全面的计算机模拟（in silico）方法，系统研究了天然黄酮类化合物黄芩素（Baicalein）与EGFR T790M的结合行为，并与第三代抑制剂奥希替尼（Osimertinib）进行了比较。进行了分子对接（Molecular Docking）、分子动力学模拟（Molecular Dynamics Simulations, MD）、MM/PBSA结合自由能计算以及ADMET预测。对接分析显示，黄芩素与EGFR结合口袋中的关键残基（包括MET793和GLY796）形成了稳定的相互作用。分子动力学模拟证明了黄芩素-EGFR复合物的结构稳定性。在应用的MM/PBSA框架内，黄芩素表现出比奥希替尼（?30.71 ± 14.14 kJ/mol）更有利的结合自由能（?43.35 ± 13.56 kJ/mol）。计算机模拟ADMET分析表明，黄芩素具有良好的药代动力学和安全性相关性质。总体而言，这些发现突出了黄芩素作为一种有前景的天然骨架，值得在EGFR T790M抑制的背景下进一步探索。
**论文解读文章**
**1. 研究背景与问题**
非小细胞肺癌（NSCLC）是全球癌症相关死亡的主要原因之一，其中表皮生长因子受体（EGFR）信号通路异常是关键驱动因素。第一、二代酪氨酸激酶抑制剂（TKIs）虽有效，但EGFR T790M突变（第790位苏氨酸被甲硫氨酸取代）导致获得性耐药，成为临床治疗瓶颈。第三代药物奥希替尼（Osimertinib）通过不可逆结合靶向T790M突变，被列为一线用药，但其使用伴随继发性耐药及严重副作用，如心脏毒性、呼吸系统和肝胆系统疾病。因此，亟需探索具有更低毒性的替代天然化合物。黄芩素（Baicalein）是从黄芩根中提取的天然黄酮类化合物，具有抗氧化、抗炎、抗肿瘤等活性，可抑制包括EGFR在内的酪氨酸激酶受体，并调控PI3K/Akt和MAPK/ERK通路。本研究旨在通过计算机模拟（in silico）方法，系统比较黄芩素与奥希替尼与EGFR T790M的分子相互作用，评估黄芩素作为天然抑制剂的可能性。该论文发表在《Results in Chemistry》。
**2. 关键技术方法**
研究人员采用以下主要技术方法（未涉及具体试剂、培养及质粒操作）：
- **分子对接**：使用AutoDock Vina和UCSF Chimera 1.17，将EGFR T790M受体（基于PDB ID: 3W2S晶体结构，手动引入T790M突变）设为刚性，配体设为柔性，基于CASTp 3.0预测活性位点。
- **分子动力学模拟**：使用GROMACS 2025.2和CHARMM36力场，对黄芩素-EGFR T790M和奥希替尼-EGFR T790M复合物进行70 ns模拟，包含NVT和NPT平衡，分析RMSD、RMSF、回转半径（Rg）和溶剂可及表面积（SASA）。
- **结合自由能计算**：通过gmx_mmpbsa工具，采用MM-PBSA方法从模拟轨迹中提取1000帧计算结合自由能（ΔGbind）。
- **ADMET预测**：使用SwissADME、pkCSM和ProTox 3.0工具预测药代动力学、药物相似性、心脏毒性、肝毒性、致突变性等参数。
**3. 研究结果**
**3.1. 分子对接研究**
通过分子对接比较黄芩素与奥希替尼与EGFR T790M的结合能力。结果显示，奥希替尼对接分数为?7.8，黄芩素为?7.3，差异较小。结构分析表明，黄芩素与活性位点关键残基MET793A和GLY796A形成三个稳定氢键，并与LEU718A、VAL726A、LEU844A产生四个疏水相互作用。奥希替尼则与LEU718A和CYS797A形成两个氢键，与七个疏水残基相互作用，并存在ASP800A的盐桥。这说明黄芩素作为天然化合物能与EGFR T790M有效结合，使其成为有潜力的候选物。
**3.2. 分子动力学模拟**
对两种复合物进行70 ns MD模拟，评估结构稳定性。RMSD分析显示，两种复合物在最初几纳秒后均达到稳定状态，表明黄芩素和奥希替尼均稳定了蛋白质结构。RMSF结果指出，在部分柔性区域，奥希替尼复合物波动更高，表明蛋白质在奥希替尼存在下更灵活。回转半径（Rg）在黄芩素复合物中数值更低，说明蛋白质在黄芩素存在下保持更紧凑稳定的结构。总能量分析显示两种复合物均处于热力学稳定水平。SASA分析表明，奥希替尼复合物在模拟后期SASA略微下降，相对稳定性略高。非共价相互作用分析显示，黄芩素在整个模拟过程中接触更均匀，形成2-3个稳定氢键，而奥希替尼的氢键更少、更不稳定。径向分布函数（RDF）结果显示，黄芩素在短距离内峰值更高、曲线下面积更大，表明其与结合位点关键残基的相互作用更强。总体而言，黄芩素在维持蛋白质紧凑性、形成稳定接触和氢键方面优于奥希替尼。
**3.3. 通过MM-PBSA估算结合自由能**
MM-PBSA计算提取了范德华能、静电能、极性溶剂化能、非极性SASA能及总结合能。数据表明，奥希替尼的范德华能（?160.21 ± 3.05 kJ/mol）和静电能（?99.54 ± 13.38 kJ/mol）更负，疏水和静电作用更强。但黄芩素的极性溶剂化能（108.24 ± 3.60 kJ/mol）和非极性SASA能（74.18 ± 2.76 kJ/mol）均低于奥希替尼（122.42 ± 3.26和106.65 ± 1.58 kJ/mol），表明黄芩素在结合位点中更紧凑，不利溶剂效应更小。最终总结合自由能显示，黄芩素（?43.35 ± 13.56 kJ/mol）比奥希替尼（?30.71 ± 14.14 kJ/mol）更负，说明黄芩素具有更有利的热力学结合稳定性。
**3.4. 药代动力学和毒理学预测**
通过SwissADME、pkCSM和ProTox 3.0综合分析。黄芩素具有高胃肠吸收、良好水溶性，满足所有Lipinski类药性规则（0项违反），不是P-糖蛋白（P-gp）底物，不抑制CYP3A4。奥希替尼虽高吸收，但为P-gp底物且抑制CYP3A4，可能引起耐药和代谢干扰。毒性方面，黄芩素无非预期hERG I/II抑制、无肝毒性、无免疫毒性、无神经毒性，但预测有致突变性和呼吸道毒性。奥希替尼预测有hERG II抑制（心脏毒性风险）、肝毒性、免疫毒性和神经毒性。因此，黄芩素在药代动力学和安全性方面具有显著优势，但致突变性和呼吸道毒性需进一步实验验证。
**4. 讨论与结论**
讨论部分总结：分子对接和MD模拟结果与以往研究一致，黄芩素通过氢键和疏水相互作用稳定结合EGFR T790M，其多酚结构和芳香环在ATP结合位点中发挥关键作用。MM-PBSA结果进一步支持黄芩素具有更有利的结合自由能，可能源于更紧凑的排列和更合适的氢键模式。ADMET分析显示黄芩素相比奥希替尼在心脏毒性、肝毒性方面更安全，但致突变性需关注。研究强调，这些计算机模拟发现为黄芩素作为EGFR T790M天然抑制剂奠定了基础，但任何关于抑制效力的结论均需实验验证。
**研究结论翻译**：
总体而言，本研究的计算机模拟结果表明，黄芩素具有作为靶向EGFR T790M的天然抑制剂的潜力，其预测的安全性和稳定性特征良好；然而，在得出关于比较抑制效力的任何结论之前，这些计算机模拟观察结果需要实验验证。此外，如果未来的实验研究证实其有效性，黄芩素与第三代EGFR抑制剂（如奥希替尼）的联合使用可能提供一种新的治疗策略，以克服耐药性并降低全身毒性。特别是，在细胞和体内模型中研究这两种化合物的潜在协同效应，可能为设计针对T790M突变NSCLC患者更安全、更有效的联合疗法开辟新途径。
