《Frontiers in Plant Science》：Structure-guided identification of catalytic determinants in ZmGA20ox enables rational CRISPR target selection for gibberellin engineering in maize
编辑推荐：
赤霉素20-氧化酶是控制赤霉素生物合成的关键酶，也是改造玉米植株结构和耐旱性的重要靶点。本研究开发了一种结构引导的计算工作流，整合基于人工智能的结构预测、分子对接、分子动力学模拟和体外CRISPR/Cas9引导RNA切割分析，研究了玉米赤霉素20-氧化酶3（Z
赤霉素20-氧化酶是控制赤霉素生物合成的关键酶，也是改造玉米植株结构和耐旱性的重要靶点。本研究开发了一种结构引导的计算工作流，整合基于人工智能的结构预测、分子对接、分子动力学模拟和体外CRISPR/Cas9引导RNA切割分析，研究了玉米赤霉素20-氧化酶3（ZmGA20ox3）中功能重要的残基。结构分析揭示了保守的Fe2?结合催化三联体（His147-Asp149-His166）被由Phe151、Trp157和Phe163组成的芳香族底物识别口袋所包围。对接分析表明，天然赤霉素底物和共底物2-氧戊二酸占据了预测的催化空腔，而调环酸在所检测配体中表现出最优的预测对接分数。计算突变分析显示，催化位点替换H147A、D149A和H166A的预测对接分数普遍较不利，而芳香族底物识别残基的替换则产生配体依赖性效应。随后的100 ns分子动力学模拟揭示了不同的配体依赖性构象行为，其中H147A影响预测催化环境的稳定性，而W157A主要影响底物结合口袋内的配体容纳能力。MM/GBSA分析进一步表明，有利的预测配体相互作用并不一定意味着具有催化能力的活性位点构型，突出了高效活性位点组织的重要性。最后，体外CRISPR/Cas9切割实验证明，靶向含有H147和W157的ZmGA20ox3区域的引导RNA具有序列特异性切割活性，支持在后续基因组编辑实验中靶向这些位点的可行性。总体而言，本研究为ZmGA20ox3中关键残基的潜在作用提供了结构引导的计算学见解，并建立了一个用于优先选择未来生化验证和基因组编辑候选残基的预测框架。
本文从结构生物学与计算生物学整合视角，系统解析了玉米（Zea mays L.）赤霉素20-氧化酶3（ZmGA20ox3）的催化机制，并建立了结构引导的CRISPR靶点优先化策略。研究背景方面，赤霉素（GA）是调控植物茎伸长、节间扩展、成花转变及生殖发育的关键二萜类激素，而GA20-氧化酶（GA20ox）是GA生物合成途径中的限速酶。在水稻与小麦中，通过半矮化育种（如水稻sd1基因、小麦Rht位点）已证实定向扰动GA稳态可显著改良株型与抗倒伏性；然而，玉米中精准设计的等位基因仍相对缺乏，且完全破坏GA生物合成常导致严重矮化、生殖适应性受损等不良多效性效应。因此，有必要在保留部分GA功能的前提下实现可调控的酶活性削减。ZmGA20ox3属于2-氧戊二酸依赖性双加氧酶（2-ODD）超家族，其催化核心依赖Fe2?和2-氧戊二酸（2OG）辅因子，是玉米中开展结构引导基因编辑的理想靶标。然而，作物GA20氧化酶的结构信息仍然匮乏，阻碍了对底物识别与催化机制的深入理解。为此，研究人员整合AlphaFold结构预测、活性位点表征、计算突变、分子对接、分子动力学（MD）模拟及体外CRISPR/Cas9切割实验，对ZmGA20ox3的催化体系进行了系统研究。
研究人员主要采用了以下关键技术方法。首先，从Ensembl Plants数据库获取ZmGA20ox3蛋白序列（UniProt A0A1D6N5E5），并使用AlphaFold数据库的预测结构（AF-A0A1D6N5E5-F1）作为模型；利用PROCHECK、ERRAT和Verify3D验证模型立体化学质量。其次，通过PyMOL Mutagenesis Wizard生成H147A、D149A、H166A、W157A、F151A、F163A六个丙氨酸替换突变体，并用GROMACS进行局部能量最小化。随后使用AutoDock Vina（经CB-Dock2平台）对GA12、GA53、GA20、2OG和调环酸（prohexadione）进行分子对接。然后选取野生型（WT）、H147A和W157A与2OG或GA20的六个复合物进行100 ns MD模拟（GROMACS，OPLS-AA/L力场），分析RMSD、Rg、RMSF、氢键及相互作用能。此外用gmx_MMPBSA计算结合自由能。最后，针对H147和W157编码区设计两条引导RNA（g1、g2），利用重组SpCas9进行体外切割实验。
研究结果分述如下。第一，ZmGA20ox3的结构特征：AlphaFold预测结构呈2-ODD家族典型的双链β-螺旋（DSBH）核心，中央存在溶剂可及的催化空腔；保守的Fe2?结合催化三联体His147-Asp149-His166构成催化中心，芳香族残基Phe151、Trp157、Phe163在周边形成疏水底物识别口袋。多序列比对显示这些残基在植物GA20氧化酶中高度保守。第二，结构验证：Ramachandran图显示90.2%残基位于最适区，ERRAT质量因子为89.78，Verify3D通过率91.9%；全局平均pLDDT为87.44，关键残基pLDDT均较高（W157为94.62，F163为96.75，H166为93.62），PAE值较低，表明结构可靠。第三，催化与底物识别残基鉴定：对接结果表明GA12、GA20、GA53及2OG均占据同一催化空腔，且邻近Fe2?中心；调环酸同样占据该空腔并与底物通道重叠，提示其竞争性抑制作用。第四，计算突变对接：丙氨酸替换后各突变体整体折叠保持，催化残基突变（H147A、D149A、H166A）主要破坏Fe2?配位环境，底物识别残基突变（W157A、F151A、F163A）则改变疏水口袋；对接分数呈配体依赖性变化。第五，蛋白-配体相互作用分析：WT与配体间形成氢键、疏水接触及Fe2?配位；催化残基突变减少氢键并导致配体位移增大，芳香族残基突变则减弱疏水稳定性。第六，MD模拟：WT–GA20复合物RMSD最稳定，H147A–GA20在约50 ns后构象波动显著增大，W157A居中；2O结合状态下H147A表现出更低RMSD和更致密结构。氢键分析显示H147A与2O形成持久氢键网络，但可能为非生产性结合；WT与GA20形成最持久氢键网络。相互作用能分析表明H147A与2O的相互作用能最负，而WT与GA20的相互作用能最有利。第七，MM/GBSA自由能计算：H147A–2O的ΔTOTAL最有利（?16.98 kcal/mol），其次为W157A–2O（?14.68）和WT–2O（?11.42）；而W157A–GA20的ΔTOTAL最有利（?17.50 kcal/mol），WT–GA20为?10.21，H147A–GA20最弱（?5.03）。这表明有利的结合自由能并不等同于催化活性，活性位点的生产性组织至关重要。第八，体外CRISPR/Cas9切割：设计的两条gRNA可引导SpCas9对含有H147和W157的基因组区域进行序列特异性切割，而无gRNA的对照无切割，证明这些位点在体外可被有效靶向。
讨论部分指出，催化三联体HXD…H基序对Fe2?配位和2OG氧化脱羧至关重要，破坏该基序会导致催化失活；而芳香族残基主要通过疏水堆积和π相互作用稳定底物方向。MD模拟中Fe2?未显式参数化，且仅对两个代表性突变体各生成一条100 ns轨迹，因此结果应视为构象动态的比较性证据而非直接催化结论。本研究建立的流程可将结构信息与功能预测整合到CRISPR靶点筛选中，弥补了传统仅依据核苷酸序列选择靶点的不足。结论部分综合认为，ZmGA20ox3的保守催化三联体（His147-Asp149-His166）是维持催化活性构象所必需，Trp157等芳香族残基主要介导底物识别与稳定；体外切割实验验证了H147和W157作为基因组编辑候选靶点的可行性。该工作提供了一种可推广的结构引导植物基因组编辑流水线，能够优先选择功能重要的编辑靶点，提高实验效率，并适用于缺乏实验结构的多种作物基因和代谢途径。
