Rieske二加氧酶（ROs）的理性工程改造：基于结构解析拓展毒性芳香族化合物生物降解的底物谱与催化效率

《FEBS Letters》：Rational engineering of Rieske dioxygenases for biodegradation of toxic aromatic compounds

　　

毒性芳香族化合物日益加剧的污染、持久性与难降解性，亟需能够优于天然微生物代谢的高效生物催化解决方案。Rieske（2Fe-2S）型二加氧酶（Rieske (2Fe-2S)-type dioxygenases, ROs）是一类有前景的酶，已知可通过羟化芳香环降解多种毒性芳香族化合物，从而促进其开环。ROs包括联苯二加氧酶（biphenyl dioxygenase, BPDO）、邻苯二甲酸二加氧酶（phthalate dioxygenase, PDO）、对苯二甲酸二加氧酶（terephthalate dioxygenase, TPDO）和间苯二甲酸二加氧酶（isophthalate dioxygenase, IPDO），这些酶将各自底物降解为相应的顺式-二氢二醇（cis-dihydrodiols）。尽管其具有较窄的底物特异性限制，但可通过重新设计活性位点结构（active-site architecture）和底物进入通道（substrate-access tunnels）对其进行定制，以生成具有更宽底物范围和更高催化效率的变体。IPDO的理性工程改造已证明，引入靶向突变可将底物范围扩展至多种邻苯二甲酸区域异构体（phthalate regioisomers），且不损害催化性能。类似地，BPDO的突变研究表明，替换相互作用残基可增加构象柔性，并容纳结构多样的底物，如多氯联苯（polychlorinated biphenyls, PCBs）。本综述基于结构研究分析多种Rieske二加氧酶的理性工程改造，以拓展其用于毒性芳香族污染物生物降解的催化范围。

引言
生物修复利用微生物及其氧化还原酶、水解酶和转移酶降解污染物，但天然途径受底物范围窄、酶效率低和环境适应性差限制；天然酶多针对天然化合物而非合成异源物。早期改造依赖定向进化和随机突变，近期结构机制研究促进基于结构的理性设计。Rieske二加氧酶（Rieske dioxygenases, ROs）为含[2Fe-2S]中心的非血红素铁加氧酶，可催化毒性芳香族化合物二羟基化；其改造靶点包括活性位点、亚基或结构域、底物进入通道和柔性环，结构功能数据不足限制非常规热点识别。

Rieske二加氧酶的特征
ROs将O₂两个氧原子掺入芳环形成顺式-二氢二醇（cis-dihydrodiols），促进开环。系统由还原酶（含FAD或FMN，接受NAD(P)H电子）、铁氧还蛋白（含[2Fe-2S]，传递电子）及加氧酶（含Rieske [2Fe-2S]和单核非血红素Fe(II)中心）组成。家族催化机器保守，功能多样性主要来自活性位点和进出通道；底物特异性由口袋形状、大小和静电性质决定。PDO、TPDO与IPDO比较显示活性位点残基细微变化即可改变底物识别，提示可改造可塑性。

Rieske二加氧酶的基本结构
Rieske型[2Fe-2S]由两个Cys和两个His配位，Fe1由Cys硫配位，Fe2由His配位，主要参与电子传递而不结合底物或氧。单核Fe(II)由两个His、一个Glu/Asp及一个水分子配位，参与氧活化、底物氧化和产物形成。加氧酶可呈异源六聚体α₃β₃、同源六聚体α₃α₃或α₆及同源三聚体α₃。异源六聚体中Rieske中心和2His-1Glu/Asp三联体位置保守；电子自Rieske簇Fe1经配位Cys/His及水网络传至三联体Asp，再传至Fe(II)。

Rieske二加氧酶的催化机制
底物进入靠近Fe(II)的疏水口袋，由催化三联体和水分子稳定。NAD(P)H电子经还原酶FAD、铁氧还蛋白[2Fe-2S]、加氧酶Rieske簇[2Fe-2S]传至Fe(II)；Rieske簇由[2Fe-2S]²⁺还原为[2Fe-2S]¹⁺，O₂结合后Fe(II)转为Fe(III)。连续电子传递和质子化生成活性中间体，攻击芳环并形成顺式-二氢二醇；追加电子并释放产物后Fe(II)再生。该机制高度保守，为改造底物特异性提供基础。

理性酶工程中的多种方法与新兴工具
ROs催化限速首步，催化中心保守，底物结合残基微小修饰即可拓展底物范围。理性酶工程整合结构、生化与计算分析，较定向进化和随机突变更具靶向性；结构生物学、计算建模、人工智能（artificial intelligence, AI）和合成生物学可提升活性、稳定性和特异性。基于X射线晶体学、冷冻电镜或核磁共振（nuclear magnetic resonance, NMR）的结构，以及同源建模和AlphaFold、RoseTTAFold、SWISS-Model、Modeller等预测结构，可识别活性位点、结合口袋、进入通道和出口环。多序列比对有助于定位保守催化残基和金属结合基序。活性位点工程针对底物相互作用残基、氢键和静电环境；通道工程调控底物进入、产物释放和溶剂可及性，计算工具可识别限制性空间屏障。

Rieske二加氧酶的理性设计
活性位点突变常显著影响效率，但也可能破坏功能，因此需关注活性位点外的柔性环、通道和门控元件。催化腔尺寸决定可结合底物范围，改变腔壁残基可拓展底物谱。近年已针对多氯联苯（polychlorinated biphenyls, PCBs）、DDT、多环芳烃（polyaromatic hydrocarbons, PAHs）、PET和增塑剂等开展ROs工程。

