IF50.5｜Rubisco单点突变生化景观图谱

大规模平行筛选揭示Rubisco单点突变对CO2亲和力的影响，少数突变可显著提升亲和力。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

Rubisco是地球上最主要的CO2固定酶，但其催化速率缓慢，且易发生与O2的副反应，被认为是光合作用的关键限速步骤。长期以来，研究者推测Rubisco的生化参数（如催化速率kcat、CO2亲和力KC、CO2/O2特异性SC/O）之间存在权衡，限制了其进化与工程改造。尽管已有数十年研究积累了几百个天然及突变Rubisco的生化数据，但主要集中在维管植物，细菌Rubisco的多样性远未充分采样。近年来，高通量筛选与机器学习为蛋白质工程提供了新工具，但如何在高通量下解析单个催化参数仍是挑战。

科学问题

卡在哪里： 此前对Rubisco的工程改造大多失败，部分原因是缺乏高通量、精确测量其生化参数的方法。传统体外生化测量通量低，难以系统探索序列空间。此外，Rubisco的序列-功能关系尚不明确，尤其是哪些突变能够改善CO2亲和力，以及这种改善是否以牺牲催化速率为代价。

本文要回答： 本文旨在建立一种高通量筛选方法，系统绘制Rubisco单氨基酸突变体的序列-功能图谱，并从中鉴定出能够改善CO2亲和力的突变，为后续酶工程提供基础。

技术路线

作者主要用了：深度突变扫描、大肠杆菌生长筛选、CO2浓度滴定、Michaelis-Menten拟合、体外生化验证。

作者利用工程化大肠杆菌Δrpi菌株，该菌株在甘油为唯一碳源时依赖Rubisco活性才能生长。他们构建了覆盖Rhodospirillum rubrum Rubisco（Form II）99%以上单点突变的条形码文库，转化入Δrpi菌株，在不同CO2浓度下进行竞争性生长筛选。通过比较筛选前后条形码丰度计算每个突变体的适应度，并利用Michaelis-Menten方程拟合不同CO2浓度下的适应度数据，推断每个突变体的表观最大速率和CO2半饱和常数。随后对部分突变体进行体外纯化和动力学测定以验证筛选结果。

图 1　Fig. 1。图1展示了筛选系统原理和整体数据质量。a-b显示Δrpi菌株生长依赖于Rubisco通量；c为文库筛选流程；d显示生物学重复间高度一致；e显示适应度与文献kcat值线性相关；f显示适应度双峰分布；g为全突变体适应度热图，顶部为序列保守性，底部为活性位点残基。

核心亮点

亮点 1｜筛选覆盖99%单点突变

作者构建了包含8,760个单氨基酸突变体（占设计文库99%以上）的条形码文库，在Δrpi菌株中进行了9次重复筛选，平均Pearson相关系数达0.98。突变体适应度呈双峰分布，约72%的突变有害，仅0.14%的突变适应度高于野生型。适应度与文献报道的77个kcat值及35个体外测量值呈线性相关，验证了筛选的可靠性。

图 2　Fig. 2。图2将适应度映射到Rubisco结构上。a显示同源二聚体表面着色，活性位点标记星号；b比较不同结构区域的突变效应；c为活性位点和Loop 6放大图，雷达图显示各位置突变效应；d为位置平均适应度与系统发育保守性的散点图，突出显示异常值如G186和M215-H257相互作用。

亮点 2｜结构区域突变敏感性差异

通过将适应度映射到Rubisco同源二聚体结构上，发现溶剂暴露区域对突变更耐受，而活性位点和埋藏残基高度敏感。值得注意的是，Loop 6区域尽管参与催化，但部分残基（如E331、E333）高度耐受突变，而活性位点K329则极为敏感。此外，序列保守性与平均适应度呈负相关，但存在多个高度保守却耐受突变的位置（如G186）。

图 3　Fig. 3。图3展示CO2滴定实验设计和结果。a为滴定示意图；b为不同CO2浓度下突变体适应度变化；c-d为两个突变体的适应度随CO2浓度变化及Michaelis-Menten拟合；e-f为体外动力学测量验证；g为所有突变体的KC估计分布；h显示Vmax和KC估计的反相关。

亮点 3｜CO2滴定推断亲和力

通过在不同CO2浓度下进行文库筛选，作者拟合了每个突变体的表观Vmax和KC。65%的突变体（5,687个）获得了可靠的KC估计。体外验证显示，部分突变体的实际KC低于预测值，表明筛选可能低估了亲和力改善。四个突变（A102Y、V266T、A289C、A289T）表现出显著降低的KC，其中V266T和A102Y的体外KC分别约为80 μM和更低，远低于野生型的150 μM。

图 4　Fig. 4。图4聚焦于亲和力改善突变。a为KC与效应大小的散点图，突出四个低KC突变；b显示这些突变在结构上的位置，位于二聚体界面；c为体外kcat和KC测量，显示这些突变落在植物Rubisco分布边缘。

亮点 4｜亲和力改善以催化速率为代价

对V266T和A102Y的体外动力学分析表明，它们提高CO2亲和力的同时，kcat成比例下降，但CO2/O2特异性SC/O未受影响。这些突变位于同源二聚体界面附近，而非活性位点，暗示其可能通过改变活性位点构象或电子性质影响CO2结合。这些突变体的kcat和KC值已超出此前所有细菌Form II Rubisco的范围，接近植物和藻类Rubisco的分布边缘。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： 深度突变扫描、大肠杆菌生长筛选、CO2浓度滴定、Michaelis-Menten拟合、体外生化验证
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事酶工程、合成生物学或光合作用研究的科研人员，尤其是对高通量功能筛选和蛋白质进化感兴趣者。

带走一句

高通量生长筛选结合CO2滴定可同时估计大量突变体的酶动力学参数，为解析序列-功能关系提供了高效框架，并揭示了Rubisco亲和力改善的可行路径。