吲哚基甲基硫代葡萄糖苷与吲哚-3-甲醇抑制艰难梭菌二元毒素的机制解析

《Scientifica》：Mechanistic Insights Into the Inhibition of Clostridioides difficile Binary Toxin by Indolylmethyl Glucosinolate and Indole-3-Carbinol

　　

艰难梭菌感染（Clostridioides difficile infection，CDI）仍然是重大医疗卫生挑战，其主要原因在于毒素介导的致病机制可导致严重胃肠道并发症。本研究采用计算方法与实验方法相结合的策略，考察来源于甘蓝（Brassica oleracea L.）的生物活性化合物对二元毒素CDTa亚基的抑制潜力。分子对接分析鉴定出吲哚基甲基硫代葡萄糖苷和吲哚-3-甲醇为最优候选分子，其与CDTa蛋白分别表现出较高结合亲和力，结合能为?9.1?kcal/mol和?8.75?kcal/mol。随后通过100?ns分子动力学（molecular dynamics，MD）模拟进一步评估配体-蛋白复合物的稳定性及动态行为，结果证实二者均可形成稳定相互作用。基于MM-GBSA（Molecular Mechanics Generalized Born Surface Area）计算的热力学评估显示其结合自由能具有有利特征，支持其作为有效抑制剂的潜力。与此同时，体外抑菌实验验证了这些化合物的抗菌效力：在100?μg/mL浓度下，吲哚基甲基硫代葡萄糖苷和吲哚-3-甲醇的最大抑菌圈分别为23.33?±?0.67?mm和22.67?±?0.33?mm。在最低抑菌浓度（minimum inhibitory concentration，MIC）和最低杀菌浓度（minimum bactericidal concentration，MBC）实验中，两种化合物均显示出显著抗菌活性，其中吲哚基甲基硫代葡萄糖苷的效力略高于吲哚-3-甲醇，前者MIC为10.33?±?0.72?μg/mL、MBC为23.33?±?1.36?μg/mL，后者MIC为11.33?±?0.27?μg/mL、MBC为26.67?±?3.60?μg/mL。此外，主成分分析（principal component analysis，PCA）和动态交叉相关矩阵（dynamic cross-correlation matrix，DCCM）结果表明，配体结合后CDTa呈现出不同的构象稳定化模式及相关运动特征，这支持其可能存在抗毒力相互作用。进一步的ADMET（absorption, distribution, metabolism, excretion, and toxicity）分析揭示了两种化合物在药代动力学与毒性特征方面的差异，为其治疗潜力提供了更多认识。总体而言，这些结果表明，来源于B. oleracea的化合物不仅对艰难梭菌（C. difficile）具有直接抗菌活性，而且在计算分析支持下可抑制CDTa，提示其还可能通过补充性的抗毒力机制促进其在CDI治疗中的应用潜力。

本文发表于《Scientifica》，围绕艰难梭菌（Clostridioides difficile）感染（CDI）的毒力干预与天然产物药物发现展开。艰难梭菌是临床重要的革兰阳性厌氧致病菌，常与医疗相关性腹泻、伪膜性结肠炎及重症胃肠道并发症相关。其致病过程高度依赖毒素，尤其是二元毒素CDT，其中CDTa为酶活性亚基，可对宿主细胞单体G-肌动蛋白Arg177进行ADP-核糖基化修饰，从而抑制肌动蛋白聚合，破坏F-肌动蛋白骨架，导致上皮屏障受损、细胞圆缩及炎症加重。尽管现有治疗方案如万古霉素、非达霉素和甲硝唑可用于CDI，但仍存在复发率高、菌群扰动明显、成本高或疗效下降等问题。因此，开发既具抗菌作用又可能削弱毒力的新型天然化合物，成为该领域的重要研究方向。

在这一背景下，研究人员选取甘蓝（Brassica oleracea）来源植物化学成分作为候选分子，聚焦艰难梭菌二元毒素的CDTa亚基，尝试从分子层面筛选潜在抑制剂，并结合体外抑菌实验验证其生物学活性。研究结果显示，吲哚基甲基硫代葡萄糖苷和吲哚-3-甲醇是两种最具代表性的候选分子。前者在计算模拟中表现出更稳定、更有利的结合特征，并且在体外实验中展现出略优的抑菌效果。研究整体提示，B. oleracea来源化合物可能兼具直接抗菌与补充性抗毒力双重作用，为CDI干预提供了新的候选方向。

