超越组蛋白去甲基化：非血红素Fe(II)/2-氧戊二酸依赖型加氧酶KDM6B催化N-烷基连续氧化的机制

《JACS Au》：Beyond Histone Demethylation: Mechanisms of N-Alkyl Consecutive Oxidations by the Non-Heme Fe(II)/2-Oxoglutarate Oxygenase KDM6B

　　

非血红素Fe(II)/2-氧戊二酸(2-oxoglutarate, 2OG)依赖型组蛋白去甲基化酶KDM6B(JMJD3)对组蛋白H3肽段中Nε-烷基化赖氨酸残基表现出多样的氧化转化能力，然而相较于典型的KDM催化的去甲基化反应，此类脱烷基作用的机制尚不清楚。研究人员采用分子动力学(molecular dynamics, MD)与量子力学/分子力学(quantum mechanics/molecular mechanics, QM/MM)方法，研究了KDM6B对H3肽链中不同N-烷基化形式的赖氨酸K27——即Nε,Nε-甲基乙基赖氨酸(Lys(Me/Eth))和Nε-异丙基赖氨酸(Lys(iPr))——进行连续氧化反应所采用的催化策略。针对生成醇、醛及随后羧酸产物的连续氧化结果表明，不同N-烷基化基团构象定位的差异由第二配位层(second coordination sphere, SCS)相互作用和长程关联运动所介导。具体而言，不同N-烷基化基团向活性Fe(IV)=O中间体接近以实现羟基化的过程，受SCS相互作用网络调控，尤其是涉及N344和Y239的相互作用，这一点通过对N344A和Y239A突变体的MD与QM/MM计算也得到了证实。醇至醛衍生物及酸衍生物的后续氧化同样受羟基化/醛取代基构象定位的引导。QM/MM计算预测区域与化学选择性氧化可通过涉及σ或π机制的氢原子转移(hydrogen atom transfer, HAT)引发。上述见解可为设计具有非天然催化活性和改变底物选择性的Fe(II)/2OG酶提供指导，同时所揭示的机制特征可被用于设计肽类药物选择性官能团化的生物催化平台。

《超越组蛋白去甲基化：非血红素Fe(II)/2-氧戊二酸依赖型加氧酶KDM6B催化N-烷基连续氧化的机制》论文解读

研究背景与立项依据

C–H键活化是生物与合成催化中的核心挑战，非血红素Fe(II)/2-氧戊二酸(2-oxoglutarate, 2OG)依赖型加氧酶可通过高价Fe(IV)=O（ferryl）中间体实现高选择性C–H氧化。Jumonji C(JmjC)结构域组蛋白去甲基化酶（lysine demethylases, KDMs）通常催化组蛋白H3上N^ε-甲基赖氨酸的去甲基化，但已有实验表明人源KDM6B（JMJD3）还能催化N^ε-烷基化赖氨酸（如N^ε,N^ε-甲基乙基赖氨酸Lys(Me/Eth)和N^ε-异丙基赖氨酸Lys(iPr)）的脱烷基及连续氧化（醇→醛→羧酸）。与TET2或AlkB家族Fe(II)/2OG酶的逐步氧化已有较深入计算研究不同，KDM6B催化非天然N-烷基底物连续氧化的构象动力学与电子转移机制尚不明确。本文拟通过计算揭示KDM6B对N-烷基赖氨酸底物的区域/化学选择性连续氧化机制及其第二配位层(second coordination sphere, SCS)残基调控作用，为酶工程及肽类药物生物催化修饰提供理论基础。本研究发表于《JACS Au》。

主要关键技术方法

研究人员以PDB ID 5oy3的人源KDM6B晶体结构为基础，将活性中心建模为Fe(IV)=O与琥珀酸配位的ferryl态，并将H3K27二/三甲基赖氨酸修改为Lys(Me/Eth)、Lys(iPr)及相应的羟基化/醛基化中间体。力场参数通过Metal Center Parameter Builder(MCPB.py)及GAFF分配，Zn(II)结合域用Zinc Amber Force Field(ZAFF)处理。进行全原子分子动力学(MD)模拟（Lys(Me/Eth)及羟基化变体模拟2 μs，其余1 μs，NPT系综，300 K，1 bar），经K-means聚类取代表性快照进行QM/MM计算。QM区域含facial triad残基(H250/E252/H330)、Fe(IV)=O、琥珀酸及截短的烷基化赖氨酸侧链，以高自旋五重态(quintet)用B3LYP/def2-SVP优化，单点能校正用def2-TZVP，零点能校正来自频率计算，过渡态用dimer法确认。辅以主成分分析(principal component analysis, PCA)、动态互相关分析(dynamic cross-correlation analysis, DCCA)、核密度估计(kernel density estimate, KDE)、自旋自然轨道(spin natural orbital, SNO)分析、碎片分子轨道(fragment molecular orbital, FMO)分析及能量分解分析(energy decomposition analysis, EDA)；并对Y239A与N344A突变体做对照MD与QM/MM计算。

