来源：Binding pocket optimization of Pueraria mirifica type II chalcone isomerase synergizing with a necrosis-attenuating host transcription factor for persistent daidzein biosynthesis（Plant Physiology and Biochemistry）| 编译：DNA Lab Space | 原文许可：elsevier-subscription
导读
在非豆科植物底盘（本氏烟草）中重建豆科异黄酮合成通路，长期受制于异源酶催化效率低与宿主组织坏死两大瓶颈。本研究从葛根（Pueraria mirifica）中鉴定出查尔酮异构酶PmCHI II型，发现其对异甘草素具有独特的碱性pH偏好（pH 7.5–8.0）；结合500 ns分子动力学模拟与结构导向突变，获得I195V变体，体外异甘草素选择性比值达1.63。作者进一步用P2A核糖体跳跃肽构建多顺反子多基因盒，并以宿主来源的NbMYB60转录因子替换拟南芥MYB12，成功减轻组织坏死，将平台有效工作时间延长至9天，侵染后6天大豆苷元滴度达1.5 mg/g DW，个别高表现重复达2.5 mg/g DW。
研究背景
异黄酮是豆科植物特有的次生代谢物，参与结瘤共生与防御信号，又因结构与哺乳动物雌激素相似而备受生物医学关注。大豆苷元及其衍生物治疗潜力突出，是代谢解析的重点目标。葛根（White Kwao Krua）原产东南亚，是研究这类高价值异黄酮的重要天然模型与资源库。然而，野生群体中含量低、地理变异大、依赖季节，直接提取受到严重限制；化学合成则因立体专一性复杂、反应条件苛刻而难以成为可持续替代方案。用合成生物学平台重建这些通路，既可绕开天然代谢约束，也能提供可系统优化、可解析通路的可控模型。异黄酮核心从头合成以苯丙烷途径提供的对香豆酰辅酶A与丙二酰辅酶A为起点，经查尔酮合酶（CHS）缩合生成查尔酮，再由查尔酮异构酶（CHI）催化立体专一性环化，进入下游黄酮与异黄酮骨架。CHI因此被视为通路的关键调控闸门。豆科植物在演化中通过基因复制与新功能化产生了II型CHI，能够接受脱氧查尔酮异甘草素，这是豆科异黄酮分支得以延伸的分子基础。本研究即以葛根CHI为切入点，试图同时解决酶效率与宿主坏死两个瓶颈。
研究方法
2.1 植物材料与培养
葛根块根采自泰国Phichit农业研究与发展中心（16.4362° N，100.2844° E）并完成鉴定。切片块根经液氮速冻后保存于−80 °C。本氏烟草种子按Methods S1所述，在25 ± 2 °C、16 h光照/8 h黑暗光周期下培养。取四週龄植株用于瞬时表达试验。
【数据】温度：25 ± 2 °C；光周期：16 h光照/8 h黑暗；植株龄期：4周；保存温度：−80 °C；采集地坐标：16.4362° N，100.2844° E
2.2 序列注释与系统发育分析
对葛根转录组（登录号SRR6917866）进行加工、注释与翻译。系统发育树用MEGA X构建，蛋白基序鉴定用MEME套件完成，参数详见Methods S2。
【数据】转录组登录号：SRR6917866；系统发育软件：MEGA X；基序鉴定：MEME suite；参数：原文未详述（见Methods S2）
2.3 分子建模与分子动力学模拟
依据各自氨基酸序列，用SWISS-MODEL生成PmCHI I型与II型的同源模型，模板分别为PDB ID 6MS8（I型）与1EYP（II型）。将模型蛋白叠合到与异甘草素结合的南极发草CHI1晶体结构（PDB ID: 5YX4）上构建柚皮素查尔酮结合复合物，随后手动为配体添加6′-羟基。PmCHI II型–异甘草素复合物由同一结构参考生成，不作配体修饰。蛋白与配体质子化状态分别用PDB2PQR与ChemAxon指定。每个复合物在AMBER 24中按既定方案于NPT系综下进行500 ns分子动力学模拟。详细模拟参数与突变体建模流程见Methods S3。
【数据】模拟时长：500 ns；软件：AMBER 24；系综：NPT；模板PDB：6MS8（I型）、1EYP（II型）、5YX4（配体参考）；质子化工具：PDB2PQR、ChemAxon
2.4 分子克隆与载体构建
用序列特异性引物（表S3、S4）从葛根块根提取总RNA、合成cDNA，并完成天然与突变PmCHI及通路酶片段的组装。重组片段分别克隆入pET21b与pEAQ-HT-DEST1载体，用于大肠杆菌与本氏烟草表达。亚细胞定位载体通过将PmCHI与GFP/tGUS-GFP标签融合构建，并转化农杆菌GV3101。完整克隆试剂盒、载体骨架与组装参数见Methods S4。
【数据】载体：pET21b、pEAQ-HT-DEST1；菌株：农杆菌GV3101；标签：GFP、tGUS-GFP；引物：表S3、S4；参数：原文未详述（见Methods S4）
2.5 细菌表达与酶纯化
工程大肠杆菌菌落于LB培养基中培养，用IPTG诱导重组His标签PmCHI蛋白表达。收集细胞经冻融循环与超声裂解，目标蛋白用镍柱纯化。纯化组分经SDS-PAGE验证，浓缩、脱盐后于−80 °C保存。详细培养参数、裂解缓冲液与洗脱曲线见Methods S5。
【数据】培养基：LB；诱导剂：IPTG；标签：His；纯化介质：镍柱；保存温度：−80 °C；验证方法：SDS-PAGE
2.6 体外底物特异性与底物竞争实验
标准酶促反应使用纯化PmCHI变体与相应底物，在HEPES缓冲液（pH 7.5）中进行。野生型酶的底物特异性以柚皮素查尔酮、异甘草素、紫铆因与卡达莫宁评估；底物竞争实验使用柚皮素查尔酮与异甘草素的二元混合物。反应经HPLC分析。完整缓冲液组成、反应时间与提取参数见Methods S6。
【数据】缓冲液：HEPES；pH：7.5；底物：柚皮素查尔酮、异甘草素、紫铆因、卡达莫宁；分析：HPLC；反应时间：原文未详述（见Methods S6）
2.7 pH优化与浓度依赖性活性谱
在25 °C下用96孔板读数仪进行分光光度测定，监测385 nm吸光度下降。在此条件下，用覆盖pH 5.0–11.0的多缓冲液体系确定最适pH谱。浓度依赖性活性谱中，取线性期初始速度（ν0），在GraphPad Prism 8.0中按Michaelis–Menten方程非线性回归拟合，估算Km与kcat，详见Methods S7。
【数据】温度：25 °C；检测波长：385 nm；pH范围：5.0–11.0；拟合软件：GraphPad Prism 8.0；参数：Km、kcat（原文未给出具体数值）
2.8 农杆菌介导的本氏烟草瞬时表达与代谢物提取
携带相应重组载体的农杆菌GV3101菌株经培养后，重悬于优化侵染缓冲液，每个构建体终OD600调至0.3。用无针注射器将悬液侵染四週龄本氏烟草叶片背面。处理叶片于侵染后6天（dpi）采收，速冻、冻干，经超声提取及酸水解后分析。详细侵染缓冲液组成、孵育方案与代谢物提取步骤见Methods S8。
【数据】菌株：农杆菌GV3101；OD600：0.3/构建体；植株龄期：4周；采收时间：6 dpi；提取：超声+酸水解
2.9 LC-MS/MS靶向蛋白定量
叶片蛋白提取后，用Amicon Ultra-0.5过滤器（10 kDa MWCO）脱盐，胰蛋白酶消化。肽段在Agilent 1290 Infinity II LC与6545XT Q-TOF质谱仪上分析。MaxQuant 2.6.3针对目标数据库处理原始数据，三重复间采用Sample Minimum Imputation（SampMin），见Methods S9。
【数据】过滤器截留：10 kDa MWCO；酶：胰蛋白酶；仪器：Agilent 1290 Infinity II LC + 6545XT Q-TOF；软件：MaxQuant 2.6.3；插补方法：SampMin；重复数：3
2.10 基于UHPLC–qTOF/MS的异黄酮定量
植物化学提取物在Shimadzu qTOF 9030系统上以负ESI模式、Luna C18(2)柱（Phenomenex，美国）进行轮廓分析与定量。目标（异）黄酮经多点校准曲线定量，原始轮廓用MS-DIAL 5.5提取（Methods S10）。
【数据】仪器：Shimadzu qTOF 9030；离子模式：负ESI；色谱柱：Luna C18(2)；定量方式：多点校准曲线；软件：MS-DIAL 5.5
2.11 酶活分析的HPLC-DAD检测
assay底物与产物在Ultimate-3000 HPLC系统（Thermo Scientific，美国）上、30 °C、Kinetex C18柱（Phenomenex，美国）定量。于285 nm对照标准曲线定量，Chromeleon处理数据（Methods S11）。
【数据】仪器：Ultimate-3000 HPLC；柱温：30 °C；色谱柱：Kinetex C18；检测波长：285 nm；软件：Chromeleon
2.12 亚细胞定位分析
于3 dpi从本氏烟草取圆形叶碟，立即用ZEISS LSM 980共聚焦显微镜配合ZEN 3.13软件观察。完整激发/发射参数、检测窗口与封片剂组成见Methods S12。
【数据】取样时间：3 dpi；显微镜：ZEISS LSM 980；软件：ZEN 3.13；参数：原文未详述（见Methods S12）
2.13 电解质渗漏实验
用打孔器从侵染植株取叶碟。每个重复将叶碟浸入去离子水，用电导率仪（Mettler Toledo，美国）测定初始电导率（EC0）。详细取样间隔、孵育参数与计算流程见Methods S13。
【数据】仪器：Mettler Toledo电导率仪；初始指标：EC0；取样间隔与孵育参数：原文未详述（见Methods S13）
研究结果
3.1 葛根查尔酮异构酶的鉴定与系统发育分类
对葛根转录组挖掘，鉴定出三个推定的查尔酮异构酶编码基因（PmCHIs）。为确定其演化归属，进行了最大似然系统发育分析（图S1A）。结果显示三个候选基因分属不同功能分支：一个PmCHI基因聚在普遍存在的CHI I型分支内，第二个基因归入豆科特有的CHI II型群（图S1A），第三个unigene被归为非催化型CHI IV型增强子。结构序列比对进一步表明，PmCHI II型活性裂隙内具有特异的苏氨酸-甲硫氨酸（TM）残基组合，而I型CHI为丝氨酸-异亮氨酸对（图S1B）。基于这些结构-化学特征，PmCHI II型被确定为下游通路重建与后续结构导向理性设计的候选酶。
【数据】候选基因数：3个；分支归属：1个I型、1个II型、1个IV型；II型活性位点残基：Thr-Met（TM）；I型对应残基：Ser-Ile（SI）
3.2 野生型PmCHI同工酶的酶学表征、pH谱与亚细胞定位
重组野生型PmCHI I型与II型采用低温诱导策略表达纯化以保持可溶性。针对四种查尔酮评估底物特异性：柚皮素查尔酮、异甘草素、紫铆因与卡达莫宁。HPLC图谱显示两种同工酶功能明显分化。对柚皮素查尔酮，两种同工酶均具催化能力，但仅缓冲液对照也出现大量自发环化（图2A）。相反，异甘草素在缓冲液中保持稳定，功能分离清晰：I型对该底物无催化活性，II型则高效催化其转化为甘草素（图2B）。对结构类似物紫铆因，I型无活性，II型高效催化其转化为紫铆因苷元（butin）（图2C）。对卡达莫宁，II型产物形成可忽略，I型亦无活性（图2D）。
连续分光光度法追踪385 nm底物吸光度显示，纯化的I型与II型酶均呈陡峭的钟形pH–活性曲线（图2E）。对柚皮素查尔酮，两种同工酶最大催化速度均位于7.2–7.5的中性pH窗口。对豆科特异底物异甘草素，野生型PmCHI II型最大催化速率略向碱性范围偏移至7.5–8.0（图2E）。因此选择pH 7.5的缓冲体系作为后续所有活性实验的操作基线。为确定细胞微环境，用拟南芥组蛋白H2B-mCherry核标记在本氏烟草瞬时共表达监测体内亚细胞定位（图2F）。共聚焦激光扫描显微镜显示，野生型PmCHI I型与II型均分布于胞质与细胞核（图2F）。为检验入核是否由被动扩散驱动，作者构建了大型大分子串联标签：将PmCHI II型与截短β-葡萄糖醛酸酶（tGUS）及eGFP融合，使总分子量升至约117 kDa。共聚焦成像显示，该117 kDa融合变体的荧光信号严格局限于细胞质，完全排除于核外（图2F）。
【数据】底物数：4种；检测波长：385 nm；柚皮素查尔酮最适pH：7.2–7.5；异甘草素最适pH：7.5–8.0；操作pH：7.5；融合标签分子量：约117 kDa；核标记：H2B-mCherry；观察时间：3 dpi
3.3 分子动力学模拟与计算机突变筛选
为阐明野生型PmCHI I型与II型底物特异性差异的分子决定因素，对三个复合物进行500 ns MD模拟：PmCHI I型–柚皮素查尔酮、PmCHI II型–柚皮素查尔酮、PmCHI II型–异甘草素（Methods S3）。RMSD分析表明三个复合物均快速达到平衡并在整个模拟中保持结构稳定。三个复合物的平均蛋白–配体氢键数分别为3 ± 1、1 ± 1与1 ± 1。原子接触分析同样显示相互作用随时间持续，三个复合物的平均接触数分别为204 ± 20、240 ± 21与198 ± 21。因此，尽管两个PmCHI II型复合物的氢键数略少于I型复合物，三者均维持稳定的蛋白–配体结合。
为表征平衡态结合模式，提取每条轨迹最后50 ns进行详细蛋白–配体相互作用分析。用MM/PB(GB)SA方法进行逐残基自由能分解，以≤ −0.5 kcal/mol为能量截断标准识别底物识别关键残基（图3A）。对PmCHI I型–柚皮素查尔酮复合物，若干疏水残基（F50、T51、I53、V98、Y109、V113、S193、M194与V200）对配体结合有有利贡献，其中S193与V200稳定作用最强，提示C端区域可能参与柚皮素查尔酮在结合口袋内的容纳。L41、I43、L104、Q108与K112的额外贡献通过预测的范德华与极性接触进一步稳定配体。相反，PmCHI II型中柚皮素查尔酮结合主要由R39、K112、V113、N116与M194稳定，而II型对应位置193与200（T193与V200）贡献可忽略，尽管其I型对应位置（S193与V200）是该同工酶中最强的两个贡献者。这表明，尽管催化架构保守，两种同工酶招募的残基网络不同。对PmCHI II型–异甘草素复合物，主要有利能量贡献来自R39、L41、F50、T51、L104、E108、Y109、K112、V113、N116与M194，其中V113稳定效应最大，其次为R39与Y109。该谱提示异甘草素主要通过其芳香核心周围的疏水堆积稳定，并辅以羰基与羟基的极性相互作用。总体而言，PmCHI II型的底物识别似乎由以R39、K112、V113、N116与M194为中心的局部网络支配，与PmCHI I型明显不同。能量筛选突出显示超过11个具有有利分解自由能（≤−0.5 kcal/mol）的关键残基，支持特定序列分歧可调节活性位点几何构型的观点。在这些热点残基中，R39因是II型查尔酮异构酶中高度保守、催化必需的残基而被排除在后续计算突变筛选之外，仅保留在实验突变面板中作为机制对照以验证其催化必需性。为评估突变其余活性位点残基的热力学影响，选取异甘草素结合口袋周围七个候选位置（T51、L104、E108、K112、V113、N116与M194）进行计算机突变筛选（图3B），并相对野生型PmCHI II型评估相对结合自由能差（ΔΔG bind）。在T51位，芳香族替换T51Y与T51F相对野生型给出不利的较高ΔΔG bind值，而不带电极性替换T51N与T51Q给出更负的有利值。在L104位，L104F的ΔΔG bind明显比野生型更负，而位于新……（原文节选至此结束）
【数据】模拟时长：500 ns；复合物数：3；氢键数：3 ± 1（I型–柚皮素查尔酮）、1 ± 1（II型–柚皮素查尔酮）、1 ± 1（II型–异甘草素）；接触数：204 ± 20、240 ± 21、198 ± 21；能量截断：≤ −0.5 kcal/mol；分析轨迹：最后50 ns；筛选位置：7个（T51、L104、E108、K112、V113、N116、M194）；热点残基数：>11
讨论与解读
CHI的功能分化是植物特化代谢的主要调控闸门，决定中间查尔酮向下游黄酮核心结构的立体专一性环化方向。本研究中PmCHIs的系统发育分群直接反映葛根在豆科中的分类位置，并与其积累多种异黄酮的代谢能力相符。从演化角度看，普遍存在的I型CHI是陆生植物基础催化模板，依靠保守氢键网络与活性位点坐标稳定配体中间体；特化的II型分支则在豆科中经基因复制与新功能化产生，以满足不同的特化代谢需求。转录组挖掘还鉴定出属于非催化IV型演化谱系的第三个unigene，其功能为辅助增强子而非直接生物催化剂。结构分析强调，PmCHI II型的双底物能力——即在经典柚皮素查尔酮之外获得接受脱氧查尔酮异甘草素的能力——由精确的演化开关驱动（原文节选至此结束）。
编译者解读：该研究把"酶口袋改造"与"宿主工程"两条线拧成一股：I195V只是保守替换，体外选择性比值1.63并不惊艳，且定量蛋白组显示体内酶丰度无显著差异，说明增益来自催化偏好而非表达量——这对合成生物学是重要提醒，酶工程收益常被体内表达波动掩盖。用P2A多顺反子简化多基因堆叠、用宿主源NbMYB60替代AtMYB12缓解坏死，是更实用的底盘层面贡献。局限在于提升倍数批间波动大，I195V的体内机制仍待更严格验证。
参考来源
期刊：Plant Physiology and Biochemistry, 2026
DOI: 10.1016/j.plaphy.2026.111883
DOI: 10.1016/j.plaphy.2026.111883