铜绿假单胞菌弹性蛋白酶（LasB）与合成聚合物结构片段的分子模拟相互作用分析

《Applied Microbiology》：In Silico Interaction Profiling of Pseudomonas aeruginosa Elastase (LasB) with Structural Fragments of Synthetic Polymers Afrah I. Waheeb, Saleem Obaid Gatia Almawla, Mayada Abdullah Shehan, Sameer Ahmed Awad, Mohammed Mukhles Ahmed and Saja Saddallah Abduljaleel

　　

背景：合成塑料在环境中的持久性积累及微塑料的广泛分布，迫切需要探索能够与聚合物材料发生相互作用甚至实现降解的生物机制，具有催化灵活性的微生物酶是该类研究的重点候选对象。目的：本研究旨在解析铜绿假单胞菌弹性蛋白酶（LasB）与典型合成塑料代表性结构片段的分子相互作用模式及动态稳定性，评估该酶对聚合物衍生底物的适配性。方法：研究人员从蛋白质数据库（Protein Data Bank, PDB）中获取LasB晶体结构（PDB ID: 1EZM），利用AutoDock4.2.6工具完成蛋白预处理；从PubChem数据库中检索聚酰胺（PA）、聚氯乙烯（PVC）、聚碳酸酯（PC）、聚对苯二甲酸乙二醇酯（PET）、聚甲基丙烯酸甲酯（PMMA）和聚氨酯（PUR）对应的聚合物衍生配体，采用MMFF94力场完成几何优化；以催化口袋（含Zn2+离子）为中心设置网格盒，利用AutoDock Vina开展分子对接模拟，结合能计算单位为kcal/mol；采用PyMOL和BIOVIA Discovery Studio分析结合构象、相互作用残基及分子间作用类型，以金属蛋白酶已知抑制剂磷酰胺素（Phosphoramidon）作为阳性对照验证对接方案可靠性；进一步采用Desmond引擎、OPLS4力场及显式溶剂模型开展50 ns分子动力学（MD）模拟，通过均方根偏差（RMSD）、均方根波动（RMSF）、回转半径（Rg）、氢键分析及溶剂可及表面积等参数评估复合物结构稳定性。结果：对接实验显示不同测试聚合物的结合亲和力存在差异，聚碳酸酯（?5.774 kcal/mol）和聚氨酯（?5.707 kcal/mol）与LasB催化口袋的相互作用最强，聚酰胺（?5.277 kcal/mol）和PET（?4.483 kcal/mol）次之，PMMA和PVC亲和力较弱；关键参与相互作用的残基包括Glu141、His140、Val137、Arg198、Tyr114和Trp115，主要通过疏水相互作用、π-阳离子相互作用和范德华力维持复合物稳定。分子动力学模拟显示所有复合物均达到稳定状态，RMSD值处于1.0–3.2 ?的可接受范围，与稳定蛋白-配体体系的特征一致；阳性对照磷酰胺素的RMSD为1.205 ?，处于该稳定区间内。主成分分析（PCA）显示LasB的不同配体结合构象存在分化：PA对应紧凑构象，PC和PVC对应中间态，PUR、PMMA和PET对应伸展构象，反映LasB结合口袋的结构稳定性和适配性。结论：研究结果表明LasB具有结构灵活的催化口袋，可容纳多种聚合物衍生配体，为分子层面解析酶与聚合物的识别机制提供了参考，同时提示LasB可能参与环境中微生物与合成塑料的相互作用过程。

铜绿假单胞菌弹性蛋白酶与合成聚合物结构片段的分子模拟研究解读

研究背景与意义

合成塑料的大规模生产与使用导致其在土壤和水生生态系统中持续累积，传统不可降解塑料不仅降低土壤肥力，还与内分泌干扰、神经毒性及肝肾损伤等严重健康风险相关。当前全球塑料年产量预计2060年将超过12亿吨，其高分子量、强疏水性及缺乏天然降解途径的特性，使得常规废物管理策略难以应对。在此背景下，基于微生物酶的生物催化降解成为塑料污染治理的重要研究方向。铜绿假单胞菌弹性蛋白酶（LasB）属于热解素样M4家族锌依赖性金属内肽酶，具有广泛的底物特异性和催化可塑性，能够水解弹性蛋白、胶原蛋白等多种宿主蛋白，但其与合成聚合物键的相互作用尚未得到明确验证。本研究通过分子模拟手段，首次系统评估LasB对工业常见聚合物结构片段的结合能力，为理解微生物-塑料界面相互作用提供分子层面的证据。