研究结果

3.1. 构象动力学如何促进KDM6B–Fe(IV)=O与N-烷基化H3-赖氨酸底物形成反应复合物

MD显示Lys(Me/Eth)甲基比乙基更常靠近Fe(IV)=O（分别占65%和35%），与实验优先脱甲基相符；乙基Cα（羟化后导致脱乙基产乙醛）与Cβ（连续氧化起始）均可在ferryl氧可及范围，略偏Cα。Lys(iPr)的异丙基Cα（叔碳）比甲基Cβ/Cβ′更常靠近Fe(IV)=O，与实验优先脱异丙基一致。SCS残基N344与底物N^ε-H形成氢键稳定甲基或异丙基Cα朝向构象；Y239与底物N^ε-H形成氢键稳定乙基朝向构象。PCA与DCCA显示不同底物引起KDM6B柔性区（β15–β17连接区及152–184无序环）及与双链β螺旋(double-stranded beta helix, DSBH)折叠、Zn(II)结合域间的关联/反关联运动差异，影响底物定位与产物分布。

3.2. KDM6B催化羟基化机制

从MD快照做QM/MM：Lys(Me/Eth)甲基HAT能垒23.0 kcal/mol低于乙基Cβ-HAT（28.1 kcal/mol），而乙基Cα-HAT（当乙基靠近ferryl时）仅17.5 kcal/mol，说明甲基优先去甲基化，乙基Cα羟化导致脱乙基。Lys(iPr)异丙基Cα-HAT能垒20.2 kcal/mol低于甲基Cβ-HAT（26.6 kcal/mol），与优先脱异丙基的实验一致。SNO与FMO分析表明HAT电子转移可经σ(z2)（σ-途径，∠Fe–O–H≈120°–180°）或π(xz/yz)（π-途径，∠Fe–O–H≈90°–120°），取决于底物取向与Fe–O–C角度分布（KDE分析证实）。E252非键合氧与底物N^ε-H的氢键（甲基朝向时）对甲基HAT至关重要，破坏后乙基HAT发生。

3.2.1.1. Y239A和N344A重塑底物取向与反应倾向

N344A突变使甲基靠近ferryl的比例升至95%，甲基HAT能垒降至21.1 kcal/mol，强化去甲基化偏好。Y239A突变使乙基靠近ferryl的比例升至88%（多为Cα朝向），乙基Cα-HAT能垒22.7 kcal/mol低于甲基HAT（30.2 kcal/mol），转向偏好脱乙基。Q237在两种突变体中部分替代Y239/N344发挥"门控"氢键作用。

3.2.3. Lys(Me/Eth)和Lys(iPr)羟基化的电子与构象特征

Lys(Me/Eth)甲基HAT走π-途径（π(yz)），乙基Cα/Cβ-HAT走σ-途径（σ(z2)）；Lys(iPr) Cα-HAT走σ-途径，Cβ-HAT走π-途径。EDA指出N191、D291、K241、S258、N253、Q237等不同组合稳定特定HAT过渡态，这些残基的运动与DSBH及Zn结合域呈特征性相关/反相关。

3.3. 羟基化烷基底物在KDM6B–Fe(IV)=O态的构象动力学

羟基化变体Lys(Me/Eth–OH)和Lys(iPr–OH)的羟基化侧链Cβ更易靠近Fe(IV)=O；N344与Lys(iPr–OH)底物N^ε-H保持氢键，但与Lys(Me/Eth–OH)作用减弱，使甲基难以接近ferryl。无序区(152–184)柔性增加并与DSBH呈强反相关，利于后续氧化步骤的底物结合。

3.4. KDM6B中羟基化烷基底物的后续氧化机制（醇→醛）

