基于宏基因组挖掘与结构导向计算设计的新型耐热3-奎宁环酮还原酶ScQR的发现与催化表征

《Applied and Environmental Microbiology》：Genomic prospecting and biochemical characterization of a novel thermostable 3-quinuclidinone reductase from hot spring metagenomes for efficient biocatalysis

　　

本研究报道了通过温泉环境宏基因组挖掘（metagenomic mining）发现的一种新型嗜热3-奎宁环酮还原酶（3-quinuclidinone reductase，ScQR）的鉴定与表征。ScQR属于短链脱氢酶/还原酶（short-chain dehydrogenase/reductase，SDR）超家族，在大肠杆菌（Escherichia coli）中异源表达（heterologous expression），并对其催化性质进行了系统表征。该酶表现出卓越的热稳定性，在70°C条件下孵育48小时后仍保留86%的活性。此外，K⁺和Mg²⁺离子在特定浓度下显著增强了ScQR的活性。结构分析揭示ScQR采用典型的SDR折叠，具有保守的催化三联体（catalytic triad，S141-Y155-K159），且为NAD(H)依赖型。酶活性测定表明，ScQR对(R)-3-奎宁环醇（(R)-3-quinuclidinol）具有高度立体选择性（stereoselectivity），对其对映体(S)-3-奎宁环醇无催化活性。该酶在pH 9和85°C条件下表现出最优活性，使其成为需要高热稳定性工业应用的有前景的候选酶。分子动力学模拟（molecular dynamics simulations）进一步揭示，ScQR在高达360 K的温度下保持整体结构完整性，而更高温度则诱导结构失稳，主要集中在C末端区域和第95–100位残基。此外，通过LigandMPNN和UniKP实现的结构导向计算设计（structure-guided computational design）产生了三个ScQR变体，这些变体在维持整体折叠和功能的同时，表现出改善的底物亲和力（substrate affinity）和催化效率（catalytic efficiency）。该工作凸显了宏基因组学与结构驱动蛋白质设计相结合在从极端环境发现具有独特催化性质的新型酶方面的强大能力，并将ScQR确立为生物技术（biotechnological）和制药（pharmaceutical）应用中极具前景的生物催化剂（biocatalyst）。

**基于宏基因组挖掘与结构导向计算设计的新型耐热3-奎宁环酮还原酶ScQR的发现与催化表征**

**一、研究背景与意义**

自然界中绝大多数微生物在标准实验室条件下仍不可培养，这极大限制了人们对微生物功能基因资源的挖掘与利用。随着大规模测序技术的飞速发展与测序成本的显著降低，宏基因组学（metagenomics）已成为识别新型功能基因——特别是来源于极端环境古菌（archaea）的功能基因——的核心手段。这些极端微生物构成了一个尚未充分开发的遗传资源库，其编码的酶类往往具有卓越的稳定性和独特的催化特性，对合成生物学和生物制药开发具有重要推动作用。

短链脱氢酶/还原酶（short-chain dehydrogenase/reductase，SDR）超家族是一个高度多样化且在进化上保守的酶类群体，广泛分布于古菌、细菌和真核生物三大生命域中。SDR参与类固醇激素、脂质、维生素、糖类及碳源代谢等多种基础代谢过程。其中，3-奎宁环酮还原酶（3-quinuclidinone reductase，QR，EC 1.1.1.B51）催化3-奎宁环酮（3-quinuclidinone）立体选择性还原为(R)-3-奎宁环醇（(R)-3-quinuclidinol），后者是合成多种手性药物——包括瑞伐托酯（revatropate）、索利那新（solifenacin）和他沙利定（talsaclidine）——的关键手性中间体。尽管(R)-3-奎宁环醇的生物合成和代谢尚未在活体生物中被观察到，但已在多种原核和真核微生物中鉴定出3-奎宁环酮还原酶。尤其值得关注的是，来自高温环境（如温泉）古菌的SDR酶通常具有较高的热稳定性，这使其在需要高温条件的生物制药应用中具有不可替代的价值。

然而，SDR超家族成员间的序列同源性普遍较低（通常为20%–30%），即使具有相同化学和区域特异性的成员之间，序列一致性也可能低于20%。例如，来自同一菌株Clostridium sardiniense的7α-羟基类固醇脱氢酶（7α-HSDH）与7β-羟基类固醇脱氢酶（7β-HSDH）之间的序列一致性仅为17.8%，尽管二者共享相同的化学和区域特异性并受同一启动子调控。因此，仅通过序列比对进行功能预测极具挑战性，SDR成员的精确催化功能必须依赖广泛的实验筛选来最终确认。在此背景下，研究人员开展了此项研究，旨在通过宏基因组学策略从温泉环境中发现新型耐热3-奎宁环酮还原酶，并整合结构导向的计算设计方法对其进行工程化改造，以期为高温生物催化提供有力的酶学工具，并建立"宏基因组挖掘—结构解析—计算设计"的酶发现与工程化新范式。