改造联苯二加氧酶（BPDO）以降解PCB同系物
BPDO决定PCB同系物是否进入分解代谢。BphAELB400中Thr335、Phe336和Asn337影响区域特异性，其中Thr335和Asn337不直接接触底物但调节周转。BphAEp4中Thr335Ala或Phe336Met可容纳2,6-二氯联苯；Thr335Ala还赋予Gln322、Gly321和Val320构象自由。BPDORR41的Asn338Gln、Ile341Val和Leu409Phe可降解2-氯二苯并呋喃。BphAEII9以BPDOB356的333GINTIRT339替换335TFNNIRI341，形成更柔性的入口环，对多种氯代联苯及DDT活性优于亲本。

改造NarA2B2用于PAHs生物降解
高分子量PAHs水溶性低、毒性高。嗜热菌NarA2B2可降解低分子量PAHs和芘但效率低。基于AlphaFold2结构预测及与NarAaAb比较，鉴定Val236、Tyr300、Trp316和Leu375造成空间阻碍；替换后获得Trp316Ile、Tyr300Phe-Trp316Ile及Val236Ala-Trp316Ile-Leu375Phe等变体，在保留低分子量PAHs降解能力的同时，对高分子量PAHs效率提高2—4倍，并扩展至荧蒽、苯并[a]蒽、苯并[b]荧蒽和苯并[a]芘等。

改造间苯二甲酸二加氧酶（IPDO）以降解邻苯二甲酸区域异构体
IPDO催化间苯二甲酸生成间苯二甲酸-3,4-顺式-二氢二醇。晶体结构显示底物C3/C4靠近催化铁，Arg234和His257与羧基形成盐桥，Val178、Phe249和Ile294产生空间限制。与PDO和TPDO比较，IPDO缺乏保守羧基结合残基；PDO的Ser182和Arg207在IPDO中对应Val178和Leu198，极性/带电残基变为非极性残基可解释其不能降解邻苯二甲酸。理性构建的IPDO Val178Ala/Phe249His双突变体可羟化三种邻苯二甲酸区域异构体，且保留天然活性；Val178Ala减小空间阻碍，Phe249His引入极性羧基相互作用。

萘1,2-二加氧酶：改造区域选择性和立体选择性
NDO可催化萘生成(+)-cis-(1R,2S)-1,2-二羟基-1,2-二氢萘，并具催化混杂性。围绕单核铁的Asn201、Phe202、Val260、Trp316、Thr351、Phe352、Trp358和Met366的突变显示Phe352是产物选择性关键；Phe352Val和Phe352Thr可逆转联苯cis-3,4-二氢二醇的对映选择性，主要生成(3S,4R)而非野生型偏好的(3R,4S)。晶体学分析表明该残基大小和几何改变底物取向，从而重定向区域和立体选择性。

甲苯二加氧酶：改造底物范围和选择性
TDO是芳香族不对称氧官能团化及手性顺式-二氢二醇合成平台，但天然活性位点限制大体积底物。半理性工程构建覆盖活性位点和进入通道27个位点的单/双突变库，筛选176个变体；Met220、Ala223和Phe366影响活性及化学、区域和对映选择性。优化变体可生成(+)-(R)-1,2,3,4-四氢喹啉-4-醇（71%分离产率、94% ee）及cis-二羟基化2-苯基吡啶衍生物（60%分离产率、98% ee），并提升对酯基官能化底物的活性。

咔唑1,9a-二加氧酶：改造底物取向和反应选择性
CARDO可催化咔唑和二苯并-p-二噁英的角式二加氧、蒽的侧向二加氧及芴和二苯并噻吩的单加氧。围绕Ile262、Phe275、Gln282和Phe329构建17个定点变体；Ile262Val和Gln282Tyr增强咔唑侧向二加氧，Phe275Trp促进4-羟基芴形成，Phe275Ala改变荧蒽氧化，Ile262Ala/Ile262Trp改变二苯并噻吩氧化，而Phe329突变多使活性丧失。结构分析显示咔唑在Ile262Val中旋转约15°，在Phe275Trp和Gln282Tyr中约25°，表明疏水腔微小改变可重定位底物并切换氧化路径。

对甲苯磺酸甲酯单加氧酶（TsaM）：单加氧与二加氧的理性切换
野生型TsaM偏好对甲苯磺酸和4-甲基苯甲酸甲位取代基的单加氧。Met230和Thr232位于活性位点顶部，控制底物与单核铁的距离和取向；引入Met230Phe/Trp、Thr232Phe/Trp及Met230Trp/Thr232Ile等大体积残基可转向芳环二加氧并专一形成顺式-二氢二醇；Met230Gly/Thr232Gly增加柔性则促进连续单加氧。由此提出“分子尺”模型，即底物与催化铁距离决定单加氧、连续单加氧或二加氧。

结论与展望
ROs具有保守的[2Fe-2S]Rieske中心、Fe(II)催化中心及电子传递通道。改造柔性环有望提高环境耐受性，而调控催化腔结构是拓展底物范围的主要因素。关键残基微小变化即可改变活性；应靶向决定口袋理化与结构性质的残基，并考虑柔性环和二级结构对口袋尺寸及底物进入的调节。定点突变是基本手段，基于结构识别热点后结合随机突变亦适用于柔性口袋。分子动力学模拟和通道分析可揭示动态与通道柔性对底物选择的影响。ROs多组分体系的长程结构耦合及与铁氧还蛋白/还原酶的电子传递尚不清楚，突变可能远距离影响催化效率；其对NAD(P)H的依赖增加反应复杂性和成本，未来需发展辅因子再生或替代电子供体。机器学习和AI辅助设计可预测突变效应并加速定向进化，但现有工具多预测稳定性、折叠或口袋几何，难以直接预测催化效率和底物特异性；生成模型仍需实验验证。对ROs这类多组分体系，应将计算预测与实验结构、靶向突变和动力学验证结合。