FAD/NAD(P)H依赖型脱氢酶中配体选择性性的结构基础：来自由锥虫硫酮还原酶(Trypanothione Reductase, TR)与II型NADH脱氢酶(Type II NADH Dehydrogenase, NDH2)的启示

《Protein Science》：Structural basis of ligand selectivity in FAD/NAD(P)H-dependent dehydrogenases: insights from trypanothione reductase and type II NADH dehydrogenase

　　

FAD/NAD(P)H依赖型脱氢酶(Flavoprotein dehydrogenase)构成结构保守的氧化还原酶家族，参与寄生虫及高等生物关键代谢过程。其中锥虫硫酮还原酶(Trypanothione Reductase, TR)是杜氏利什曼原虫(Leishmania infantum, Li)氧化还原代谢的关键组分，是抗利什曼药物开发的重要靶标。但因若干黄素蛋白(Flavoprotein)具相似折叠及辅因子结合架构，TR抑制剂的选择性在早期药物发现中仍是关键问题。研究人员通过整合序列分析、结构模建、对接模拟及体外生化验证，探讨了婴儿利什曼原虫中FAD/NAD(P)H依赖脱氢酶的结构图谱。互逆序列搜索鉴定出11种与TR结构相关的寄生黄素蛋白，含二氢硫辛酰胺脱氢酶(Dihydrolipoamide Dehydrogenase, DLD)样蛋白、II型NADH脱氢酶(Type II NADH Dehydrogenase, NDH2)样蛋白及二烯酰辅酶A还原酶(Dienoyl-CoA Reductase, deCoAR)样蛋白。跨寄生虫及哺乳动物同源蛋白的比较对接分析可考察这些酶潜在交叉反应模式。重组LiTR与栖热碱芽孢杆菌(Caldalkalibacillus thermarum)NDH2(CtNDH2)酶活测定证实选定配体可产生酶依赖性效应：金诺芬(Auranofin)及硝基呋喃腙(Nitrofurazone)抑制LiTR；三联吡啶-铂(Terpyridine–Pt)衍生物强烈抑制LiTR同时激活CtNDH2活性。结果表明结构相关黄素蛋白可容纳共同配体并产生不同催化结局。该整合计算与生化流程为评估FAD/NAD(P)H依赖脱氢酶家族配体选择性提供实用框架，可支持靶向寄生虫氧化还原代谢选择性调节剂的开发。

论文解读：《Structural basis of ligand selectivity in FAD/NAD(P)H-dependent dehydrogenases: insights from trypanothione reductase and type II NADH dehydrogenase》发表于Protein Science

一、研究背景与立项依据

利什曼病(Leishmaniasis)由杜氏利什曼原虫(Leishmania infantum)等引起，现有疗法(五价锑剂、两性霉素B、米替福星等)存在毒性大、费用高、需长期给药及耐药性问题。锥虫硫酮还原酶(Trypanothione Reductase, TR，EC 1.8.1.12)催化氧化型锥虫硫酮(Trypanothione disulfide, T(S)₂)还原为还原型(Trypanothione, T(SH)₂)，是利什曼原虫特有抗氧化系统核心酶，而宿主细胞使用谷胱甘肽还原酶(Glutathione Reductase, GSR)，故TR被认为是高价值物种特异性药靶。然而TR与哺乳动物及其他寄生虫的FAD/NAD(P)H依赖脱氢酶——如人源谷胱甘肽还原酶(GSR)、硫氧还蛋白还原酶(Thioredoxin Reductase, TrxR)、凋亡诱导因子(Apoptosis-Inducing Factor, AIF)、二氢硫辛酰胺脱氢酶(Dihydrolipoamide Dehydrogenase, DLD, 人源及寄生虫源)——共享FAD及NAD(P)H结合折叠，TR抑制剂极易发生脱靶(off-target)交叉反应，因此配体选择性评价是早期抗利什曼药物设计瓶颈。本研究旨在建立计算—生化联用流程，系统评估LiTR抑制剂在寄生虫内源黄素蛋白酶组及人源同源酶中的选择性谱，为合理设计高选择性TR调节剂提供依据。

二、主要关键技术方法