**二、关键技术与方法概述**

在方法学层面，研究人员从NCBI序列读段归档（Sequence Read Archive，SRA）数据库获取温泉宏基因组测序数据（登录号SRR10095340），经组装后以已知QR（登录号OYR08778.1）为诱饵蛋白进行BLASTP分析，筛选序列相似性≥25%的开放阅读框（open reading frames，ORFs）构建候选文库，利用EFI-EST生成序列相似性网络并通过Cytoscape可视化。候选酶基因按大肠杆菌密码子偏好合成后克隆至pGEX-6p-1载体，转化大肠杆菌BL21进行异源表达，经GST亲和层析和PreScission蛋白酶柱上切除N端GST标签获得纯酶。酶活性通过340 nm处分光光度法测定还原型辅酶的吸光度变化。三维结构采用AlphaFold2深度学习框架预测。分子动力学（molecular dynamics，MD）模拟使用AmberTools24准备体系（ff19SB和gaff2力场，TIP3P水盒子），经能量最小化和多阶段升温后，在GROMACS2022中执行100 ns模拟，分析均方根偏差（root mean square deviation，RMSD）、均方根波动（root mean square fluctuation，RMSF）、回旋半径（radius of gyration，R_g）和溶剂可及表面积（solvent-accessible surface area，SASA）等参数。酶计算设计以结构为导向，采用LigandMPNN图神经网络对结合口袋及周边残基进行靶向优化，在保持天然三维结构的前提下增强底物识别能力，并以UniKP预测变体的米氏常数（K_m）、转换数（k_cat）和催化效率（k_cat/K_m）等动力学参数。

**三、主要研究结果**

**异源表达**：DNA测序确认ScQR正确插入pGEX-6p-1载体的多克隆位点（登录号OM654425.1）。ScQR由259个氨基酸组成，理论分子量28.14 kDa，注册于NCBI数据库。通过500 mL重组大肠杆菌BL21培养体系进行表达，经PreScission蛋白酶切除N端GST标签后，获得2.32 mg纯酶蛋白（浓度1.16 mg/mL），SDS-PAGE分析显示分子量与理论值一致。

**序列与结构比较**：多序列比对显示，ScQR与来自高温环境的另外两个QR及六个中温环境QR相比，序列一致性为27.5%–54.9%，但均保守地含有SDR超家族的特征基序，包括N端GIG基序、催化三联体（S141-Y155-K159）及辅酶偏好性位点。D36残基的存在提示ScQR为NAD(H)依赖型。AlphaFold2预测的三维结构呈现典型SDR架构，以七条平行β链形成的紧密β片层为核心，外围由α螺旋环绕，经环区连接构成紧凑的球状罗斯曼样折叠（Rossmann-like fold）。与其他SDR家族成员（PDB 4O0L、5EPO、5GT9）的结构比对TM分数为0.76–0.90，表明高度结构保守。

**酶活性**：研究人员对包括(R)-3-奎宁环醇、(S)-3-奎宁环醇、牛磺鹅去氧胆酸（taurochenodeoxycholic acid，TCDCA）、牛磺熊去氧胆酸（tauroursodeoxycholic acid，TUDCA）、葡萄糖、20-羟基蜕皮酮、乙醇和β-雌二醇在内的多种底物及NAD⁺/NADP⁺辅酶进行了系统筛选。结果表明ScQR的唯一底物为(R)-3-奎宁环醇，且表现出严格的立体选择性，对其对映体(S)-3-奎宁环醇无任何催化活性。实验同时确认该酶为NAD(H)依赖型。动力学参数为：K_m = 0.236 ± 0.004 mM，k_cat = 0.762 ± 0.016 s^?1，k_cat/K_m = 3.232 ± 0.087 s^?1 mM^?1。

**最适pH、温度和热稳定性**：ScQR在pH 9的Gly-NaOH缓冲体系中活性最高，最适温度约为85°C，表明其为典型的嗜热酶。热稳定性测试中，纯化ScQR在25°C孵育48小时后活性无任何变化；在70°C孵育48小时后仍保留86%的活性，展现出在SDR超家族中极为突出的热稳定性。

