AlphaFold3预测人类血红素结合蛋白结构的评估

《International Journal of Molecular Sciences》：Evaluation of AlphaFold3 for Predicting Human Heme-Binding Protein Structures

　　

血红素是一种关键的辅因子，参与多种生物学过程，如电子传递、催化和氧结合。理解血红素与血红素结合蛋白（HBPs）的结合对于阐明其功能和分子机制，以及用于酶工程和疗法开发至关重要。AlphaFold3（AF3）是一种基于人工智能（AI）的大分子预测工具，已应用于包括HBPs在内的多种蛋白质的结构建模。然而，AF3是否为HBPs提供可靠的结构信息尚未得到研究。为了确定AF3对HBPs的预测效果，研究人员检查了AF3生成的四种人类HBP（包括亚硫酸盐氧化酶的细胞色素b₅结构域（SO-b₅）、NADH细胞色素b₅氧化还原酶的细胞色素b₅结构域（Ncb5or-b₅）、细胞色素b₅类型B（CYB5B）和神经红蛋白（NGB））的脱辅基（apo）和全蛋白（holo）状态，并将其与实验HBP结构进行比较。血红素分子的总体位置被很好地对接到了HBPs的血红素结合口袋中；然而，在血红素结合构型上存在差异，包括口袋几何形状和配位环境，这些对功能解释至关重要。实验NGB结构在血红素结合区域附近含有一个二硫键，而AF3预测模型缺少这个键，导致局部折叠的差异，从而影响血红素结合环境。分子动力学模拟表明，与实验结构相比，AF3预测的NGB结构表现出不同的分子构象和柔韧性。这些数据表明了使用AF3预测结构模拟HBPs血红素结合状态的潜力和局限性。

以下是根据原文内容撰写的论文解读文章：

**研究背景与问题**

血红素（heme）是一种含铁的卟啉化合物，作为关键辅因子参与氧气运输、电子传递和氧化还原反应等众多生命过程。血红素结合蛋白（Heme-binding proteins, HBPs）在生命科学、医学、环境技术和工业中应用广泛，例如细胞色素P450（CYP）参与药物代谢，肌红蛋白和血红蛋白用作氧传感器等。理解HBPs的结构对于阐明其分子机制以及进行理性设计以提升工业应用性能至关重要。然而，许多HBPs缺乏实验解析的结构，常依赖计算预测模型，但这类模型（尤其是血红素配位等细节）的可靠性尚未明确。近年来，基于人工智能（AI）的蛋白质结构预测工具如AlphaFold3（AF3）在预测蛋白质-配体相互作用方面展现出高精度，但其对HBPs的预测准确性尚未被系统评估。因此，研究人员开展本项研究，旨在评估AF3预测人类HBPs结构（包括血红素结合状态）的可靠性与局限性。

**研究内容与结论**

研究人员选取了四种高质量晶体结构的人类HBPs：亚硫酸盐氧化酶细胞色素b₅结构域（SO-b₅，PDB: 1MJ4）、NADH细胞色素b₅氧化还原酶细胞色素b₅结构域（Ncb5or-b₅，PDB: 3LF5）、细胞色素b₅类型B（CYB5B，PDB: 3NER）和神经红蛋白（NGB，PDB: 7VQG）。利用AF3预测这些蛋白的脱辅基（apo）和全蛋白（holo）状态模型，并与相应的实验结构进行全面比较。研究结论主要包括：（1）AF3预测的HBP整体结构具有可靠的置信度评分（pTM、ipTM、pLDDT），血红素分子能成功对接至结合口袋，但在血红素结合口袋几何形状、配位距离、血红素构象（如丙酸基取向）等细节上与实验结构存在差异；（2）对于NGB，AF3未能预测到实验结构中存在的Cys46-Cys55二硫键，导致局部折叠显著不同，并影响了血红素结合环境；（3）分子动力学模拟显示，AF3预测的NGB与实验NGB在二硫键区域和血红素结合区域的分子柔韧性存在显著差异。这些发现表明了AF3在预测HBPs整体结构上的潜力，同时也揭示了其在精细结构特征（特别是二硫键和血红素配位细节）方面的局限性。该论文发表在《International Journal of Molecular Sciences》。

**主要关键技术方法**

研究人员使用的方法主要包括：（1）蛋白质结构预测：将选定HBPs的氨基酸序列提交至AF3服务器，生成apo和holo状态的结构模型（默认参数下的单体预测）；（2）结构比较：使用PyMOL等软件对AF3预测模型与实验晶体结构（来源于蛋白质数据库PDB）进行叠加并计算均方根偏差（RMSD），分析血红素结合口袋体积（使用CB-Dock2计算）、血红素-蛋白质结合亲和力（使用PRODIGY-LIGAND网络服务器计算）以及铁-组氨酸配位距离；（3）分子动力学模拟：以AF3-NGB和实验NGB结构为初始模型，使用GROMACS（版本2025）结合CHARMM36力场，在310K（人体生理温度）下进行200 ns的全原子模拟，分析RMSD、回转半径（Rg）、溶剂可及表面积（SASA）和均方根波动（RMSF）等指标。未涉及试剂、培养、质粒构建或样本队列来源。