研究人员以NCBI nr数据库对婴儿利什曼原虫(Leishmania infantum, taxid:5671)做BLASTP搜索，查询序列来自已解析的人/大肠杆菌/酿酒酵母等FAD/NAD(P)H依赖脱氢酶晶体结构；经互逆BLASTP确认获得11条全长黄素蛋白序列。采用pGenThreader与I-TASSER进行折叠识别(Fold Recognition)，选取人源DLD(PDB 1zmd)、酵母NDH2(PDB 4g6h)、大肠杆菌deCoAR(PDB 1ps9)为模板，用Swiss PDB Viewer(SPDBV)构建Li DLD样(XP_001468025.1)、Li NDH2样(XP_001469921.1)、Li deCoAR样(XP_001468165.1)同源模建(Homology Modeling)并能量最小化。从PDB收集LiTR已知抑制剂(如Auranofin/TS8、RDS 777、JV0、MWT、H6H、WPF、BVN)及辅因子(FAD、NADPH、Trypanothione GCG、Ubiquinone衍生物UQ5/DCQ)，在LiTR晶体结构(PDB 4apn)上行重对接(Redocking)验证并定义各结合口袋格点盒(Gridbox)。将同源模建与人和寄生虫晶体结构(人AIF 4bur、人GSR 1grb、人TrxR 2j3n、人DLD 1zmd、CtNDH2 5wed)叠加后执行分子对接(Molecular Docking, AutoDock 4.2.6)。重组LiTR与CtNDH2经pET-21b(+)原核表达、Ni-NTA亲和纯化后，用紫外分光光度法监测NADPH/NADH氧化速率测定酶动力学参数(K_m、V_max)及测试化合物(100 μM，硝基呋喃腙50 μM)对酶相对活性的影响。

三、研究结果

3.1 Flavoprotein homologous sequences in L. infantum

经BLASTP及互逆BLASTP，从L. infantum基因组鉴定出11条全长FAD/NAD(P)H依赖脱氢酶序列，分属TR(XP_001462998.1)、DLD(如XP_001468025.1、XP_001466710.1、AKGDH)、NDH2(XP_001469921.1)、deCoAR/nCoAR(XP_001468165.1、XP_001463176.1)、NADH氧化酶/烯还原酶等，确认与已解析晶体结构同源。

3.2 Fold recognition tools suggest the presence of a natural grouping of flavoproteins in L. infantum

pGenThreader与I-TASSER折叠识别将利什曼黄素蛋白自然聚类为TR、DLD、NDH2、混合蛋白(deCoAR/nCoAR含NDH2亚基+小还原酶亚基)及小还原酶(吗啡酮还原酶等)；PyMOL结构叠合证实NDH2与deCoAR的NDH2亚基RMSD≤4 ?，人源无混合蛋白及小还原酶直系同源，利什曼DLD/TR与人GSR/TrxR/DLD相似度最高。

3.3 Sequence conservation of flavoproteins in L. infantum and other species

多序列比对显示TR与DLD辅因子(FAD、NAD(P)H)结合区残基保守性最高，NDH2次之，deCoAR更低；抑制剂结合区(距晶体配体5 ?内残基)保守性差异更大——TR与DLD可区分，NDH2介于TR/DLD与deCoAR间，小还原酶及α-磷酸甘油氧化酶明显偏离，被排除后续分析。

3.4 Analysis of TR structure, cofactors and inhibitor binding regions

LiTR(PDB 4apn)具4个辅因子结合区——FAD结合袋、NAD(P)H结合袋、锥虫硫酮(Trypanothione, GCG)结合袋(近二聚体界面)、泛醌(Ubiquinone)类似物袋；已知TR抑制剂结合于三个口袋：Site 1(近FAD但不重叠FAD袋，含JV0、RDS、MWT、H6H、TS8、WPF)、Site 2(远离催化口袋，含JV0*、RDS、TS8)、Site 3(近NADPH腺嘌呤部分，含RDS、BVN)。

3.5 Homology models of DLD, NDH2, and deCoAR

以人DLD(1zmd, 55%相同, 97%覆盖)模建Li DLD样；以酵母NDH2(4g6h, 28%相同, 81%覆盖)模建Li NDH2样；以E. coli deCoAR(1ps9, 44%相同, 91%覆盖)模建Li deCoAR样。叠合模板后RMSD分别0.641 ?、0.329 ?、1.467 ?(<2 ?)，模型质量可靠。

3.6 Redocking analysis on LiTR

FAD、NADPH重对接refRMS<2.5 ?验证格点盒合理性；氧化型GCG在单体可获合理pose(refRMS≈2.57 ?)，还原型GCG因柔性过大难复现(refRMS>4 ?)，改用T. cruzi TR(1bzl)作参比；部分大体积抑制剂(JV0 Site2、RDS Site1/Site2、BVN Site3)需在二聚体LiTR上重对接改善pose(RMSD<2.65 ?)，证明部分配体结合涉及单体—单体界面。

3.7 Preliminary selectivity analysis on human FAD/NAD(P)H-dependent dehydrogenases

对接人DLD(1zmd)、AIF(4bur)、GSR(1grb)、TrxR1(2j3n)：辅因子结合能相近；抑制剂中RDS、JV0、WPF在人DLD也获稳定pose提示潜在交叉反应；JV0、H6H、MWT、WPF预测ΔBE(人?TR)>0偏向LiTR选择性结合；TS8(Auranofin)、BVN ΔBE接近零或负提示可能与人黄素蛋白交叉反应。比较结合能差值(ΔBE = BE_human? BE_LiTR)可定性判断选择性趋势。

3.8 Selectivity and docking analysis on homology models of L. infantum FAD/NAD(P)H-dependent dehydrogenases structurally related to LiTR