**离子对酶活性的影响**：研究人员考察了不同浓度金属离子盐酸盐（K⁺、Na⁺、Mg²⁺、Mn²⁺和Cu²⁺）及金属螯合剂乙二胺四乙酸（ethylenediaminetetraacetic acid，EDTA）对ScQR活性的影响。K⁺和Mg²⁺对酶活性具有最显著的增强作用，分别在300 mM和150 mM时达到约20%和50%的最优增幅。Na⁺和Mn²⁺在测试浓度范围内影响极小，而Cu²⁺和EDTA仅表现出抑制作用。EDTA可能通过螯合必需金属辅因子而抑制ScQR活性。

**酶失活机制**：通过100 ns分子动力学模拟系统分析了ScQR在不同温度下的分子结构稳定性。RMSD结果表明，酶结构在300 K、330 K和360 K下保持稳定，但当温度升至390 K时整体分子结构无法维持稳定性，结构偏差达2.667 ?。R_g和SASA分析也一致证实了390 K时的相对热不稳定性。RMSF分析进一步揭示，C末端区域（残基252–259）柔性和结构不稳定性显著增加，同时残基95–100区域的不稳定性也显著升高。综合这些结果，ScQR在高温下的结构失稳主要归因于羧基末端和第95–100位氨基酸片段的去稳定化。

**ScQR定向进化**：基于ScQR与底物的结合构象，研究人员运用LigandMPNN——一种在大规模蛋白质结构数据集上训练的图神经网络——实现了从目标结构到高亲和力序列的逆向设计。利用UniKP计算酶动力学参数进行筛选，获得三个高底物亲和力变体：redesign1（K_m = 0.084 mM，k_cat = 12.755 s^?1，k_cat/K_m = 151.845 s^?1 mM^?1）、redesign2（K_m = 0.079 mM，k_cat = 10.525 s^?1，k_cat/K_m = 133.228 s^?1 mM^?1）和redesign3（K_m = 0.066 mM，k_cat = 7.375 s^?1，k_cat/K_m = 111.742 s^?1 mM^?1）。三维模型结构比较显示三个变体与ScQR的RMSD分别为0.501、0.486和0.478 ?（均小于2 ?），表明整体构象高度保守。序列比对进一步揭示L193和L194在定向进化过程中发挥关键作用。该结构导向策略——利用LigandMPNN优化结合位点残基同时保持全局蛋白质稳定性——有效增强了底物亲和力和催化效率。

**四、讨论与结论**

讨论部分深入阐释了SDR超家族功能注释的固有挑战：尽管基因组和宏基因组分析能够发现大量SDR成员的开放阅读框，但由于家族成员间序列同源性普遍较低，基因功能预测无法直接确定底物类型，最终催化功能的确认必须依赖广泛的实验筛选和验证。古菌在极端温度、pH和压力环境中繁衍生息，其产生的酶类在这些条件下表现出卓越的结构和功能稳定性，具有重要的工业应用价值。随着测序成本的显著降低，宏基因组学已成为从极端环境中发现新型生物催化剂的强大工具。蛋白质工程研究虽已揭示蛋白质热稳定性的多重结构基础（如氢键网络、盐桥、疏水核心堆积等），但全面理解蛋白质热适应机制可能需要整合人工智能方法。机器学习技术则为新型酶的从头设计提供了前所未有的可能性。

在本研究中，分子动力学模拟精确揭示了ScQR在高温下C末端和残基95–100区域的结构破坏模式。基于此结构洞察，研究人员整合LigandMPNN和UniKP两种前沿深度学习框架，成功实现了ScQR的定向进化。LigandMPNN将蛋白质残基和配体表示为图模型，通过原子间距离和元素特征编码分子相互作用，精准捕获氢键和疏水效应等关键互作特征，从而实现对结合口袋残基的理性优化。值得注意的是，ScQR的热稳定性优于研究人员此前发现的其他SDR超家族耐热成员。

研究结论：该研究利用温泉环境宏基因组分析，通过同源比对、异源表达和特异性筛选，鉴定出一种新型耐热3-奎宁环酮还原酶ScQR。在70°C条件下孵育48小时后，ScQR仍保留86%的活性，彰显其卓越的热稳定性。此外，K⁺和Mg²⁺在特定浓度下显著增强了酶活性。ScQR的工业应用潜力包括在高温反应条件下参与瑞伐托酯、索利那新和他沙利定等手性药物的生物合成。更重要的是，该工作凸显了宏基因组学与结构驱动计算设计的协同整合作为一种强大且可推广的策略，可用于从极端环境中发现和工程化具有独特且具有工业价值的催化特性的新型生物催化剂，为生物催化与蛋白质工程领域提供了重要的方法论参考。