《Bacteria》：Integrated Structure-Based Virtual Screening and Molecular Dynamics Simulations Identify Marine Natural Products as Potential Inhibitors of AmvR from Acinetobacter baumannii
编辑推荐：
抗微生物耐药性（AMR）已成为全球公共卫生关注的新兴问题，亟需鉴定新型药物靶点与生物活性分子以对抗耐药细菌病原体。本研究将AmvR（TetR家族转录调控因子）作为推定的潜在药物靶点加以研究，因其在细菌适应性及相应抗微生物耐药机制中发挥关键作用。研究人员采用结构
抗微生物耐药性（AMR）已成为全球公共卫生关注的新兴问题，亟需鉴定新型药物靶点与生物活性分子以对抗耐药细菌病原体。本研究将AmvR（TetR家族转录调控因子）作为推定的潜在药物靶点加以研究，因其在细菌适应性及相应抗微生物耐药机制中发挥关键作用。研究人员采用结构基础虚拟筛选方法（SBVS），从综合海洋天然产物数据库（CMNPD）和海藻代谢物数据库（SWMD）中鉴定海洋来源的新型抑制剂。首先以共结晶配体亚精胺（spermidine）验证对接方案（RMSD = 0.12 ?），随后针对AmvR结合位点筛选约32,200种海洋天然产物（MNPs）。基于对接评分（范围为?11.4至?20.3 kcal/mol）及对接阶段MMGBSA结合能估算（范围为?89.7至?104.4 kcal/mol），与对照配体相比，四种MNPs（CMNPD9680、CMNPD23795、CMNPD5249和SWMD-BS008）被鉴定为排名最高的候选物。药代动力学与毒性分析表明所选MNPs具有有利的类药性和安全性特征。通过500 ns分子动力学（MD）模拟进一步评估apo AmvR及配体结合复合物的结构稳定性。对RMSD、RMSF、Rg、SASA、氢键动力学、动力学后MM/PBSA结合自由能、主成分分析（PCA）、自由能景观（FEL）映射及基于DSSP的二级结构评估的分析一致表明，蛋白–配体相互作用在整个模拟期间保持稳定。在所鉴定的MNPs中，CMNPD9680表现出最高的结合亲和力、最有利的对接阶段MM/GBSA结合自由能以及在AmvR结合口袋内稳定的相互作用网络。基于Morgan指纹的结构相似性分析揭示，入围的海洋化合物与针对鲍曼不动杆菌（Acinetobacter baumannii）的实验验证抗菌化合物具有共同的保守特征，使其成为有前景的抗菌先导物。总体而言，这些发现证明海洋天然产物作为有前景的AmvR抑制剂的潜力，并鉴定CMNPD9680作为针对鲍曼不动杆菌进一步实验研究的有前景的计算候选物。
**研究背景与问题**
抗微生物耐药性（AMR）的出现在全球范围内构成重大公共卫生威胁，对感染性疾病的预防和管理造成巨大挑战。世界卫生组织（WHO）数据显示，细菌AMR导致全球数百万人死亡，并给医疗系统带来经济负担。多重耐药（MDR）细菌感染的增加降低了标准治疗药物的效率，表明需要开发新的治疗策略。在各类病原微生物中，不动杆菌属（Acinetobacter）因其耐受多样环境条件并快速获得抗生素耐药性的能力而备受关注，可导致肺炎、血流感染、伤口感染和尿路感染等多种感染类型，尤其在免疫功能低下患者中。不动杆菌属的卓越适应性主要归因于其复杂的调控网络，协调应激反应、毒力、代谢适应和抗微生物耐药性。AmvR属于TetR家族转录抑制因子，该家族是最大且研究最广泛的细菌转录调控因子家族之一，调控与多重耐药、膜转运、应激反应、生物膜形成、代谢和毒力相关的基因表达。近期研究揭示，破坏TetR介导的信号传导可能导致细菌适应性和耐药机制显著降低。因此，TetR家族调控因子已成为开发下一代抗菌药物的有吸引力的治疗靶点。与聚焦直接抑制细菌生长的传统抗生素不同，抑制转录调控因子的策略以更间接的方式抑制细菌病原体，不施加选择压力促使抗生素耐药性的发展。此外，AmvR晶体结构的可用性为结构基础药物发现提供了优良的结构框架。天然产物长期以来是治疗剂的重要来源，海洋生物尤其栖息于独特生态环境中，促进结构多样、具有广泛生物活性的次级代谢产物的生物合成。综合海洋天然产物数据库（CMNPD）和海藻代谢物数据库（SWMD）作为化学多样性海洋来源化合物的宝贵储库，适用于药物发现应用。计算方法通过整合结构基础虚拟筛选、分子对接、药代动力学分析和分子动力学模拟等技术，在药物发现中变得越来越重要，允许快速筛选大量化合物而不产生高昂的实验成本和时间。
**研究内容与结论**
研究人员采用集成计算方法，通过结构基础虚拟筛选、类药性与毒性评估、分子动力学模拟、主成分分析、自由能景观映射和二级结构分析，从CMNPD和SWMD数据库中鉴定AmvR的海洋来源抑制剂。研究得出以下结论：从约31,000种海洋天然产物中鉴定出四种有前景的候选物（CMNPD9680、CMNPD23795、CMNPD5249和SWMD-BS008），均表现出优于对照配体亚精胺的结合亲和力；所有化合物满足Lipinski五规则且具有有利的ADMET特征；500 ns MD模拟表明AmvR-MNP复合物在整个模拟期间保持结构完整性；PCA和FEL分析证实形成了能量上有利的构象状态；CMNPD9680表现出最强的结合亲和力和最有利的MM/GBSA结合自由能；基于Morgan指纹的结构相似性分析表明入围的海洋化合物与实验验证的抗鲍曼不动杆菌抗菌化合物共享保守结构特征。这些计算发现证明海洋天然产物是新型AmvR抑制剂的宝贵资源，为针对鲍曼不动杆菌的未来实验验证和结构指导的发现提供了计算基础。
**关键技术方法**
研究人员采用的主要技术方法包括：结构基础虚拟筛选（SBVS），使用Glide模块的虚拟筛选工作流（VSW），采用高通量虚拟筛选（HTVS）、标准精度（SP）和额外精度（XP）分级对接协议；对接阶段MMGBSA结合能估算；使用QikProp模块进行药代动力学预测和ProTox-3.0服务器进行毒性预测；使用GROMACS v2022进行500 ns分子动力学模拟（蛋白拓扑采用CHARMM27力场，配体拓扑采用SwissParam服务器生成）；动力学后MM/PBSA结合自由能计算（使用g_mmpbsa工具，取最后100 ns轨迹的1000帧）；主成分分析（PCA）基于Cα原子位置波动的协方差矩阵；自由能景观（FEL）映射通过将MD轨迹投影到PC1-PC2平面；基于DSSP算法的二级结构分析；以及基于Morgan圆形指纹（1024位）的Tanimoto相似性分析和K-means无监督聚类（使用R程序version 4.4.1的rcdk包）。参考数据集为ChEMBL数据库中350个实验验证的抗鲍曼不动杆菌化合物。样本队列来源为CMNPD数据库（约31,000种独特海洋来源化合物）和SWMD数据库（约1190种化合物），靶蛋白晶体结构来源于RCSB蛋白质数据库（PDB ID: 8WP7）。
**研究结果**
3.1 使用对照分子亚精胺验证对接方案：通过重新对接共结晶配体亚精胺验证Glide模块的对接方案，获得的RMSD值为0.12 ?，远低于公认的验证阈值2.0 ?，证实对接方案准确重现了亚精胺的实验观察结合模式，证明了整个对接程序的可靠性。
3.2 通过虚拟筛选鉴定AmvR的潜在海洋来源抑制剂：使用Glide虚拟筛选工作流筛选CMNPD和SWMD数据库中的MNPs，鉴定出四种海洋来源命中候选物，均表现出高于共结晶对照配体亚精胺的结合亲和力。CMNPD9680为最有前景的命中化合物，Glide对接评分为?20.3 kcal/mol，MMGBSA结合自由能为?104.4 kcal/mol，与Leu87、Glu106和Glu164形成氢键，与Trp102和Trp126形成π-π堆积相互作用，并与Glu106和Glu164形成盐桥。分析表明Trp102、Trp126、Glu106、Glu164、Leu87和Ile85残基在配体识别和结合中发挥重要作用。
3.3 所选海洋天然产物的药代动力学与毒性评估：所有所选化合物分子量在305.50至336.40 Da之间，满足Lipinski五规则无违反。CMNPD23795和CMNPD5249的Caco-2细胞通透性预测值分别高达748.9和1160.1 nm/s。毒性预测显示CMNPD5249具有最有利的安全性特征（LD50为3500 mg/kg，毒性等级V），CMNPD23795次之（LD50为1600 mg/kg，毒性等级IV）。CMNPD9680和SWMD-BS008的LD50分别为652 mg/kg和920 mg/kg，具有中等毒性但在可接受水平内。毒性热图显示CMNPD23795对免疫毒性终点有活性（预测评分0.95），其余化合物所有毒性终点均显示无活性。
3.4 Apo和配体结合AmvR复合物的分子动力学模拟分析：500 ns MD模拟显示所有系统RMSD达到平衡并保持在可接受范围内（<0.30 nm）。AmvR-CMNPD23795和AmvR-SWMD-BS008复合物RMSD波动最小。配体RMSD显示CMNPD23795具有最低值。RMSF分析显示配体结合未显著改变AmvR的内在动态行为，大多数残基波动小于0.15 nm，170–180区域对应柔性环区出现显著峰值。Rg分析显示所有复合物维持约1.85至1.95 nm的稳定值。SASA分析显示所有复合物具有稳定曲线。氢键动力学显示CMNPD9680在整个模拟期间维持持久氢键相互作用。
3.5 动力学后结构稳定性分析：通过叠加100、200、300、400和500 ns的代表性结构，发现所有模拟系统维持二级和三级结构排列。对照复合物叠加RMSD值为1.45 ?，AmvR-CMNPD23795复合物显示最低叠加RMSD值（1.09 ?），其次为AmvR-CMNPD5249（1.13 ?）和AmvR-SWMD-BS008（1.38 ?）。AmvR-CMNPD9680显示相对较高的叠加RMSD值（2.42 ?），但整体蛋白结构保持良好，无结构失稳或去折叠证据。
3.6 主成分分析（PCA）与自由能景观（FEL）分析：PCA显示apo AmvR探索相对广泛的构象空间，配体结合约束了过度原子运动导致更紧凑的构象采样。CMNPD9680、CMNPD5249和SWMD-BS008显示高度紧凑的构象聚类。FEL分析显示apo蛋白具有较宽的低能盆地，所有配体结合复合物显示更窄更深的能量极小值。AmvR-CMNPD9680和AmvR-SWMD-BS008具有高度集中的低能盆地和明确的全局极小值。配体结合系统未出现多个孤立能量盆地，表明500 ns模拟期间未发生显著构象转变。
3.7 基于MM/PBSA的结合自由能分析：所有四种所选分子结合复合物显示比对照复合物（?30.59 ± 65.31 kJ/mol）更负的平均MM-PBSA结合能。CMNPD9680显示最有利的结合自由能（?216.977 ± 20.77 kJ/mol），其次为CMNPD23795（?150.064 ± 20.72）、CMNPD5249（?138.506 ± 48.60）和SWMD-BS008（?134.482 ± 13.87）。范德华贡献对所有四种MNP复合物均比对照复合物更有利，表明疏水和紧密接触相互作用对MNP结合起重要贡献。
3.8 二级结构稳定性分析：DSSP分析显示AmvR的整体二级结构组成在所有模拟系统中保持良好保守，α-螺旋构成最丰富的结构组分。在转角和无规卷曲区域观察到微小波动，主要位于柔性环段。配体结合后未观察到明显的螺旋到卷曲转变或其他重大二级结构变化，表明所选海洋天然产物不影响蛋白质的结构完整性。
3.9 所选海洋化合物与实验验证抗鲍曼不动杆菌化合物的结构相似性分析：Tanimoto相似性分析显示所有四种海洋化合物与实验验证抗菌化合物具有结构相似性，最大Tanimoto相似性评分分别为0.975、0.982、0.976和0.986。每种海洋化合物与至少一种实验验证的ChEMBL化合物共享共同化学骨架或关键结构片段。
3.10 基于Morgan指纹的结构聚类：K-means聚类分析显示入围的海洋化合物与结构相关的ChEMBL化合物聚类在一起，表明海洋化合物具有与实验验证抗鲍曼不动杆菌化合物相似的结构。气泡图直观展示了每种海洋化合物与其最相似ChEMBL类似物之间的化学关系。
**讨论与结论**
多重耐药鲍曼不动杆菌加剧了对参与调控细菌存活和适应性耐药机制的替代潜在治疗靶点的探索。TetR家族转录调控因子成员因其调控抗生素耐药、毒力、代谢和环境适应相关基因表达而备受关注。AmvR作为这些调控因子之一，显示出作为治疗靶点的巨大潜力。阻断AmvR功能可能干扰细菌调控网络，使鲍曼不动杆菌对不利条件的适应性降低，从而对抗生素更敏感。海洋生态系统是化学多样性天然产物的显著来源，具有广泛抗菌特性。本研究调查AmvR作为海洋天然产物基础发现的潜在靶点，靶向AmvR的理由不仅在于抑制细菌生长的效果，还在于其破坏靶蛋白调控功能的能力，可能导致细菌适应性和抗生素治疗适应性降低，并可能抑制耐药相关表型的发展。约31,000种海洋天然产物的结构基础虚拟筛选产生了四种有前景的候选物，展现不同的结构架构。对接方案通过重新对接共结晶亚精胺成功验证（RMSD = 0.12 ?）。相互作用分析表明Trp102、Trp126、Glu106和Glu164等残基通过氢键、π-堆积和静电相互作用在配体识别中发挥重要作用。所有所选海洋来源化合物满足Lipinski五规则且显示有利的ADMET特征。MD模拟过程中，CMNPD23795显示最小RMSD值（1.09 ?），CMNPD9680显示相对较高的RMSD值（2.42 ?），但未出现结构崩溃或去折叠。动力学后MM/PBSA评估证实了CMNPD9680的高结合亲和力，与对接结果一致。PCA和FEL映射提供了支持蛋白–配体复合物稳定性的额外证据。DSSP分析揭示AmvR二级结构在模拟过程中保持主要保守。基于化学相似性的分析显示发现的海洋命中化合物与抗鲍曼不动杆菌分子具有高结构相似性。研究局限性包括：每个系统使用单一轨迹分析，需要独立重复模拟以提供更可靠的构象采样和统计验证；缺乏实验验证，所鉴定化合物被视为计算命中候选物而非验证的治疗先导物；CMNPD9680显示预测QPlogHERG值为?6.7，表明潜在的hERG相关风险值得进一步研究。未来研究将聚焦于通过AmvR结合、功能抑制和抗鲍曼不动杆菌抗菌敏感性测定进行化合物验证。
研究结论部分翻译：本研究采用由结构基础虚拟筛选、分子对接、MM/GBSA结合自由能计算、药代动力学与毒性评估、分子动力学模拟、PCA、FEL映射和基于DSSP的二级结构分析组成的集成计算流程，从海洋来源天然产物中鉴定AmvR调控蛋白的潜在抑制剂。多级虚拟筛选鉴定出CMNPD9680、CMNPD23795、CMNPD5249和SWMD-BS008为有前景的候选物，表现出有利的结合亲和力、稳定的蛋白–配体相互作用和理想的类药性特征。分子动力学模拟研究表明AmvR-MNP复合物在采样的500 ns轨迹内维持整体结构完整性，PCA和FEL分析证实了能量上有利的构象状态的形成。在所鉴定的化合物中，CMNPD9680表现出最强的结合亲和力和最有利的MM/GBSA结合自由能，强调其作为靶向AmvR的有前景的计算命中候选物的潜力。使用Morgan指纹的结构相似性分析表明入围的海洋化合物与实验验证的抗鲍曼不动杆菌抗菌化合物共享保守结构特征。总体而言，本研究的计算发现证明海洋天然产物代表新型AmvR抑制剂的宝贵资源，为针对鲍曼不动杆菌的未来实验验证和结构指导的发现提供了计算基础。
