文章标题： Connecting single-cell transcriptomes to projectomes in the mouse visual cortex
中文标题： 破译神经元远程线路：从单细胞转录组连接全脑投射图
关键词： 单细胞转录组、Patch-seq、神经元投射组、多模态细胞分型、视觉皮层
这篇文章通过整合单细胞转录组、电生理、局部形态和全脑轴突投射数据，探索了小鼠视觉皮层兴奋性神经元的分子身份、功能特征与远距离连接模式之间的对应关系，这对于建立真正具有生物学解释力的神经元分类体系、理解大脑皮层信息传递规律具有重要意义。
近年来，单细胞RNA测序已经建立了覆盖全脑的大规模转录组细胞分类图谱，但某种基因表达类型的神经元究竟投射到哪里仍然是脑科学中的关键空白。Patch-seq能够在同一细胞中同时获得转录组、电生理和局部形态信息，却通常无法完整保留长距离轴突；全神经元形态成像能够显示一条神经元在全脑范围内的轴突走向，却缺少对应的转录组信息。因此，作者将两类互补数据连接起来：一类是1,528个视觉皮层兴奋性Patch-seq神经元，另一类是341个具有完整轴突重建信息的全神经元形态样本。
作者首先根据Patch-seq数据定义了17种形态—电生理—转录组类型，即MET类型。随后利用Patch-seq神经元与全神经元形态数据共有的树突特征，构建分层分类器，将没有转录组信息的全神经元形态样本映射到相应的MET类型。由此，研究实现了从基因表达身份到完整投射路线的跨模态连接。
研究发现，转录组差异并不是与神经元形态、电活动和连接模式彼此独立的。相反，无论是在不同MET类型之间，还是在同一类型内部，基因表达的连续变化都能够对应树突结构、动作电位模式及投射靶区的变化。细胞在视觉皮层中的解剖位置与转录组潜在因子相结合，还可以预测其是否投射到丘脑、高级视觉区、纹状体或上丘等特定脑区。
文章尤其揭示了第五层端脑内投射神经元和端脑外投射神经元的复杂异质性。例如，大多数端脑内投射类型将主要轴突保留在同侧初级视觉皮层，而L4/L5 IT及L5/L6 IT Car3等类型表现出更强的跨半球投射，并在对侧皮层中形成不同的层特异性终止模式。L5 ET-1 Chrna6神经元则在转录组、形态、电生理和投射数量上均表现出明显独特性。
这项工作的重要贡献不只是增加了若干新的神经元类型，而是提出了一套可扩展的跨数据集、跨模态细胞身份映射框架。它表明，神经元类型不能只依赖某一个标志基因或某一种形态来定义，而应同时考虑分子组成、电生理状态、局部结构、解剖位置和全脑连接。未来，这一方法可扩展到其他脑区、其他物种以及突触连接数据，为建立全脑多模态细胞图谱提供基础。
研究使用出生后45—70天的转基因小鼠，重点采集初级视觉皮层及高级视觉区的兴奋性神经元。Patch-seq实验首先在急性脑片中对单个神经元进行全细胞膜片钳记录，给予标准化的超极化和去极化电流刺激，提取动作电位、输入阻抗、膜时间常数、放电适应等电生理特征；随后吸取细胞核和胞质进行单细胞RNA测序，并用生物胞素填充细胞，对树突和局部轴突进行三维重建。
研究人员将Patch-seq细胞映射到既有的转录组T-type参考体系，再结合局部形态和电生理特征，将相关T-type整合为17个MET类型。研究使用降维、聚类和交叉验证分析不同模态间的一致性，并通过稀疏降秩回归提取能够同时解释转录组变化与形态、电生理变化的潜在因子。
在平行的全神经元形态实验中，作者通过转基因标记、病毒稀疏标记和全脑荧光成像获得单个神经元的完整树突与长距离轴突结构，并将其注册到Allen小鼠脑公共坐标框架。图1所示的核心流程先依据丘脑、脑干和皮层投射模式，将全神经元划分为IT、ET、CT等广义投射类别，再利用树突深度、分支范围、顶树突形态和局部轴突特征训练亚类特异性分类器，把341个全神经元映射到MET类型。
最后，作者使用逻辑回归模型，将MET类型、胞体位置、树突特征及转录组潜在因子作为输入，预测单个神经元是否投射到特定脑区，并分析各类特征对预测结果的相对贡献。
Patch-seq数据最终支持17种兴奋性神经元MET类型。不同类型在皮层层位、顶树突结构、放电形式和转录组空间中形成相对稳定但并非完全割裂的分布，说明神经元身份兼具离散类别与连续变化两种属性。稀疏降秩回归进一步显示，转录组潜在因子能够解释多项电生理和形态特征，证明同一类型内部的分子差异也可能对应功能差异。
形态分类器将341个全神经元样本连接到MET体系后，研究获得了不同转录组类型的完整投射图谱。多数L2/3、L4、L5 IT-2和L6 IT神经元约有80%—90%的轴突保留在同侧初级视觉皮层。相比之下，L4/L5 IT和L5/L6 IT Car3神经元超过一半的轴突离开同侧初级视觉皮层，尤其偏向对侧皮层，但二者在对侧终止层位上不同：前者更偏向第一层及第四、第五层交界，后者更多进入第二、第三层和上部第四层。
高级视觉区神经元与初级视觉皮层同类型神经元具有相似的靶区选择，但高级视觉区神经元局部轴突较少、远距离投射较多，提示细胞所处区域会重新调整局部加工与跨区域通信之间的平衡。
L5 ET-1 Chrna6类型在全部模态中均表现出独特性。这类细胞较少出现爆发式放电，投射靶区数量也少于其他L5 ET类型，说明传统的“端脑外投射神经元”类别内部仍包含功能不同的回路单元。模型还发现，细胞在视觉皮层中的位置对高级视觉区和尾状壳核等靶区的预测较重要，而转录组潜在因子对外侧膝状体和外侧后核等靶区更有解释力。例如，Rxfp1表达可能帮助识别投射到外侧后核的部分L5 ET神经元。
首次系统连接转录组与完整投射组。 通过树突形态这一共享信息桥梁，把Patch-seq和全脑单神经元形态数据真正整合起来。
建立多模态神经元分类体系。 不再把标志基因、形态或投射方向中的任何一项视为唯一分类依据。
揭示同一转录组类型内部的连续功能差异。 证明转录组潜在变化可对应树突、电生理及投射靶区变化。
具有较强扩展性。 该框架原则上可应用于全脑、其他物种及突触连接数据。
全神经元形态样本仅341个，部分MET类型样本数量较少或缺失，可能影响稀有类型投射规律的稳定性。
全神经元样本的MET身份是通过形态分类器间接推断，而不是在同一细胞中直接获得转录组与完整轴突。
研究主要集中于成年小鼠视觉皮层兴奋性神经元，尚不能直接推广到抑制性神经元、其他脑区或人类大脑。
投射存在并不等于形成有效突触，研究没有直接测量突触强度、信息方向及神经元在行为中的因果作用。
文章标题： Competing programs shape cortical sensorimotor–association axis development
中文标题： 大脑皮层如何分区？两套竞争程序塑造感觉—联想轴
关键词： 皮层发育、感觉运动—联合轴、丘脑皮层投射、MIND模型、轴突导向
这篇文章通过整合人类、猕猴和小鼠的发育转录组、单细胞数据、空间表达图谱、遗传操作与神经回路示踪，探索了大脑皮层感觉运动区与高级联合区如何在发育过程中形成相互分离又彼此联系的空间结构，这对于理解人类高级认知网络的起源及神经发育障碍具有重要意义。
大脑皮层可以沿感觉运动—联合轴组织：一端是视觉、听觉、躯体感觉和运动等初级区域，另一端是支持抽象认知、社会行为和跨模态信息整合的前额叶、顶叶和颞叶联合网络。过去的研究较清楚地解释了初级感觉皮层如何受到丘脑输入和感觉活动的塑造，却难以解释联合皮层为何具有广泛分布、渐变边界和长距离连接。
作者提出网络发育中的多节点诱导—排斥模型，简称MIND模型。该模型认为，大脑皮层的宏观分区不是由单一分子梯度决定，而是由两套空间方向相反的发育程序竞争产生。第一套是周中央程序，从额叶和颞叶两端附近产生，逐渐向皮层几何中心推进，赋予神经元高级联合皮层特征；第二套是中央程序，在未来的初级感觉和运动区域形成，主要由第一阶感觉运动丘脑皮层输入诱导，并主动排除周中央程序。
多物种转录组分析表明，联合皮层相关基因模块首先在额颞两极及海马、杏仁核等区域出现，随后沿额—颞方向向内扩展；感觉运动相关模块则在未来初级区域形成多个局部中心。随着胎儿发育推进，两类模块逐渐形成显著负相关，最终对应成年大脑的感觉运动—联合轴。
两套程序涉及的分子功能也不同。联合程序富集轴突导向、细胞黏附、视黄酸信号、WNT信号、跨突触通信及边缘系统发育基因；感觉运动程序则富集动作电位、长时程增强、SMAD信号和神经突生长相关基因。多种自闭症风险基因也分布于这些模块，提示皮层空间组织异常可能成为神经发育疾病的共同基础。
作者进一步利用减少丘脑皮层输入或改变皮层投射神经元命运的遗传模型检验MIND模型。当第一阶丘脑皮层输入减少，中央感觉运动程序的范围收缩，联合皮层分子特征和连接模式则向原本属于初级区域的空间扩展。由此说明，初级感觉运动区不仅被丘脑输入诱导出来，同时还通过占据空间和建立排斥性边界，限制联合程序继续扩张。
文章还确定了两种代表性分子：中央程序中的SEMA7A和周中央程序中的PLXNC1。二者呈互补空间分布，并参与初级区域与高级联合区域轴突之间的排斥。相关基因敲除后，感觉运动—联合轴上的皮层连接出现异常。作者因此将初级区域比喻为分布式联合网络“海洋”中的若干“岛屿”，认为诱导与排斥是共同建立皮层层级结构的两项基本原则。
研究首先整合三套独立的灵长类发育脑转录组数据，包括人类脑组织微阵列、人类RNA测序和猕猴RNA测序数据。作者选择四个未来初级区域和七个联合皮层区域，并加入海马和杏仁核，依据胎儿发育阶段中的区域共表达关系构建基因模块。联合模块被划分为额叶和颞叶两部分，感觉运动模块则对应运动、躯体感觉、视觉和听觉四个中心。只有在三套数据中重复出现的基因才被保留为严格的共享模块。
作者利用基因集富集分析识别各模块涉及的信号通路，并结合胎儿人脑单细胞RNA测序数据，判断这些基因在兴奋性投射神经元、抑制性神经元及不同皮层层次中的表达。随后使用原位杂交、RNAscope和免疫染色等方法，在小鼠及其他物种组织中验证候选基因的空间分布。
为了检验丘脑输入的作用，研究分析了Celsr3、Gbx2等条件性敲除模型，这些模型具有不同程度的丘脑皮层轴突减少；同时使用Fezf2、Satb2、Zbtb18等皮层发育相关遗传操作，使初级感觉运动区域发生收缩或身份改变。研究比较SEMA7A、PLXNC1及其他区域标志物的表达边界，并通过顺行和逆行示踪评估联合皮层连接是否向初级区域扩展。
最后，作者构建Sema7a和Plxnc1单独或联合缺失的小鼠，在初级躯体感觉皮层和内侧前额叶注射不同颜色的逆行AAV，对沿感觉运动—联合轴的投射神经元进行定量，从分子和回路层面检验两种信号之间的排斥作用。
多套人类和猕猴数据均显示，联合皮层模块与感觉运动模块在胎儿中晚期逐渐形成两个方向相反的表达梯度。额叶联合程序先在前部达到高表达，再向顶叶和颞叶延伸；颞叶联合程序则从颞叶向内扩展。感觉运动模块随后在未来的运动、躯体感觉、视觉和听觉皮层形成局部表达中心。到发育后期和成年期，两类模块呈显著负相关，证明它们不是简单并行表达，而是存在空间竞争关系。
功能富集结果显示，联合程序偏向轴突导向、细胞黏附、视黄酸、WNT和突触形成，感觉运动程序则偏向神经元兴奋性、长时程增强和感觉回路成熟。单细胞分析表明，这些程序并不限于单一细胞类型，而是贯穿多类兴奋性和抑制性神经元，但具有明显的区域和层次差异。
在丘脑皮层输入减少的Celsr3及Gbx2相关模型中，中央SEMA7A表达区域明显缩小，原本受限于联合皮层的PLXNC1及其他周中央特征向感觉运动区扩展。改变皮层投射神经元命运、导致初级区域缩小的遗传模型也出现类似结果：感觉运动岛屿收缩，联合皮层连接和分子标志占据更多皮层空间。这些结果支持中央程序依赖丘脑输入诱导，并对联合程序产生排除作用。
SEMA7A和PLXNC1在正常小鼠中沿感觉运动—联合轴形成互补梯度。联合敲除相关信号后，初级躯体感觉皮层与前额叶之间的投射分布发生改变，原本受到限制的轴突进入不适宜区域，说明这一配体—受体系统参与建立不同层级皮层网络之间的边界。
提出MIND统一模型。 用多节点诱导—排斥同时解释初级皮层、联合皮层及层级网络的形成。
证据跨越物种和技术层级。 从人类胎儿转录组到小鼠遗传操作，再到神经回路示踪，形成较完整的证据链。
强调联合皮层具有主动发育程序。 联合区不再被视为初级区域之间被动剩余的空间。
发现SEMA7A—PLXNC1轴突排斥机制。 将宏观皮层分区与具体分子和连接机制联系起来。
MIND模型主要由多种静态发育阶段数据重建，尚未直接实时观察两套程序在同一个体中的推进和竞争过程。
不同物种的数据平台、取样区域和发育阶段并不完全一致，跨物种对应关系仍存在一定推断成分。
多种遗传操作可能产生细胞命运、轴突生长或整体脑发育等多重效应，不能完全归因于单一诱导—排斥过程。
自闭症风险基因的富集主要属于关联性证据，尚未证明这些基因通过MIND机制直接导致具体临床表型。
文章标题： Chromatin landscape and epigenetic heterogeneity of acute myeloid leukaemia
中文标题： 越过基因突变：染色质图谱重构急性髓系白血病分型
关键词： 急性髓系白血病、染色质可及性、ATAC-seq、表观遗传异质性、精准治疗
这篇文章通过对1,563名初诊急性髓系白血病患者进行大规模ATAC-seq和多组学分析，探索了染色质可及性如何定义白血病细胞状态、临床预后和药物敏感性，这对于突破仅依赖基因突变的AML分型体系、发展表观遗传指导的精准治疗具有重要意义。
现有AML分类主要依据融合基因、驱动突变和染色体异常，但具有相同遗传改变的患者在细胞分化状态、疾病进展和治疗反应方面仍可能截然不同。作者由此提出，遗传异常并不能完整解释AML异质性，而染色质开放状态可能记录了驱动突变、细胞起源、分化阶段和环境影响共同形成的最终调控结果。
研究建立了由瑞典和日本患者组成的eCHROMA队列，在全部1,563例样本中进行染色质可及性分析，并整合驱动基因测序、拷贝数变化、RNA测序、全基因组测序、DNA甲基化、组蛋白修饰、转录因子网络、单细胞RNA/ATAC和离体药物筛选等数据。研究从不同患者中识别出176,853个重复出现的开放染色质区域，其中大量差异最明显的位点位于非启动子区域，提示增强子等远端调控元件是AML表观遗传异质性的主要来源。
根据染色质可及性，作者将AML划分为16个ATAC亚群。部分亚群与PML::RARA、RUNX1::RUNX1T1、CBFB::MYH11和CEBPA等经典遗传亚型高度对应，但多数亚群无法由单一突变或突变组合准确界定。同一个遗传亚型可能分布到多个ATAC亚群，而同一ATAC亚群也可能包含不同遗传背景。这表明，染色质状态不是基因突变的简单重复，而是一个独立的疾病分层维度。
不同ATAC亚群具有特定的造血分化阶段、基因表达、DNA甲基化、超增强子和转录因子活性。染色质开放程度与DNA甲基化总体呈反向关系，亚群特异性超增强子则连接了关键转录因子和白血病细胞身份。单细胞多组学进一步表明，一个患者内部虽然存在处于不同分化位置的白血病细胞，但这些细胞通常共享所属亚群的核心染色质指纹。这意味着ATAC亚群不是由少数正常细胞污染造成，而是白血病克隆广泛具有的稳定调控状态。
在临床层面，ATAC亚群与白细胞计数、原始细胞比例、分化状态和生存结局密切相关。即使将WHO、ICC或其他基因组分类纳入模型，ATAC分型仍然具有独立预后价值。离体药物筛选也显示，不同ATAC亚群对特定靶向药、信号通路抑制剂和细胞毒药物存在不同敏感性。
因此，这项研究将AML分型从有哪些突变推进到这些突变最终建立了怎样的调控网络。染色质图谱有望作为基因组信息的补充，用于识别传统分型中被混合在一起的患者、预测疾病行为，并为药物选择提供新的生物标志物。
研究纳入1,563名初诊AML患者，样本来自瑞典和日本多个临床中心。所有患者均接受ATAC-seq，以测定全基因组开放染色质区域；同时使用覆盖331个已知或候选驱动基因的靶向捕获测序，并结合SNP位点分析拷贝数异常。部分样本进一步接受RNA-seq、全基因组测序、DNA甲基化芯片、ChIP-seq以及单细胞RNA/ATAC联合分析。
作者对ATAC峰进行统一质控和跨队列整合，共获得176,853个重复开放区域。随后依据患者间差异较大的可及性峰开展无监督聚类，建立16个ATAC亚群。通过重抽样和共聚类分析评价稳定性，平均调整兰德指数为0.79，共聚类准确率约为0.97，但O和P亚群稳定性相对较低。
研究使用正常造血细胞参考图谱对bulk ATAC数据进行去卷积，推断各亚群对应的造血分化阶段；再将亚群与WHO、ICC遗传分型、FAB形态分类、驱动突变、融合基因、拷贝数变化及临床参数比较。
ChIP-seq用于检测H3K27ac、H3K27me3、CTCF、SMC1和RNA聚合酶Ⅱ等信号，以识别增强子、超增强子和调控边界。转录因子基序富集和基因调控网络分析用于寻找各亚群的关键调控因子。单细胞RNA/ATAC数据用于重建白血病分化轨迹，并检验同一患者不同细胞是否共享亚群特异性染色质状态。最后，研究通过多变量生存分析和离体药物筛选评估ATAC亚群的独立预后价值与药物反应差异。
ATAC-seq将AML稳定划分为16个亚群，各亚群在开放染色质、DNA甲基化、基因表达、造血分化状态和超增强子组成方面均具有独特特征。非启动子开放区域的患者间变异明显高于启动子区域，说明增强子重塑是AML异质性的核心来源。
四个亚群可由经典融合基因或CEBPA特定突变较好定义，但其余亚群无法通过已知驱动突变或突变组合完全重建。即使进一步分析全基因组测序数据，也未找到能够单独解释多数ATAC亚群的新遗传事件。ATAC分型与遗传分型因此既相互联系，又提供互补信息。
与传统基因亚型相比，ATAC亚群与白血病细胞分化程度、白细胞计数、原始细胞比例和表达状态的相关性更强。某些具有相似NPM1或KMT2A相关改变的病例，会因染色质状态和分化轨迹不同而进入不同亚群，反映其调控网络和临床行为存在实质差异。
单细胞分析发现，同一患者内不同成熟程度的白血病细胞仍保留一致的亚群特异性染色质指纹；与此同时，染色质可及性和基因表达沿异常分化轨迹协同变化。超增强子分析进一步指出，各亚群由不同造血关键转录因子网络维持，而亚群特异性超增强子是稳定白血病身份的重要调控枢纽。
临床分析显示，ATAC亚群可独立于现有基因组风险分类预测患者结局。药物筛选则发现不同亚群呈现可重复的药物敏感性模式，说明染色质状态可能用于识别潜在治疗弱点。不过，这些结果目前主要来自离体实验，仍需前瞻性临床验证。
队列规模大。 在1,563例AML患者中系统开展ATAC-seq，显著超过以往小样本表观基因组研究。
建立16类染色质分型。 揭示传统遗传分类无法解释的表观遗传异质性。
多组学证据完整。 同时整合突变、表达、甲基化、组蛋白标记、超增强子和单细胞轨迹。
具有临床转化潜力。 ATAC亚群具有独立预后价值，并与离体药物敏感性相关。
RNA-seq、甲基化、ChIP-seq及单细胞数据只覆盖部分患者，不同分析层面的有效样本量并不一致。
队列主要来自日本和瑞典，尚需在更多族群、儿童AML和不同医疗体系中验证。
研究集中于初诊样本，未系统追踪化疗、复发和克隆演化过程中ATAC亚群是否发生转变。
O、P等部分亚群稳定性较低；临床应用还需要解决样本处理、测序成本、分类阈值及跨实验室标准化问题。
药物敏感性主要来自离体实验，尚不能直接等同于真实患者的临床疗效。
文章标题： Targeted enzyme discovery using metal-coordination mining
中文标题： 从金属配位到新酶：结构挖掘发现广谱自由基卤化酶
关键词： 酶发现、金属配位、AlphaFold、自由基卤化酶、生物催化
这篇文章通过将金属酶催化机理转化为原子尺度的三维配位特征，在AlphaFold蛋白结构数据库中定向挖掘稀有自由基卤化酶，探索了如何摆脱对序列相似性和偶然发现的依赖，这对于高效发现工业生物催化剂及理解蛋白质序列—结构—功能关系具有重要意义。
传统酶挖掘主要依赖序列同源性，但同一蛋白超家族可能共享相似折叠，却催化完全不同的反应；反过来，功能相同的酶也可能因长期进化而具有很低的序列相似性。随着AlphaFold等结构预测数据库快速扩大，研究重点开始从“寻找相似序列”转向“寻找相似活性位点”。然而，仅比较整体结构仍难以区分同一折叠中的羟化酶、去饱和酶、差向异构酶和卤化酶。
作者选择Fe²⁺/α-酮戊二酸依赖酶作为验证对象。这类酶均形成高活性的Fe(IV)-氧物种并从非活化C(sp³)—H键夺取氢原子，但后续反应结果取决于铁中心的配位环境。多数羟化酶由两个组氨酸和一个天冬氨酸或谷氨酸配位铁；自由基卤化酶则缺少第三个酸性配体，从而为空间中的氯离子保留配位位置。也就是说，2His加开放配位位点本身就是从反应机理推导出的卤化酶结构指纹。
作者将两个组氨酸之间的空间距离、所在β链的骨架氢键关系及第三配位残基类型编码为最小三维搜索模体。该方法不需要对所有蛋白进行两两结构比对，因此计算成本较低。研究从超过2.2亿个AlphaFold蛋白条目中逐级筛选含cupin结构域的蛋白，定位潜在2His金属位点，再根据第三配位位置是Asp/Glu还是Gly/Ala，区分常见非卤化酶与候选卤化酶。
筛选最终获得946个高优先级自由基卤化酶候选。序列相似性网络不仅回收了已知自由基卤化酶，还识别出约70个此前没有卤化酶注释的候选簇。作为对比，过去约20年通过天然产物研究仅发现少数主要家族，显示原子级结构挖掘能够显著拓展可探索的进化空间。
作者选择一个新候选簇进行实验验证，发现AspX和BtnX两种新酶。AspX能够对游离L-天冬氨酸进行立体选择性氯化，生成3S-氯-L-天冬氨酸。BtnX则以生物素为天然底物，在特定位置生成2R-氯代生物素，并能够使用溴离子和叠氮离子形成相应产物。
更引人关注的是，BtnX可接受多种具有羧酸头基的底物，表现出自由基卤化酶中罕见的底物宽容性。1.20 Å高分辨率晶体结构显示，生物素的羧酸基通过Gln119、Ser120和Arg248附近的氢键被精确定位，而其余部分位于开放、疏水且与溶剂连通的通道中。这种反应位点严格识别、底物其余部分宽松容纳的结构模式解释了BtnX的广谱性。
研究证明，酶发现可以从化学机理出发，在预测结构海洋中寻找决定功能的最小原子特征。该框架还可扩展到其他金属配位模体，甚至用于识别参与分子识别和信号转导的非酶蛋白功能位点。
作者首先从Fe²⁺/α-酮戊二酸依赖酶的反应机理出发，定义自由基卤化酶所需的最小结构模体。该模体包括位于相邻β链上的两个组氨酸、约束二者骨架位置的氢键以及适合铁配位的组氨酸侧链距离。第三配位位置若为Asp或Glu，则更可能属于羟化、去饱和或其他反应酶；若为Gly或Ala并留下开放位置，则被归为卤化酶候选。
研究从AlphaFold蛋白结构数据库检索超过2.2亿个蛋白，利用InterPro等注释提取约180万个含cupin结构域的蛋白，并进一步过滤低质量片段、重复序列和不完整结构。算法定位潜在2His金属位点并分类其配位环境，最终形成候选卤化酶集合。
作者利用序列相似性网络分析候选蛋白的进化聚类，并结合基因邻域判断其可能底物。例如，邻近酰基载体蛋白基因提示酶可能作用于载体结合底物，而邻近氨基酸转运、生物素转运或代谢基因则可提示自由小分子底物。
候选AspX和BtnX经异源表达和纯化后，在Fe²⁺、α-酮戊二酸、氧气、卤离子和抗坏血酸条件下进行体外酶反应。研究使用液相色谱—质谱、同位素峰和产物结构分析确认氯化、溴化及叠氮化产物，并测定稳态动力学和底物范围。最后通过X射线晶体学解析BtnX与铁、氯离子、α-酮戊二酸和生物素的复合物结构，并利用定点突变检验关键氢键残基。
金属配位挖掘在复杂cupin超家族中准确回收了全部主要已知自由基卤化酶类型，并从低序列相似度区域发现大量新候选。筛选获得946个高置信候选，其序列网络包含约70个此前没有被标记为自由基卤化酶的簇，覆盖细菌、真菌和植物等广泛系统。对于传统同源搜索只能找到少数成员的稀有真核卤化酶，结构方法可扩展出包含数十个候选的完整簇。
实验验证确认AspX能够氯化游离L-天冬氨酸，形成3S构型产物，代表一类作用于α-羧酸底物的新自由基卤化酶。BtnX则将生物素转化为2R-氯代生物素，其动力学参数为k_cat约0.96 min⁻¹、K_m低于2 μM，显示对生物素具有较高亲和力。BtnX还能够转移Br⁻和N₃⁻，生成溴代及叠氮化衍生物。
底物筛选显示，BtnX可以接受多种不同长度和结构的羧酸底物，底物宽容性显著高于多数已知自由基卤化酶。晶体结构显示，底物羧酸头基通过特定氢键被固定，而主体位于开放通道中。Ser120及周围氢键网络突变会同时降低活性和化学选择性，验证了该结构模型。BtnX处在载体蛋白依赖卤化酶与作用于游离小分子卤化酶之间的进化分支点，可能是一种更容易通过定向进化改造的酶骨架。
从反应机理而非序列出发。 将开放金属配位位点转化为可计算搜索的三维功能指纹。
计算效率高。 不需要大规模两两结构比对，适合持续扩张的蛋白结构数据库。
显著拓展酶家族空间。 一次性发现约70个此前未被识别的候选卤化酶家族。
完成实验闭环。 从数据库筛选到AspX、BtnX生化验证、动力学和晶体结构解析。
发现高可塑性BtnX。 兼具宽底物范围以及氯化、溴化和叠氮化能力，具有生物催化潜力。
946个候选中只有少量代表得到实验验证，其余候选仍可能包含错误预测或无活性蛋白。
方法依赖AlphaFold静态结构，难以处理活性位点需要大幅构象变化、辅助蛋白或翻译后修饰才能形成的酶。
当前框架适用于具有明确几何配位特征的功能，对缺乏清晰原子模体的酶类未必有效。
基因邻域只能推测生理底物和生态功能，AspX、BtnX在天然生物体中的真实作用仍需验证。
BtnX的催化周转率距离工业生产要求仍有差距，需要通过定向进化改善稳定性、活性及底物选择性。
致谢橙子牛奶糖（陈文燕），请用参考模版：We thank the blogger (orange_milk_sugar, Wenyan Chen) for XXX
感谢小可爱们多年来的陪伴， 我与你们一起成长~
此内容由惯性聚合(RSS阅读器)自动聚合整理，仅供阅读参考。 原文来自 — 版权归原作者所有。