香芹酚(Carvacrol)与百里香酚(Thymol)烷基化衍生物靶向WSSV及致AHPND副溶血弧菌(Vibrio parahaemolyticus)的分子与结构机理研究

《Archives of Microbiology》：Molecular and structural insights into carvacrol and thymol alkylated derivatives targeting WSSV and AHPND-causing Vibrio parahaemolyticus

　　

对虾养殖业集约化发展使环境不稳定相关的传染性疾病发生率上升，尤以白斑综合征病毒(White Spot Syndrome Virus, WSSV)和由副溶血弧菌(Vibrio parahaemolyticus)引起的急性肝胰腺坏死病(Acute Hepatopancreatic Necrosis Disease, AHPND)为重，导致高死亡率和重大经济损失。鉴于现有治疗手段效果有限，寻找天然源生物活性化合物成为创新且可持续的策略。本研究通过计算筛选，考察香芹酚(carvacrol)与百里香酚(thymol)衍生物对靶蛋白PirB(V. parahaemolyticus)及胸苷酸合成酶(thymidylate synthase, WSSV wTS)的作用。腔体分析结果表明，衍生物与PirB结合能有利且复合物构象稳定性高，提示强抗菌抑制潜力；虽wTS亲和力参数表明抗病毒优先级较低，但与共结晶配体的持续重叠支持其具相关生物活性的合理性。值得注意的是，香芹酚和百里香酚衍生物CVC及TMC化合物表现更优，具更高结构稳定性和相互作用稳健性。综上，香芹酚和百里香酚衍生物——特别是CVC和TMC——因在评估靶点上兼具适中结合亲和力和一致的结构稳定性，是颇具前景的候选分子，为对虾养殖中新治疗策略的合理开发提供计算框架，有助于更高效和可持续的方法，但尚需进一步实验验证。

论文解读：香芹酚与百里香酚烷基化衍生物靶向WSSV及致AHPND副溶血弧菌的分子与结构机理研究

对虾（主要是凡纳滨对虾Litopenaeus vannamei）集约化养殖中，环境波动诱发的高发性传染病——白斑综合征病毒（White Spot Syndrome Virus, WSSV，属Nimaviridae科，双链DNA病毒）与急性肝胰腺坏死病（Acute Hepatopancreatic Necrosis Disease, AHPND，主要由携带pirA/B毒素质粒的副溶血弧菌Vibrio parahaemolyticus引起）——造成养殖群体死亡率超65%，全球累计经济损失达数十亿美元。目前尚无感染后有效抗病毒疗法，抗生素滥用亦带来耐药性与残留问题。WSSV复制依赖病毒胸苷酸合成酶（viral Thymidylate Synthase, wTS，PDB: 7XYJ）催化dUMP→dTMP以供给dTTP合成；AHPND致病依赖PirB毒素蛋白（PDB: 3×0U）介导肝胰腺损伤。天然酚类单萜烯——香芹酚（carvacrol, 5-异丙基-2-甲基苯酚）与百里香酚（thymol, 2-异丙基-5-甲基苯酚）——具广谱抗细菌及膜干扰活性，但其烷基化卤代衍生物对wTS和PirB的分子水平抑制作用尚未明确。本研究以计算药物设计手段，评估香芹酚与百里香酚烷基化衍生物（CVA、CVB、CVC系列及TMA、TMB、TMC系列）对wTS和PirB的结合潜力、动态稳定性及生态毒理学风险，为水产养殖用新型植物源防治剂开发提供理论依据。本文发表于《Archives of Microbiology》。

研究人员采用的主要关键技术方法如下：以PASS Online平台预测香芹酚/百里香酚及其衍生物的生物活性谱（Pa/Pi值）；用MarvinSketch构建烷基化衍生物结构，Avogadro（MMFF94力场）优化配体构象；从RCSB PDB获取wTS（PDB 7XYJ，共结晶配体dUMP/uMP）及PirB（PDB 3×0U，参考抑制剂庆大霉素Gentamicin/GEN）三维结构，UCSF Chimera去水及去除原配体后加极性氢；AutoDock Vina（exhaustiveness=64）进行分子对接(molecular docking)，以共结晶配体重对接(redocking)验证（RMSD≤2 ?为成功），选取最佳pose按结合自由能ΔG计算抑制常数K_i；Pharmit平台做药效团(pharmacophore)建模比对；iMODS服务器基于简正模分析(Normal Mode Analysis, NMA)做蛋白质-配体复合物分子动力学模拟，分析Cα RMSD、可变形性(deformability)及B因子；JANUS软件基于QSAR/QSPR模型预测logK_ow、水溶性、各环境区室持久性及鱼/Daphnia/藻类急慢性毒性(LC₅₀/EC₅₀/ChV)。

Bioactivity in carvacrol/thymol derivatives（香芹酚/百里香酚衍生物的生物活性）

经PASS Online预测，母体化合物及衍生物抗细菌Pa值偏低（香芹酚0.319，百里香酚0.336），系模型难以捕捉膜破坏等非特异性机制所致；部分百里香酚衍生物抗病毒Pa≈0.53。结果表明结构修饰影响预测活性谱，但属初步筛选，不替代实验验证。

Correlation between inhibition potential and affinity energy（抑制潜能与亲和能的相关性）