本研究主要采用以下关键技术方法：首先从PubChem数据库挖掘51种B. oleracea来源化合物，并对其进行能量优化；随后以PDB数据库中CDTa结构（PDB ID: 6X41）为靶标开展盲法分子对接筛选。针对优势配体，进一步进行100?ns分子动力学（MD）模拟，结合RMSD、RMSF、SASA、回转半径（Rg）、氢键、MM/PBSA、主成分分析（PCA）和动态交叉相关矩阵（DCCM）评估复合物稳定性与构象变化。同时，采用pkCSM与SwissADME进行ADMET与类药性预测。实验部分使用来源于孟加拉国University of Rajshahi实验室、初始分离自CDI患者的艰难梭菌菌株，开展纸片扩散、MIC及MBC测定。

3.1. Molecular Docking Study
研究人员首先通过分子对接评估51种甘蓝来源化合物与CDTa蛋白的结合能力。结果显示，吲哚基甲基硫代葡萄糖苷和吲哚-3-甲醇的结合能分别为?9.1?kcal/mol和?8.75?kcal/mol，均优于对照药环丙沙星的?8.2?kcal/mol。吲哚基甲基硫代葡萄糖苷与ASN262、ARG302、ASN342、ASN263、PHE343等残基形成多重氢键，并伴随C-H键、π-π堆积及电荷吸引相互作用，提示其具有较强结合稳定性。吲哚-3-甲醇则与ALA31、GLU35、GLU47和THR177形成氢键，并在LYS43处产生π-烷基作用，显示较好的结合潜力。该部分结果表明，甘蓝来源化合物能够与CDTa潜在结合位点有效结合，具有作为毒力抑制剂的可能性。

3.2. MD Simulation
为验证分子对接结果的可靠性，研究人员对优势复合物进行了100?ns MD模拟，综合评估其动态稳定性。整体结果显示，对接获得的关键配体在动力学过程中总体保留了与CDTa的相互作用模式，说明初步筛选结果具有一定稳定性和生物物理合理性。

3.3. Analysis of RMSD
RMSD（均方根偏差）分析显示，吲哚基甲基硫代葡萄糖苷-CDTa复合物平均RMSD为1.922??，与apo蛋白的1.803??及环丙沙星对照的1.729??较为接近，说明该配体结合后未引起蛋白主链显著重排，整体稳定性较高。相比之下，吲哚-3-甲醇复合物在30.75?ns后波动增加，平均RMSD为2.28??，提示其结合后伴随更明显的全局构象变化，稳定性相对较弱，但并不意味着体系未平衡，而更可能反映较高的配体-蛋白柔性。

3.4. Analysis of Radius of Gyration (Rg)
回转半径分析用于反映复合物紧致性。吲哚基甲基硫代葡萄糖苷复合物在整个模拟期间平均Rg为24.119??，与apo蛋白和环丙沙星对照接近，提示体系保持较高紧致性。吲哚-3-甲醇复合物平均Rg为24.28??，在40–80?ns间出现轻微波动，说明其整体仍维持较紧凑构象，但稳定程度略逊于前者。

3.5. SASA
溶剂可及表面积（SASA）分析显示，吲哚基甲基硫代葡萄糖苷和吲哚-3-甲醇复合物的平均SASA分别为20,011.41??²和20,178.49??²，与apo蛋白19,751.86??²接近，表明配体结合并未显著改变蛋白整体溶剂暴露状态或引发大规模构象重塑。该结果支持两种配体与CDTa的结合总体较为温和且结构可耐受。

3.6. Analysis of Hydrogen Bond
氢键分析表明，吲哚基甲基硫代葡萄糖苷、吲哚-3-甲醇、apo蛋白和环丙沙星体系在模拟期间平均氢键数分别为347.83、342.84、350.80和339.01。两种候选配体的氢键网络均较稳定，提示其在动态环境下能够维持持续的配体-蛋白相互作用，为其潜在抑制活性提供了支持。

