四种新型环境细菌分离株对黄曲霉毒素的生物降解：比较基因组学与分子对接概念验证研究

《Foods》：Biodegradation of Aflatoxins by Four Novel Environmental Bacterial Isolates: A Comparative Genomic and Molecular Docking-Proof-of-Concept Study

　　

黄曲霉毒素污染仍然是食品和饲料安全的关键挑战，需要有效的缓解策略。本研究通过鉴定四种能够降解四种主要黄曲霉毒素同类物（B₁、B₂、G₁和G₂）的新型细菌候选菌株来应对这一挑战。定量高效液相色谱（HPLC）和定性液相色谱-高分辨质谱（LC-HRMS）结合串联质谱（MS/MS）分析显示，分离的细菌菌株实现了86–94%的高降解效率，并通过独特的脱甲基和脱羧途径进行。牛津纳米孔（Oxford Nanopore）全基因组测序将这些高效降解菌鉴定为大肠杆菌（Escherichia coli）EC₂和肺炎克雷伯菌（Klebsiella pneumoniae）KP₄、KP₅和KP₆。为阐明该表型的基因组基础，研究人员对两种非降解大肠杆菌菌株进行了比较基因组分析。功能注释显示，仅在降解菌株中存在完整的芳香族降解模块和完整的酶促级联反应，证明这种生物降解能力具有高度菌株特异性。此外，对活性级联反应中特有的关键推定酶进行分子对接，结果显示其与所有四种黄曲霉毒素同类物具有非常有利的结合亲和力（高达?8.70 kcal/mol），并通过稳健的氢键和疏水网络得以稳定。这些发现为黄曲霉毒素生物转化提供了新颖的基因组和结构蓝图，突出了用于未来研究的特定推定候选酶paaH和mhpB。

黄曲霉毒素（aflatoxin）是由曲霉属（Aspergillus）真菌产生的有毒次级代谢产物，其中黄曲霉毒素B₁（AFB₁）被国际癌症研究机构（IARC）列为1类致癌物，严重威胁全球食品安全和人类健康。传统的物理和化学脱毒方法存在成本高、营养损失及产生有毒副产物等缺陷，微生物降解法因其绿色可持续而受到广泛关注。然而，尽管已分离出多种具有降解能力的微生物，其分子机制仍不明确；传统筛选依赖生理生化应激实验，难以揭示完整代谢通路，且现有数据库缺乏基因组背景。因此，研究人员构建了一个多层级发现框架，旨在系统解析黄曲霉毒素生物降解的遗传基础与关键酶。

研究人员从埃及开罗、亚历山大和达米埃塔的受污染环境中采集了12份样本，包括牛粪便、土壤、霉变谷物饲料、生牛奶和Romy奶酪。通过以香豆素为唯一碳源的选择性富集培养，分离出4株香豆素耐受菌。体外降解实验结合活菌、无细胞上清液和热灭活对照，明确了降解过程主要由胞外酶催化，而非物理吸附。采用Oxford Nanopore长读长全基因组测序完成菌株鉴定和基因组组装，并通过比较基因组学、代谢产物鉴定和分子对接，最终揭示了候选酶及其与底物的结合模式。

本研究的主要技术方法包括：以香豆素作为黄曲霉毒素结构替代物进行选择性富集；采用活菌、无细胞上清液（CFS）和热灭活对照组对比，结合HPLC定量区分生物降解与生物吸附；利用Oxford Nanopore测序获得完整基因组，并通过KEGG模块重构和CARD数据库筛查进行功能注释与安全评估；采用UPLC-HR-MS/MS结合质量偏移分析鉴定降解产物；利用ProTox-3.0预测产物毒性；最后通过Quick Vina 2分子对接评估酶-配体结合亲和力。样本来源为埃及三地的12份农业和环境样本。

### 3.1 环境候选菌株的筛选和选择
通过以0.5%（w/v）香豆素作为唯一碳源的选择性富集培养，从12份样本中筛选出4株耐受菌（编号2、4、5、6），它们能在此培养基上快速生长，而阴性对照无生长，表明这些菌株具备裂解香豆素内酯环的潜力。

### 3.2 体外生物降解试验
#### 3.2.1 香豆素耐受性筛选谱
4株菌中，菌株6耐受率最高（96%），菌株4为66.90%，菌株5为59.50%，菌株2最低（51.30%）。菌株1和3因不能在香豆素培养基上生长，被用作非耐受对照。

#### 3.2.2 活菌、无细胞上清液和热灭活生物质的比较评估
活菌处理组显示了最高的总黄曲霉毒素去除率，菌株5达到94.02%±0.35%；无细胞上清液去除率为71.04%±2.11%至81.31%±1.52%；而热灭活对照组仅去除9.14%±3.68%至16.62%±1.72%。双向方差分析显示，处理类型是主要变异来源（偏η2=0.952，p<0.001），而菌株差异不显著（p=0.149）。这证明去除过程以酶促降解为主，物理吸附贡献极小。