Li NDH2模建中仅MWT(Site1)及泛醌类似物保留微摩尔级预测K_i，多数TR抑制剂不结合；Li DLD模建中WPF(Site1)、JV0(Site2)、RDS(Site2/3)有较好预测亲和力，提示部分TR配体可与寄生虫内源DLD交叉作用；Li deCoAR模建中RDS(Site3)预测K_i<500 μM。表明寄生虫黄素蛋白酶组内也存在脱靶可能，不同抑制剂选择性谱各异。

3.9 Docking of known FAD/NAD(P)H-dependent dehydrogenase inhibitors on TR

将氯代三联吡啶铂(Chloro(2,2′:6′,2″-terpyridine)platinum, TPT类似物)、卡莫司汀(Carmustine)、硝基呋喃腙(Nitrofurazone)、精胺(Spermidine, SPD)分别对接LiTR(推定位点)与CtNDH2。硝基呋喃腙可定位于LiTR GCG袋(refRMS合理)，精胺定位于GCG袋，氯代三联吡啶铂在CtNDH2泛醌袋附近得合理pose。

3.10 Expression and purification of LiTR and CtNDH2

重组LiTR(带N端6×His)与CtNDH2(带C端6×His)分别在大肠杆菌BL21(DE3)与C41(DE3)可溶性表达，Ni-NTA纯化后SDS-PAGE定量浓度分别为3.06 μg/μL与0.35 μg/μL。

3.11 Kinetic characterization of recombinant LiTR and CtNDH2

LiTR对锥虫硫酮表观K_m=48.55 μM，V_max=56.05 μmol·min^?1·mg^?1；CtNDH2对癸基泛醌(Decylubiquinone)表观K_m=90.65 μM，V_max=254.45 μmol·min^?1·mg^?1，符合Michaelis-Menten动力学。

3.12 Modulation of enzymatic activity by selected compounds

100 μM测试下：LiTR被Auranofin(残余活性5.76%)、氯代三联吡啶铂(4.44%)显著抑制，硝基呋喃腙抑制至41.52%，卡莫司汀(Carmustine)轻微激活(132.38%)，PX-12与精胺无影响；CtNDH2被卡莫司汀抑制至66.72%，氯代三联吡啶铂激活至186.05%，Auranofin/硝基呋喃腙/PX-12/精胺影响甚微。证实配体对同一黄素蛋白家族不同成员可产生抑制或激活的酶依赖性双向调控。

四、讨论与结论翻译

研究人员指出，L. infantum存在至少11种FAD/NAD(P)H依赖脱氢酶，其中DLD样、NDH2样及deCoAR样与TR结构相关但功能各异，其保守辅因子结合区可致某些TR抑制剂发生寄生虫内脱靶交叉反应，分子对接可揭示该倾向。Auranofin(已上市抗风湿药，具抗寄生虫报道)与硝基呋喃腙可抑制LiTR；氯代三联吡啶铂(Ⅱ)强力抑制LiTR同时激活CtNDH2，体现结构相关黄素酶局部残基差异引致配体相反功能后果。结合预测结合能与对接pose可靠性(refRMS<2.65 ?)可作交叉反应实验优先级判据。比较ΔBE提示JV0、H6H、MWT、WPF对LiTR较人同源酶具相对选择性，TS8/BVN潜在交叉反应需关注。本工作所建"序列鉴定→折叠识别→同源模建→辅因子/抑制剂重对接定义格点→跨酶对接比较"流程(图示Workflow)可用于FAD/NAD(P)H依赖脱氢酶家族配体选择性初筛，不仅服务抗利什曼药物开发，也可拓展至宿主氧化还原酶(如AIF、GSR、TrxR1、DLD)修饰剂设计及相关疾病(癌症、代谢病)研究。

Conclusive Remarks翻译：

总之，本研究确立了在FAD/NAD(P)H依赖脱氢酶更广家族内评估锥虫硫酮还原酶(Trypanothione Reductase, TR)抑制剂选择性的结构知情与实验验证框架。通过整合比较对接分析与重组LiTR和CtNDH2生化验证，研究表明结构相关黄素蛋白可在保守结合区容纳共同配体并产生抑制或酶特异性调节效应。预测结合能为负且对接pose可靠(refRMS<2.65 ?)的TR抑制剂应优先进入体内源及人源黄素蛋白酶交叉反应性测试，以预判脱靶并发现该氧化还原活性酶家族的多靶标调节可能。所提计算—生化流程不仅是优化抗利什曼药物发现策略的工具，也是调查结构保守黄素蛋白配体选择性的可迁移方法。Auranofin、Nitrofurazone及Terpyridine-derived Platinum化合物作为直接酶调节剂的实验验证强化了本研究转化前景，同时强调需对人源同源酶做系统性评价。鉴于FAD/NAD(P)H依赖脱氢酶的进化保守性，这些发现可超越寄生虫系统，为感染、代谢或氧化还原相关人类疾病涉及相关酶调节提供结构依据。未来需整合分子动力学模拟及扩展酶学与细胞学验证以精确界定治疗窗与转化适用性。