3.7. RMSF
RMSF（均方根涨落）用于描述残基水平柔性。结果显示，两种复合物在大部分残基区段的波动模式相近，主要在模拟起始区段如THR1、THR2、TYR3、LYS4、ALA5、PRO6和ILE7处观察到共同波动。作者指出，RMSF较低而RMSD相对较高的现象并不矛盾，这意味着蛋白局部残基柔性保持有限，但整体可能发生刚体运动或结构域位移。该解释尤其适用于吲哚-3-甲醇复合物，其较高RMSD更多反映全局构象改变，而非局部结构失稳。

3.8. Analysis of MM–PBSA
MM–PBSA自由能计算进一步区分了两种配体的结合优劣。吲哚基甲基硫代葡萄糖苷与CDTa复合物的结合自由能为?176.14?kJ/mol，说明相互作用具有明显热力学有利性。环丙沙星对照为?189.78?kJ/mol，也表现出较强结合。相比之下，吲哚-3-甲醇复合物的结合自由能为4.03?kJ/mol，提示在该模拟条件下其结合并不有利，支持其结合稳定性低于吲哚基甲基硫代葡萄糖苷。这一结果与前述RMSD和构象漂移分析相互印证。

3.9. Analysis of Superimposition
叠合分析比较了0、25、50、75和100?ns时点的构象变化。结果显示，吲哚基甲基硫代葡萄糖苷在整个模拟期间始终位于结合口袋内，0–50?ns阶段结构特别稳定，后续虽存在轻微相互作用变化，但总体结合持续存在。吲哚-3-甲醇在50和75?ns出现短暂脱离结合口袋现象，但在100?ns重新建立相互作用，提示其结合具有可逆性和动态性。该结果再次说明，吲哚基甲基硫代葡萄糖苷的复合物稳定性优于吲哚-3-甲醇。

3.10. Analysis of PCA
PCA结果显示，apo蛋白的PC1和PC2分别解释45.27%和14.84%的方差，表明未结合状态下CDTa构象柔性最高。环丙沙星复合物的PC1和PC2下降至20.25%和12.08%，显示配体结合后动态受限。吲哚基甲基硫代葡萄糖苷复合物的PC1和PC2为26.33%和19.19%，提示中度构象稳定化；吲哚-3-甲醇复合物PC1达36.64%、PC2为12.44%，表现为更明显的单一主方向构象偏移。总体而言，配体结合均改变了CDTa的构象运动方式，并对其动态自由度产生限制。

3.11. DCCM
DCCM分析进一步揭示了配体结合对残基间相关运动的影响。与apo状态相比，配体结合后的复合物在若干残基区间出现新的正相关或负相关运动模式，提示蛋白内部动力学耦联和变构通讯发生改变。吲哚基甲基硫代葡萄糖苷复合物在50–100和200–250残基区段周围显示增强的正相关运动，而吲哚-3-甲醇复合物在150–200残基区段部分削弱了apo状态下的正相关模式。这些结果支持配体结合能够重塑CDTa的动态行为，并可能影响其功能相关构象转换。

3.12. In Vitro Antibacterial Activity
在体外抑菌实验中，两种化合物均表现出对艰难梭菌的明确抑菌活性。吲哚基甲基硫代葡萄糖苷在50、75和100?μg/disc下的抑菌圈分别为11.33?±?0.33?mm、18?±?0.57?mm和23.33?±?0.67?mm；吲哚-3-甲醇分别为13.67?±?0.33?mm、16.16?±?0.16?mm和22.67?±?0.33?mm，均呈浓度依赖性增强。MIC/MBC测定进一步表明，吲哚基甲基硫代葡萄糖苷的MIC和MBC分别为10.33?±?0.72?μg/mL和23.33?±?1.36?μg/mL，略优于吲哚-3-甲醇的11.33?±?0.27?μg/mL和26.67?±?3.60?μg/mL。值得注意的是，环丙沙星虽在计算模拟中可与CDTa稳定结合，但在体外未显示抑菌圈，说明菌株可能存在氟喹诺酮耐药性，也提示CDTa结合本身不足以直接抑制细菌生长。作者据此明确区分了“抗菌效应”与“抗毒力效应”：CDTa作为分泌性毒力因子，并非细菌生存必需靶点，因此两种化合物的抑菌作用很可能还涉及膜完整性或代谢过程干扰，而针对CDTa的作用更可能构成补充性的抗毒力机制。