Inhibition potential and affinity energy via V. parahaemolyticus（副溶血弧菌PirB靶点）：对接显示参考药GEN结合能?7.3 kcal/mol（pK_i=5.46），CVC结合能?6.3 kcal/mol（pK_i=4.64），优于其他衍生物（?5.7至?6.0 kcal/mol），所有复合物RMSD<2 ?。CVC与PirB结合能最负（除GEN外），提示较强预期相互作用。

pK_ivs. E_Avia wTS（wTS靶点）：共结晶配体uMP结合能?7.6 kcal/mol（pK_i=5.58），酚类衍生物结合能?5.2至?5.9 kcal/mol（pK_i=4.08–4.41），RMSD多<1.7 ?。CVC(?5.9 kcal/mol, pK_i=4.41)在衍生物中结合略优，但整体亲和力弱于uMP，提示较低抗病毒优先级。

Molecular docking results（分子对接结果）

Redocking and structural evaluation of protein–ligand interactions between uMP and GEN（uMP与GEN重对接及评价）：uMP–wTS重对接RMSD低，可重现原结合模式（H171/S191/Y233氢键及精氨酸盐桥）；GEN–PirB结合能?7.3 kcal/mol，与Y201/R236形成氢键及疏水堆积，验证对接流程可靠。

Protein–ligand interaction analysis of carvacrol/thymol derivatives with V. parahaemolyticus（衍生物与PirB相互作用分析）：香芹酚仅靠疏水作用（I130/V133/L193）；CVA多氢键+Y201疏水；CVB氢键(D216/R236)+π-π堆积(F244)；CVC形成深疏水接触(F244/M247/L267)+π-π堆积(F244)，结合最稳定。百里香酚靠疏水+I130/I138/L193及Y192/N196氢键；TMA仅疏水+π-π堆积；TMB多氢键但无显著能量增益；TMC（结合能同组最优）具深疏水接触(F244/M247/L267)+π-π堆积(F244)，腔重叠与GEN部分一致，暗示潜在抑菌作用。

Interaction profiling of carvacrol/thymol derivatives with wTS（衍生物与wTS相互作用分析）：衍生物均占据wTS活性中心，主要靠疏水/π-π堆积（F200A/L196A/I83A/W84A），与原配体uMP（多重氢键+盐桥）不同；CVC和TMA与关键残基维持较稳定接触。

Pharmacophore-based analysis（药效团分析）

香芹酚衍生物与PirB作用以疏水接触（I83/L196/F200）及π-π堆积（F200）为主，芳香中心对稳定复合物关键；百里香酚衍生物呈相似模式（结构异构导致药效团特征相近）。CVC/TMC因卤代+脂肪链增强疏水接触与空腔适应，结合略优。配体与uMP占位区部分重叠，理论上具竞争/协同阻断底物进入可能。

NMA-based molecular dynamics Results（基于NMA的分子动力学结果）

MD V. parahaemolyticus between carvacrol/thymol derivatives（PirB体系）：PirB–CVC与PirB–TMC复合物Cα RMSD终值分别为1.41 ?与1.57 ?（<2 ?），可变形性峰接近GEN水平，配体未引起异常柔性改变，表明CVC/TMC与PirB形成结构稳定复合物。

MD wTS between carvacrol/thymol derivatives（wTS体系）：wTS–CVC与wTS–TMA复合物Cα RMSD终值分别约1.16 ?与1.18 ?，柔性分布与游离蛋白相似，结合未破坏蛋白骨架完整性，与对接预测一致。

Ecotoxicological analysis of carvacrol/thymol derivatives（生态毒理分析）

logK_ow介于3.38（百里香酚）–4.01（CVA/TMA），属中等亲脂；CVC/TMC因氯代长链烷基致水溶性剧降（≈3.04–3.05 mg/L）且沉积物持久性>180天（最长≈987天），水体半衰期≤19天，土壤≤70天。急/慢性毒性预测显示CVC与TMC对鱼类LC₅₀<1 mg/L、对Daphnia LC₅₀<1 mg/L、对藻类EC₅₀极低，ChV对鱼亦达高毒等级；其余衍生物多为中毒。提示卤代修饰增加环境持久性与水生毒性风险，需谨慎评估。

讨论与结论（Conclusion部分翻译浓缩）

高级计算分析表明香芹酚与百里香酚衍生物可有效作用于抗菌（PirB）及抗病毒（wTS）相关靶蛋白。CVC与TMC在抗菌靶点上展现高结构稳定性与最小构象变化；在抗病毒靶点占据与天然配体重叠的疏水区，提示潜在空间位阻阻断底物进入。生态毒理预测显示CVC/TMC因卤代致持久性与毒性升高，但该变化部分源于预测模型对卤代酚类的外推局限。综上，结构修饰的酚类单萜烯——特别是CVC（香芹酚氯代烷基衍生物）和TMC（百里香酚氯代烷基衍生物）——具作为对虾养殖中抗WSSV/AHPND生物活性候选分子的潜力。鉴于本研究为纯计算模拟，需后续体内外实验确证其调控抗菌与抗病毒活性的能力，从而为控制对虾WSSV与AHPND的治疗策略提供依据。本研究为后续验证CVC/TMC生物活性及探索增强分子相互作用的新结构修饰奠定计算基础。