编辑推荐：
本文通过计算机模拟与实验验证相结合的方法，揭示了天然抗氧化剂去甲二氢愈创木酸（NDGA）与人源花生四烯酸脂氧合酶15（ALOX15）的相互作用机制。研究聚焦于催化口袋关键氨基酸Gln595，通过同源建模、分子对接、分子动力学（MD）模拟及Gln595点突变体（Ala、Leu、Ile、Glu）的体外功能实验，证实Gln595对NDGA-ALOX15复合物的稳定性和抑制效能具有决定性影响。该研究为未来开发亚型特异性ALOX15抑制剂提供了重要的结构基础与理论依据。
1. 引言
哺乳动物花生四烯酸脂氧合酶（ALOXs）是一类含非血红素铁的脂质过氧化酶，在炎症、过度增殖及神经退行性疾病中扮演重要角色。去甲二氢愈创木酸（NDGA）是一种天然抗氧化剂，同时也是强效的脂氧合酶抑制剂。然而，NDGA与ALOX相互作用的分子基础尚未明确。在人类多种ALOX亚型中，ALOX15在铁死亡、炎症及免疫应答中处于中心地位，其异常表达与动脉粥样硬化、神经退行性变、中风、糖尿病等多种病理状况相关。尽管已发现包括NDGA在内的多种天然与合成ALOX15抑制剂，但其确切的作用机制，特别是与酶活性位点的具体相互作用细节，仍有待阐明。本研究旨在通过计算机模拟与实验相结合的策略，深入探究NDGA与人源ALOX15（h-ALOX15）结合的分子基础，尤其关注关键氨基酸残基的作用。
2. 材料与方法
研究采用了多层次的研究策略。首先，由于人源ALOX15的晶体结构尚未解析，研究者利用与其具有高度同源性（>85%）的兔源ALOX15晶体结构（PDB ID: 2P0M）作为模板，通过Schr?dinger套件的Prime模块进行了h-ALOX15的同源建模。NDGA的3D结构从ChEMBL数据库获取并优化。随后，使用GOLD软件将NDGA对接到h-ALOX15的底物结合口袋中，确定其可能的结合构象。为评估复合物稳定性，对排名最高的对接构象进行了长达500纳秒的分子动力学（MD）模拟。
基于对接和MD模拟结果，研究锁定了活性位点中的一个关键氨基酸——谷氨酰胺595（Gln595）。为验证其功能，研究通过in silico（计算机模拟）突变研究了四种点突变体：Gln595Ala、Gln595Leu、Gln595Ile和Gln595Glu，并同样进行了MD模拟分析。与此同时，在实验层面，研究者将野生型h-ALOX15及上述四种突变体在大肠杆菌中重组表达为N端六组氨酸标签融合蛋白。通过定量蛋白质免疫印迹评估表达水平，并使用花生四烯酸（AA）体外氧合酶活性实验来表征各突变体的功能特性。最后，通过测定NDGA对野生型及各突变体酶活的半抑制浓度（IC50），验证了计算预测的结果。
3. 结果
3.1. NDGA结合于人源ALOX15活性位点且Gln595具有机制相关性
分子对接显示，NDGA结合在h-ALOX15的底物结合口袋内，周围涉及包括Glu175、Phe352、Ile592、Gln595、Ile662在内的多个氨基酸残基。NDGA的两个儿茶酚环（cat-A和cat-B）通过烷基链连接。在活性位点中，cat-A环靠近催化铁，并与Ile662的骨架原子形成氢键；cat-B环则与Ile592的骨架原子及Gln595的侧链形成氢键。这些相互作用表明Gln595在NDGA结合中可能起重要作用。
3.2. ALOX15–NDGA复合物稳定性与Gln595的作用
500 ns的MD模拟表明，野生型ALOX15-NDGA复合物整体稳定。特别值得注意的是，cat-B环与Gln595侧链之间的距离在模拟过程中从3.25 ?缩短至3.05 ?，提示两者间相互作用增强。而cat-A环与Ile662的相互作用在整个模拟期间保持稳定。这进一步支持了Gln595是NDGA结合的关键残基。
3.3. Gln595的in silico突变影响部分ALOX15–NDGA复合物稳定性
对Gln595进行计算机模拟突变后的MD模拟显示，所有突变体蛋白骨架在约100 ns后均达到稳定。然而，配体（NDGA）的稳定性因突变类型而异。Gln595Ala突变由于侧链体积显著减小，在活性位点产生了一个空隙，导致NDGA在模拟前期构象漂移，甚至部分移出底物结合口袋。Gln595Leu和Gln595Ile突变对复合物稳定性的影响相对温和，NDGA最终能找到替代的稳定氢键（如与Glu356或Gln589）。Gln595Glu突变引入负电荷，但复合物稳定性在模拟后期得以恢复，NDGA通过与Ile662、Gln589和Glu356形成氢键网络而稳定。
3.4. 野生型与突变体ALOX15的成功表达与功能表征
所有h-ALOX15变体均在大肠杆菌中成功表达，表达水平相当。功能实验表明，与野生型相比，大多数Gln595突变体的催化活性略有下降，但Gln595Ile突变体的活性反而高出四倍以上。产物分析显示，野生型酶主要产生15S-羟基二十碳四烯酸（HETE）和少量12S-HETE。Gln595Ala、Leu、Ile突变减少了12-HETE的生成比例，而Gln595Glu突变则增加了12-HETE的比例。手性分析证实所有酶变体均高度立体选择性地产生S型对映体，表明酶制备质量良好。
3.5. Gln595Ala和Gln595Leu突变严重削弱NDGA的抑制效能
关键的抑制实验数据显示，NDGA对野生型h-ALOX15的IC50为126 nM。Gln595Ile和Gln595Glu突变体的IC50值与野生型处于同一数量级。然而，Gln595Ala和Gln595Leu突变体的IC50值分别比野生型高出40倍和20倍以上，表明这两种突变严重损害了NDGA的抑制效力。这一结果与MD模拟的预测相符：当Gln595被Ala或Leu替代后，NDGA的结合模式不稳定或发生改变；而当被Ile或Glu替代时，NDGA能通过与邻近的其他残基（如Ile662、Gln589）形成稳定的氢键网络，从而维持较高的抑制活性。
4. 讨论
4.1. 研究的新颖性与科学进展
本研究首次在分子水平上详细阐明了NDGA与人源ALOX15的相互作用模式，并明确揭示了Gln595在这一过程中的三重功能相关性：不仅参与底物脂肪酸的定位和酶别构调控，还对抑制剂结合至关重要。序列比对显示，Gln595在哺乳动物ALOX15同源物中保守，但在其他ALOX亚型中不完全保守。其附近的Gln589和Gln547具有更高的保守性，可能在Gln595突变时起到功能补偿作用。另一个有趣的发现是Gln595Ile突变显著提升了酶的催化活性并使其产物特异性更偏向15S-HpETE，这可能是由于Ile的引入优化了底物在活性位点的朝向，使AA的C13更靠近催化铁。
4.2. 研究的局限性
本研究使用的酶制剂为粗提物，尽管已对蛋白含量进行归一化，但杂质蛋白的潜在影响无法完全排除。不过，野生型与Gln595Ala/Leu突变体间高达两个数量级的IC50差异，很难归因于杂质。此外，所有计算模拟基于兔源ALOX15晶体结构构建的人源模型，尽管两者同源性和功能高度相似，但与人源真实结构仍可能存在细微差异。
5. 结论
这项结合计算与实验的比较突变分析，凸显了Gln595在调节NDGA与人源ALOX15结合中的关键作用。研究证实，Gln595Ala和Gln595Leu突变会严重削弱NDGA的抑制效力，而Gln595Ile和Gln595Glu突变的影响较小，这可能得益于Gln589等邻近残基的补偿作用。这些发现为深入理解ALOX15与抑制剂的相互作用机制提供了新见解，并对未来设计和优化亚型特异性ALOX15抑制剂具有重要的指导意义。