#### 3.2.3 黄曲霉毒素去除的时间过程分析
以菌株5为代表性菌株，在72小时内监测四种毒素的去除动态。结果表明，去除从对数生长早期（6-12小时）开始，24小时时AFB₁和AFG₁去除率分别达99.43%±0.15%和97.88%±2.32%，而AFG₂和AFB₂分别为84.67%±1.20%和74.33%±1.45%。72小时时，AFB₁和AFG₁接近完全去除（99.83%±1.04%和99.70%±1.06%），AFG₂和AFB₂达89.25%±0.45%和92.58%±0.62%。差异与呋喃环的饱和度有关。

### 3.3 菌株的鉴定和表征
根据Bergey’s Manual进行生化鉴定，结合基因组测序，将菌株2鉴定为大肠杆菌EC₂，菌株4、5、6鉴定为肺炎克雷伯菌KP₄、KP₅、KP₆。

### 3.4 比较基因组组装、注释和生物安全分析
基因组组装和KEGG注释显示，降解菌株EC₂和KP₄-KP₆拥有完整的苯甲酸降解模块（M00551）、儿茶酚邻位裂解模块（M00568）和苯乙酸降解模块（M00878），而非降解菌株EC₁和EC₃缺少或仅有部分模块。CARD数据库筛查发现，这些降解菌株携带多种内在多药外排系统和抗生素耐药基因（如SHV-27），提示它们属于条件致病菌，不宜直接应用于食品环境。

### 3.5 推定生物转化代谢物的质谱分析
通过UPLC-HR-MS/MS分析降解菌株CFS中AFB₁的转化产物，检测到脱甲基产物P1（m/z 299.0558，[C₁₆H₁₁O₆]⁺）和脱羧产物P2-P4（m/z 269.0809，[C₁₆H₁₃O₄]⁺），表明降解途径包括脱甲基和脱羧反应。

### 3.6 推定降解产物的毒理学和细胞毒性预测
ProTox-3.0预测结果表明，亲本AFB₁具有高致突变性（0.95）和致癌性（0.68），而脱甲基产物P1致突变概率降至0.56，脱羧终产物P2预测无致突变性、无线粒体毒性且不被细胞色素P450激活，说明生物转化显著降低了毒性。

### 3.7 分子对接
对候选酶与四种黄曲霉毒素进行分子对接，苯乙酸途径中的3-羟基丁酰辅酶A脱氢酶（paaH）表现出最高亲和力，结合能范围为?8.00至?8.70 kcal/mol，其中对AFG₂亲和力最高（?8.70 kcal/mol）；2,3-二羟基苯丙酸1,2-双加氧酶（mhpB）结合能为?6.0至?6.7 kcal/mol。来自非降解菌株的mhpE作为阴性对照，显示不利的正或负空间契合能。相互作用分析表明，氢键和疏水接触是稳定复合物的主要作用力。

讨论部分指出，香豆素筛选策略有效富集了具有内酯环裂解能力的菌株，活菌与CFS的高效降解和热灭活对照的低活性证实了酶促机制。比较基因组分析揭示，降解菌株独有的完整芳香族降解模块是其功能的关键，且该能力具有高度菌株特异性。尽管这些菌株携带AMR基因，不适合直接应用，但它们为酶挖掘提供了宝贵资源。研究结论强调，通过整合全基因组测序、代谢物谱和分子对接，本研究为黄曲霉毒素生物降解提供了基因组和结构蓝图，候选酶paaH和mhpB经异源表达后有望成为安全高效的脱毒工具。

研究结论部分翻译如下：总之，本研究通过将经典生理筛选与高通量基因组分辨率相结合，解决了黄曲霉毒素缓解中的关键瓶颈。通过结合Oxford Nanopore长读长测序、UPLC-HR-MS/MS质量偏移分析和分子对接，研究人员绕过了片段化标记基因数据库和未解析生化途径的限制，提供了生物学资源和预测性结构见解。本研究为公共数据库贡献了六个完整细菌基因组，建立了清晰的比较基因组基础，区分了表型活性黄曲霉毒素降解菌与非降解菌。活性降解菌表现出卓越的生物缓解能力，对黄曲霉毒素B₁、B₂、G₁和G₂的去除率达86–94%。计算机毒性模型表明，预测的脱甲基和脱羧代谢物比亲本毒素具有更低的致突变潜力。这些菌株中存在可转移的抗菌素耐药性（AMR）基因，转变了应用范式，使这些菌株成为有价值的非活体、酶富集资源库，用于定向生物挖掘。通过挖掘KEGG通路谱，研究人员鉴定出仅存在于活性降解菌株中的候选氧化还原酶、水解酶和双加氧酶。分子对接模拟解析了酶-配体相容性结构，产生了异常高的对接分数，表明与黄曲霉毒素具有强结合亲和力。因此，本研究为将这些预测性、结构相容的候选酶转化为异源表达、安全且高度纯化的生物催化剂奠定了基础。未来工作将集中在这些优先酶的定量表达、纯化和湿实验验证上，以释放其作为安全高效工业黄曲霉毒素去污生物制剂的全部潜力。