**研究结果介绍**

**2.1. AF3预测的HBPs**：通过对AF3生成的apo和holo状态模型进行置信度分析，所有HBP的pTM值均在0.81–0.91之间，ipTM值（除CYB5B的0.75–0.78外）在0.85–0.88之间，平均pLDDT值均在87.96–95.49之间，表明根据AF3置信度标准，大部分预测结构是可靠的，但CYB5B和Ncb5or-b₅分别显示相对较低的ipTM和pLDDT值。

**2.2. AF3预测与实验HBP结构比较**

* **2.2.1. 亚硫酸盐氧化酶细胞色素b₅结构域（SO-b₅）**：AF3-SO-b₅的holo状态与实验结构(EXP-SO-b₅)叠加RMSD为0.165–0.217 ?。血红素分子位置高度一致（差异<0.4 ?），但血红素周围残基（如Lys46）侧链构象以及血红素结合口袋体积（AF3: 648–747 ?³ vs 实验: 882 ?³)存在细微差异，结合亲和力相似（约11.2–11.3 kcal/mol）。

* **2.2.2. NADH细胞色素b₅氧化还原酶（Ncb5or-b₅）**：AF3预测的holo状态与实验结构叠加RMSD为0.131–0.286 ?。血红素铁原子位置最大偏差0.28 ?，血红素的丙酸基等基团构象存在差异。血红素结合口袋体积（AF3: 443–513 ?³ vs 实验: 439/475 ?³)相近，但结合亲和力（AF3: 10.34–10.73 kcal/mol vs 实验: 11.54–11.61 kcal/mol）存在差异。

* **2.2.3. 细胞色素b₅类型B（CYB5B）**：AF3预测模型与实验结构叠加RMSD为0.298–0.528 ?。血红素分子在结合口袋中的位置及乙烯基、甲基、丙酸基构象存在差异（部分归因于实验结构中血红素的内在柔性）。结合口袋体积（AF3: 401–453 ?³ vs 实验: 478–490 ?³）差异显著，但结合亲和力相近。

* **2.2.4. 神经红蛋白（NGB）**：AF3预测的holo状态与实验结构叠加RMSD为0.360–0.404 ?。关键差异在于AF3-NGB中未能形成Cys46-Cys55二硫键，导致Gln43-Ser59区域的折叠方式不同（AF3中该区域为环状，而实验结构中为α-螺旋）。这导致Tyr44和Asn45的侧链取向截然不同，进而改变了血红素丙酸基的构象和周围残基的位置。血红素结合口袋体积（AF3: 943–1071 ?³ vs 实验: 908 ?³）和结合亲和力（约11.3–11.4 kcal/mol）相似，但血红素环境细节差异显著。

**2.3. NGB的分子动力学模拟**：对AF3-NGB和EXP-NGB的apo和holo状态进行200 ns MD模拟。RMSD分析显示EXP-NGB波动更大（平均RMSD约2.04 ?和1.52 ?）高于AF3-NGB（约1.10 ?和1.13 ?）。RMSF分析表明EXP-NGB中Cys46-Cys55二硫键区域及其周围的血红素结合区域表现出更高的柔韧性，而AF3-NGB中对应区域由于缺乏二硫键反而相对刚性，且血红素结合使该区域柔性进一步降低。这表明AF3预测模型在下游计算（如MD）中会产生对分子柔韧性的误导性解释。

**讨论与结论总结**

讨论部分指出，AF3对SO-b₅、Ncb5or-b₅和CYB5B的预测总体上可靠，可用于识别血红素结合位点和关键相互作用残基，为蛋白质工程提供基础。但对于NGB，AF3未能预测Cys46-Cys55二硫键，导致结构和动力学特征显著偏离实验。研究人员通过测试溶菌酶和奇异果甜蛋白证实AF3能准确预测某些二硫键，表明NGB的失败是个例，凸显了AF3在预测特定二硫键上的局限性。此外，AF3对血红素分子的精细变形（如皱褶、翘曲等振动模式）的捕捉能力有限，且在预测时不能明确指定铁氧化态、电荷态或组氨酸质子化状态，因此解释血红素配位几何和功能状态时需谨慎。研究结论部分翻译如下：

**结论**：AF3-HBP模型可以提供HBPs可靠的总体结构，并成功将血红素分子对接至HBPs的血红素结合口袋中。然而，与血红素结合相关的精细结构特征（尤其在血红素结合空腔内）与实验结构中观察到的存在差异。特别地，AF3未能预测到NGB中观察到的二硫键形成。这些结果为AF3在HBPs结构分析中的应用提供了宝贵见解，并可能有助于未来AF3算法的改进。