关键技术方法

研究人员采用多尺度计算生物学策略开展研究：首先从PDB数据库获取分辨率1.50 ?的LasB晶体结构（PDB ID: 1EZM），保留催化必需的Zn²⁺离子并完成预处理；选取6种典型工业塑料（PA、PVC、PC、PET、PMMA、PUR）的低分子量结构类似物作为配体，以磷酰胺素作为阳性对照；采用AutoDock Vina开展分子对接，设置包含催化口袋的20×20×20 ?网格盒，通过结合能排序筛选最优结合模式；利用PyMOL和BIOVIA Discovery Studio解析相互作用类型；进一步采用Desmond软件开展50 ns全原子分子动力学模拟，基于OPLS4力场构建TIP3P显式溶剂环境，通过RMSD、RMSF、Rg、氢键数量及溶剂可及表面积等多参数评估复合物动态稳定性；最后通过主成分分析（PCA）解析不同配体诱导的蛋白构象变化特征。

研究结果

分子对接结果

聚碳酸酯（PC）表现出最高的预测结合亲和力（?5.774 kcal/mol），其次是聚氨酯（PUR，?5.707 kcal/mol），聚酰胺（PA，?5.277 kcal/mol）和PET（?4.483 kcal/mol）处于中等水平，PMMA（?3.051 kcal/mol）和PVC（?2.520 kcal/mol）亲和力最弱。PC可与LasB形成π-π T型堆积、π-阳离子及多重疏水相互作用，PUR则通过π-σ相互作用、π-堆积及碳氢键维持稳定结合；PA和PET分别依赖烷基/π相互作用及π-烷基接触实现结合；PVC仅存在弱疏水接触，无稳定氢键形成。所有配体的关键相互作用残基均集中在Glu141、His140、Val137、Arg198、Tyr114和Trp115，验证了LasB催化口袋的多功能性。

分子动力学模拟结果

所有LasB-配体复合物在50 ns模拟中均保持结构完整性：PC和PUR复合物的RMSD稳定在2.7–3.2 ?，PA和PMMA复合物RMSD低于1.8 ?，PVC复合物RMSD为1.5–2.2 ?，均未达到显著构象漂移阈值；回转半径（Rg）波动幅度小于0.5 ?，表明聚合物结合未引起蛋白整体折叠改变。PC和PET复合物中观察到持续的π-π堆积相互作用，为主要稳定力；PA和PVC复合物以瞬时氢键和范德华接触为主。径向分布函数（RDF）分析显示所有配体均与多个周围残基形成短程接触，而非局限于单一结合位点，证实LasB通过表面多位点吸附模式与聚合物片段相互作用。阳性对照磷酰胺素复合物的RMSD为1.205 ?，氢键占有率超过70%，验证了模拟体系的可靠性。

主成分分析结果

PCA将LasB-配体复合物分为三个构象簇：PUR、PMMA和PET对应伸展构象（PC1和PC2正值区域），PC和PVC对应中间态构象，PA对应紧凑构象（双轴负值区域）。不同构象簇反映了配体化学结构对结合口袋的动态调控作用，与对接和动力学模拟的稳定性结果一致。

讨论与结论

讨论部分明确指出，本研究仅证明LasB对聚合物片段的结构兼容性，而非证实其具有塑料降解功能——LasB作为蛋白酶，其催化机制演化以适应肽键水解，无法直接断裂聚合物主链的化学键。但研究结果仍具有重要环境意义：在自然生物膜界面，非经典塑料降解酶可能通过与风化塑料表面的氧化片段、添加剂或低分子量降解中间产物相互作用，参与微生物对塑料的定殖过程。研究同时承认局限性：低分子量配体无法完全模拟真实塑料的结晶度、表面异质性等宏观特性；50 ns模拟未能覆盖长时程构象转变；未考虑生物膜环境的复杂理化条件。

结论部分强调，LasB的催化口袋具有显著的结构灵活性，可通过疏水作用、π相互作用及范德华力稳定结合多种聚合物衍生片段。该发现为理解酶-聚合物分子识别提供了基础数据，可作为后续实验研究的参考依据。研究人员建议未来结合MM/PBSA自由能计算、体外酶活测定及聚合物降解表征，进一步明确该类相互作用的环境功能。本研究发表于《Applied Microbiology》，为塑料生物降解领域的计算生物学研究提供了新的视角。