对Lys(Me/Eth–OH)：先O_oh–H（羟基质子的HAT）能垒20.9 kcal/mol生成氧中心自由基(Fe(III)–OH)，再Cβ–HAT能垒7.3 kcal/mol得醛；也可先Cβ–HAT（24.4 kcal/mol）再O_oh–HAT（1.3 kcal/mol），两途径均可行但直接O-HAT先行稍有利。水分子介导的HAT在前步相当、在后步不利；偕二醇(gem-diol)中间体形成后脱水能垒过高，故醛生成主要通过第二次HAT而非偕二醇路径。对Lys(iPr–OH)：Cβ–HAT先行（17.9 kcal/mol）再O_oh–HAT（~0.4 kcal/mol）比先O-HAT（24.1 kcal/mol）有利，与Lys(Me/Eth–OH)顺序相反。S258与底物羰基/羟基形成氢键参与稳定。SNO显示醇氧化HAT可走σ-或π-途径取决于Fe–O–X角度。

3.5. 醛基类底物在KDM6B–Fe(IV)=O态的构象动力学

Lys(Me/Eth–AL)和Lys(iPr–AL)的醛氢最靠近Fe(IV)=O，符合下一步氧化要求；N344与Lys(iPr–AL)底物N^ε-H保持氢键，Lys(Me/Eth–AL)则视构象分别与Y239或N344作用。MM/GBSA结合能暗示醇/醛中间体比原始底物结合更强，支持进程性(processive)连续氧化。

3.6. KDM6B催化醛基取代底物的氧化（醛→羧酸）

Lys(Me/Eth–AL)醛HAT能垒17.0 kcal/mol（σ-途径），OH回弹(rebound)能垒12.6 kcal/mol生成羧酸，终态与E252非键合氧形成氢键稳定。Lys(iPr–AL)醛HAT能垒20.9 kcal/mol（π-途径），OH回弹近无垒。S258与醛羰基氧形成氢键协助定位。EDA示E104、S258、N260、N344、Q237等按底物不同参与过渡态稳定。

3.7. hTET2与KDM6B连续氧化机制比较

两者均经Fe(IV)=O介导的逐步HAT实现烷基→醇→醛→酸转化。hTET2中5hmC→5fC的HAT先后顺序存在计算模型依赖性争议；KDM6B则明确表现底物依赖性顺序：Lys(Me/Eth)醇→醛优选先O-HAT后C-HAT，Lys(iPr)醇→醛优选先C-HAT后O-HAT，源于SCS相互作用与静电环境差异导致底物取向与可接受轨道不同。

讨论与结论翻译（摘自原文浓缩）

本研究联合经典MD与QM/MM计算阐明了非血红素Fe(II)/2OG依赖型加氧酶KDM6B对合成修饰赖氨酸底物区域选择性及连续氧化的机制细节。模拟显示Lys(Me/Eth)甲基比乙基更易接近Fe(IV)=O，Lys(iPr)异丙基Cα比甲基Cβ更易接近。N344、E252和Y239决定N-烷基在Lys(Me/Eth)和Lys(iPr)中的取向。KDE分析表明较高∠Fe–O–H利于σ-途径HAT，较低∠Fe–O–H利于π-途径HAT，几何偏好反映于QM/MM优化的过渡态并由SNO与自旋密度分析证实。FMO分析σ(z2)与π(xz)间能隙有助于理解电子转移途径与活化能垒关联。EDA鉴定出N344、K241、S258、N260和D291等在过渡态稳定/去稳定中具底物特异性作用；DCCA显示这些残基与DSBH折叠相关且与Zn(II)结合域反相关，强调其别构影响。对Y239A与N344A突变体的MD与QM/MM表明这些SCS残基突变可改变活性部位底物定位及产物比例。羟基化与醛基化中间体的模拟显示酶可调节活性部位与构象动力学以支持后续氧化；S258与醛羰基形成稳定氢键促进氧化至相应酸。结果表明KDM6B催化受底物取向、电子转移机制及蛋白–底物动态相互作用共同调控，区域选择性与连续氧化进程由局部化学环境与涉及DSBH折叠及Zn(II)结合域的长程关联运动调制。鉴于KDM6B的反应多样性，其可被改造为设计修饰模拟H3尾序列合成肽的生物催化工具，该研究为拓展非天然非血红素Fe(II)/2OG酶催化谱系及指导酶工程利用KDM6B作为肽类药物选择性官能团化生物催化平台提供了机制基础。