3.13. ADMET Analysis
ADMET结果显示，两种化合物在药代和毒理学特征上存在明显差异。吲哚-3-甲醇的水溶性、CaCO₂通透性及人肠道吸收率均优于吲哚基甲基硫代葡萄糖苷，后者人肠道吸收率仅5.082%，前者为92.085%。分布方面，吲哚-3-甲醇更易穿越血脑屏障（BBB），而吲哚基甲基硫代葡萄糖苷相对较弱。代谢方面，吲哚-3-甲醇可抑制CYP1A2，而吲哚基甲基硫代葡萄糖苷不抑制CYP1A2和CYP2D6，也不抑制P-糖蛋白，提示其药物相互作用风险较低。毒性方面，两者均不表现为AMES毒性，也不抑制hERG I/II；但吲哚基甲基硫代葡萄糖苷预测存在肝毒性，而吲哚-3-甲醇预测存在皮肤致敏性。类药性分析还显示，吲哚基甲基硫代葡萄糖苷存在Lipinski规则违例，口服生物利用度评分较低，仅0.11；吲哚-3-甲醇则无违例，生物利用度评分为0.55，药物开发属性更优。

3.14. Limitation of the Study
研究人员明确指出，本研究虽然通过分子对接、MD模拟和MM–PBSA计算提供了CDTa-配体相互作用的机制线索，但尚未通过微量热泳动（MST）、等温滴定量热（ITC）或核磁共振滴定等生物物理手段对直接结合进行实验确认。现有体外实验验证的是对艰难梭菌生长的抑制效果，而非对CDTa毒力功能或分泌过程的直接抑制。因此，后续仍需开展毒素功能实验和直接结合实验证实其抗毒力潜能。

综合讨论部分可见，本研究的核心贡献在于将天然产物筛选、分子模拟和体外抑菌验证相结合，提出了甘蓝来源吲哚类化合物兼具“抗菌+抗毒力”双重潜力的概念。结果表明，吲哚基甲基硫代葡萄糖苷在CDTa靶向稳定性、自由能优势以及体外抑菌指标方面整体优于吲哚-3-甲醇，是更值得优先推进的候选分子。不过，作者也谨慎指出，计算上的CDTa结合并不能等同于对细菌增殖的直接抑制，二者应在机制上区分理解。该研究的重要意义在于，为高毒力CDI的天然产物干预策略提供了一个兼顾毒力削弱和细菌抑制的研究框架，也为后续开展药代优化、毒素功能验证和体内疗效评估奠定了基础。

研究结论部分可译为：本研究围绕艰难梭菌二元毒素CDTa亚基（PDB ID: 6X41）与所选生物活性化合物之间的相互作用进行了系统考察，并开展了100?ns分子动力学模拟。通过RMSD、Rg和SASA等指标评估了环丙沙星、吲哚-3-甲醇和吲哚基甲基硫代葡萄糖苷与蛋白的结构稳定性及结合亲和性。在所考察配体中，吲哚基甲基硫代葡萄糖苷在整个模拟过程中表现出更优稳定性和持续紧凑构象，提示其与CDTa之间形成了强而稳定的相互作用。吲哚-3-甲醇也维持了有利相互作用，尽管伴有轻微构象偏移，但未破坏蛋白整体结构。与此同时，体外抗菌实验显示，吲哚基甲基硫代葡萄糖苷对艰难梭菌具有较强生长抑制活性，MIC和MBC结果进一步支持其内在抗菌效应。重要的是，预测得到的CDTa相互作用被认为代表一种补充性的抗毒力机制，可能有助于减弱毒素介导的致病性，而并非直接解释细菌生长抑制。总体而言，这些发现既揭示了配体-CDTa相互作用的机制特征，也凸显了B. oleracea来源化合物的抗菌属性。其中，吲哚基甲基硫代葡萄糖苷可作为治疗CDI的潜在候选分子，既可作为独立抗菌化合物，也可作为靶向毒力的辅助治疗分子。未来仍需通过药代动力学研究和体内验证进一步评估其疗效、安全性及转化潜力。