使用密度泛函理论(DFT)对O6-烷基鸟嘌呤-DNA-烷基转移酶(AGT)抑制剂的硝基类前体药物在硝基还原酶(NTR)作用下的还原机理进行了研究。结果显示,在反应的六个过渡态中,生成亚硝基反应的过渡态(TS2)和生成苯胺反应的过渡态(TS5)为反应的两个控速步骤。整个还原反应过程中水参与的质子电子耦合转移在过渡态TS2和TS6降低了反应的活化能,推测水介导的质子电子耦合转移为这两步反应的优势途径。通过绘制反应势能曲线阐明了硝基芳烃类前体药物还原机理及其反应所涉及的热力学问题,为制备高效低毒的AGT抑制剂的前体药物提供理论基础。 The reduction mechanism of O6-Alkylguanine DNA Alkyltransferase inhibitor of the nitro group of the pro-drug in the Nitro Reductase (NTR) was studied by using the Density Functional Theory (DFT). The results showed that the transitions of TS2 (yielding nitroso compound) and TS5 (yielding aniline compound) were the rate-limiting steps among the six reaction transition states. Single water molecule-mediated proton coupled electron transfer can lower the reactive activation energy of TS2 and TS6; thus it is speculated that water-mediated proton coupled electron transfer is the optimal reaction pathway of TS2 and TS6. The reduction mechanism of the prodrug of nitro aromatic compound and the involved reaction thermodynamics can be clearly presented by the potential energy curve. This study can provide a theoretical basis for designing more efficient hypoxia-activated AGT inhibitors with less side effects.
肖家斌1,2*,孙国辉1,2
1北京工业大学生命科学与生物工程学院,北京
2环境与病毒肿瘤学北京市重点实验室,北京
收稿日期:2016年4月15日;录用日期:2016年5月2日;发布日期:2016年5月5日
使用密度泛函理论(DFT)对O6-烷基鸟嘌呤-DNA-烷基转移酶(AGT)抑制剂的硝基类前体药物在硝基还原酶(NTR)作用下的还原机理进行了研究。结果显示,在反应的六个过渡态中,生成亚硝基反应的过渡态(TS2)和生成苯胺反应的过渡态(TS5)为反应的两个控速步骤。整个还原反应过程中水参与的质子电子耦合转移在过渡态TS2和TS6降低了反应的活化能,推测水介导的质子电子耦合转移为这两步反应的优势途径。通过绘制反应势能曲线阐明了硝基芳烃类前体药物还原机理及其反应所涉及的热力学问题,为制备高效低毒的AGT抑制剂的前体药物提供理论基础。
关键词 :AGT抑制剂,前体药物,硝基还原酶,质子电子耦合转移
氯乙基亚硝基脲(CENUs)是一类重要的抗癌烷化剂,是目前临床上治疗恶性肿瘤的化疗药物之一。CENUs通过导致DNA生成股间交联(dG-dC交联),阻断肿瘤细胞DNA的正常复制,诱导肿瘤细胞凋亡而发挥抗癌作用 [
肿瘤细胞的快速增殖导致高氧消耗,肿瘤脉管系统的结构功能异常导致供氧不足,使得瘤内低氧是实体瘤的重要特征之一。因此,设计合成具有生物还原性的前体药物,使其依赖于肿瘤低氧的微环境,在低氧条件下激活为细胞毒素并靶向作用于肿瘤细胞,已成为提高化疗药物疗效的新思路。研究表明,生物还原性前体药物的激活主要通过单/双电子还原反应来实现。巴诺蒽醌在低氧细胞内通过细胞色素P450催化的两步电子还原反应最终转变为具有细胞毒性的DNA结合剂,可特异性地抑制DNA拓扑异构酶II [
首先,以硝基苯为模型,研究了其在低氧条件下由FMNH催化的还原反应(如图1所示)进行了研究。在此基础上,对图2所示的O6-4-硝基苄基鸟嘌呤(O6-4-NBG)在FMNH作用下分别被还原为具有AGT抑制活性的O6-4-氨基苄基鸟嘌呤(O6-4-ABG)的作用机理进行了研究。计算采用DFT方法在UB3LYP/6-31+ G(d,p)水平上进行,对所有反应物、过渡态、中间体和产物进行全几何结构优化。在相同水平上对优化得到的各驻点几何结构进行振动频率计算,确认所有稳定点无虚频;所有过渡态只含有一个虚频,并且频率的振动模式为沿着反应坐标方向的伸缩振动。计算内禀反应坐标以确认所有过渡态结构的合理性,并获得反应的最低能量途径。所有计算均采用Gaussian09程序包进行。
如图1所示,硝基苯在FMNH的作用下被还原为苯胺的反应总共需要六个电子和六个质子的过程,它可以分解为三个步骤,这意味着存在两个中间体。在第一步还原反应中,硝基基团通过2e−/2H+的转移还原成中间体I(亚硝基),同时消去一分子水。在第二步还原反应中,亚硝基通过2e−/2H+的转移过程转化为中间体II(羟胺)。最后,羟胺通过2e−/2H+的转移过程转化为苯胺,同时消去一分子水 [
图1. FMNH催化硝基苯还原的反应机理
图2. FMNH还原O6-4-硝基苄基鸟嘌呤的反应路径
FMNH-CH3),并将使用未简化和简化后的模型计算得到能量进行了比较,以确认简化模型的合理性。以FMNH-CH3为模型采用GaussView 5.0计算得到的各步骤的过渡态结构如图3所示。在六个反应过渡态中FMNH的N5位上的N-H键与初始键长(1.01 Å)相对增长,表明FMNH的N5位上的氢有断键趋势。其余各键中氧氢键的距离减小,氮氧键增长及新的氮氢键的形成,且各反应过渡态都存在唯一虚频。
在B3LYP/6-31+G(d,p)水平上以RC(0.00)为势能曲线基准采用吉布斯自由能绘制的势能曲线如图4所示。在整个反应过程中,只有第一步羟基化反应生成IC1的步骤为吸热反应,其它五个步骤均为放热反应。其中TS2反应所需跨越的能垒最高(30.11 kcal/mol),因此经第二步脱水反应生成亚硝基苯的反应为控速步骤,并且反应放热达31.56 kcal/mol。TS3反应所需过渡态能量最低,亚硝基比硝基基团更容易发生反应,因此亚硝基苯的含量很少。TS5反应所需过渡态能量次于TS2,TS5可能是羟胺还原这一步反应的控速步骤,由于TS5反应较慢,因此羟胺被大量聚集,而TS6反应所需能量较低且前一步反应提供大量的热量,反应生成苯胺的速率较快。因此整个反应主要以羟胺和苯胺形式存在,这与实验结果相一致 [
此外,考虑了溶剂效应对反应的影响,研究了水分子参与反应的过程,其势能曲线如图4所示。水参与的质子电子耦合转移 [
图3. 在B3LYP/6-31+G(d,p)水平上优化得到FMNH催化硝基苯还原反应的过渡态结构
图4. B3LYP/6-31+G(d,p)水平上计算得到的FMNH催化硝基苯还原反应的势能曲线
图5. 水分子参与下硝基苯的过渡态TS2和TS6的几何结构及HOMO示意图
径,形成质子电子耦合转移。Wilkie等人的研究表明 [
对图2所示反应途径中的各驻点在B3LYP/6-31+G(d,p)水平上进行全几何结构优化和过渡态的搜索,同时得出TS2和TS6在水参与下的质子耦合电子转移的最优反应路径,其GaussView 5.0计算得到各反应的过渡态结构及结构参数如图6;计算结果表明,O6-4-硝基苄基鸟嘌呤(O6-4-NBG)在FMNH的条件下被还原成O6-4-氨基苄基鸟嘌呤(O6-4-ABG)的反应分六个步骤。在各反应中FMNH的N5位上的N-H键与初始键长(1.01 Å)相对增长,其余各键中氧氢键的距离减小,氮氧键增长及新的氮氢键的形成,且各反应过渡态都存在唯一虚频。
在B3LYP/6-31+G(d,p)水平上以RC(0.00)为势能曲线基准的采用吉布斯自由能绘制的势能曲线如图7。其在采用B3LYP方法计算所得的结构下各过渡态反应能垒依次为:18.97、30.14、9.45、15.68、24.81、16.91 kcal/mol。O6-4-NBG的反应趋势和反应的活化能与硝基苯势能曲线一致,可以看出硝基类化合物的
图6. 在B3LYP/6-31+G(d,p)水平上优化得到FMNH催化O6-4-NBG还原反应的过渡态结构
图7. B3LYP/6-31+G(d,p)水平上计算得到的FMNH催化O6-4-NBG还原反应的势能曲线
反应机理与硝基所连的基团并没有太大影响。整个反应存在两个控速步骤分别生成亚硝基的过渡态TS2和生成氨基的过渡态TS5,整个反应主要以羟胺化合物、苯胺化合物的形式存在,亚硝基化合物被瞬间反应含量很少。其中水分子参与过渡态反应优化了反应能垒的TS2和TS6反应所需能量为29.75 kcal/mol和14.13 kcal/mol,分别降低了0.39 kcal/mol和2.78 kcal/mol,因此这两步水分子参与反应降低了反应的活化能为该步反应的优势途径。采用GaussView 5.0计算得到水参与质子转移的过渡态结构和结构参数及过渡态下的HOMO轨道示意图如图8所示。其结果表明,在TS2-H2O中,质子从FMNH的N5位转移
图8. 水分子参与下O6-4-NBG的过渡态TS2和TS6的几何结构及HOMO示意图
到水的一个氧原子上(N…H…O),然后通过水转移质子到羟基化的氧上(O…H…O),并脱去两个水分子为质子转移通道。在TS2-HOMO中,FMNH的N5的离域与硝基苯的N的离域两原子轨道重叠(N5…N),因此在FMNH的N5位和硝基苯的N之间作为电子传递途径。在TS6-H2O中,质子从FMNH的N5位转移到水的一个氧原子上(N…H…O),然后通过水转移质子到IC5的N原子上(O…H…N),并脱去一个水分子为质子转移通道。在TS6-HOMO中,FMNH的N5的离域与中间体IC5的N的离域两原子轨道重叠(N5…N),因此在FMNH的N5位和IC5的N之间作为电子传递途径,形成质子电子耦合转移。过渡态TS2和TS6降低了过渡态反应的活化能分别为这两步反应的优势路径。
综上所述,在B3LYP/6-31+G(d,p)水平上进行全几何结构优化和过渡态的搜索的计算结果表明,硝基类化合物能够在低氧条件下被FMNH催化还原生成氨基类化合物。在整个还原反应中过渡态TS2和过渡态TS5分别为生成亚硝基和生成苯胺的两个反应的控速步骤。由于水分子在生物体内大量存在,考虑水分子对整个反应的影响,得出由水介导的过渡态TS2和TS6为这两步催化还原反应的优势途径。综上所述,通过密度泛函的方法,阐明了在低氧条件下NBGs经过硝基还原酶和辅酶FMNH介导的单电子传递过程和还原反应机理,为低氧激活靶向性抗肿瘤药物的开发提供了可靠的理论依据。
肖家斌,孙国辉. FMNH催化O6-4-硝基苄基鸟嘌呤还原机理的研究 Mechanism Study of the Reduction of O6-4-Nitrobenzyl-Guanine Catalysed by FMNH[J]. 药物资讯, 2016, 05(02): 25-32. http://dx.doi.org/10.12677/PI.2016.52005