Shorea parvifolia根皮中未揭示寡聚白藜芦醇作为黄嘌呤氧化酶（XO）抑制剂的分离、体外活性、量子化学与结构计算研究

《RSC Advances》：Unrevealed oligoresveratrols from Shorea parvifolia root bark as xanthine oxidase inhibitors: isolation, in vitro bioassay, quantum studies and structure-based computational profilingOpen Access

　　

白藜芦醇寡聚物（resveratrol oligomers）来源于Shorea物种，具有多样生物活性，但决定其抑制黄嘌呤氧化酶（xanthine oxidase, XO）的结构因素仍认识不足。本研究调查了从Shorea parvifolia根皮中分离的寡聚白藜芦醇（oligoresveratrols），并考察其结构与XO抑制潜力之间的关系。研究人员分离了五种寡聚白藜芦醇，并通过光谱分析进行结构表征。研究人员评估了它们的XO抑制活性，随后对最活跃化合物进行了药代动力学（pharmacokinetics）和类药性（drug-likeness）预测、密度泛函理论（density functional theory, DFT）计算、分子对接（molecular docking）和分子动力学（molecular dynamics, MD）模拟。研究鉴定了五种化合物，包括isoampelopsin D（1）、diptoindonesin G（2）、shorealactone（3）、α-viniferin（4）和neohopeaphenol A（5）。XO抑制在化合物间差异显著，shorealactone（3）显示最强活性（IC₅₀ = 32.53 μM），优于临床参考标准allopurinol（40.1 μM）。Isoampelopsin D（1）也显示强抑制（65.15 μM），而diptoindonesin G（2）在该系列中活性最弱（161.08 μM）。分子对接识别出化合物1–3在XO钼蝶呤结构域（molybdopterin domain）内的有利相互作用，MD在100 ns模拟中提供化合物稳定性。MM-PBSA计算得到化合物2、1和3的有利结合自由能，ΔG_bind分别为?18.03、?17.17和?14.36 kcal mol^?1，而allopurinol为?9.9972 kcal mol^?1。然而，化合物3达到最有利热力学亲和力（ΔG_exp = ?6.12 kcal mol^?1），优于allopurinol（?6.00 kcal mol^?1）。其后为化合物1（?5.71 kcal mol^?1）和化合物2（?5.17 kcal mol^?1）。总体而言，这些不同的生物活性谱确定isoampelopsin D（1）和shorealactone（3）为未来开发天然XO抑制剂的有前景结构骨架。

研究背景与目的：Dipterocarpaceae科植物在东南亚低地热带雨林中具有重要生态与化学多样性，Shorea属尤其是Shorea parvifolia具有广泛分布。既往植物化学研究发现Shorea属富含白藜芦醇寡聚物（resveratrol oligomers, ROs），这类化合物由C₆–C₂–C₆单元的多种区域选择性和立体选择性氧化偶联形成，包括二氢苯并呋喃、茚满和双环骨架等，并受氧化、重排和杂合偶联进一步修饰。ROs已报道具有抗氧化、抗炎、抗癌、抗菌和抗糖尿病等活性，但其抑制XO的结构决定因素仍不清楚。XO是可催化次黄嘌呤和黄嘌呤氧化为尿酸并产生超氧阴离子自由基（O₂^˙?）和过氧化氢（H₂O₂）的金属酶，是氧化应激的重要内源性来源；多酚可通过清除自由基或结合并抑制XO催化域发挥作用。因此，本研究旨在从S. parvifolia根皮分离寡聚白藜芦醇，评价其XO抑制活性，并通过药代/类药性、DFT、分子对接和MD解析结构–活性关系，为发现天然XO抑制剂提供依据。

关键技术方法：植物材料为2023年7月采自印度尼西亚中加里曼丹West Kotawaringin Regency、Kumai Subdistrict、Keraya Village的Shorea parvifolia根皮，由Generasi Biologi Indonesia Foundation鉴定（证书号372/02.Genbinesia/2023）。研究人员以丙酮提取，经硅胶柱层析、Sephadex LH-20、制备TLC和半制备HPLC分离化合物，以NMR及文献比对鉴定结构；采用体外XO抑制实验测定IC₅₀；运用SwissADME和pkCSM进行药代与类药性预测；以DFT/B3LYP/6-311++G(d,p)计算FMO和MEP；以DOCK6对牛乳XO（PDB ID: 1FIQ）进行分子对接并以SAL重对接验证；以GROMACS 2021进行100 ns MD模拟及MM-PBSA结合自由能计算。

研究结果：分离化合物的结构鉴定：研究人员从丙酮根皮提取物中分离得到1–5，并经¹H、¹³C NMR、HSQC、HMBC、NOESY及文献比对鉴定。1为isoampelopsin D（C₂₈H₂₂O₆），属修饰茚型二聚体，含亚甲基桥；其C-7a为sp²碳（δ_C 150.4），C-7b为sp³亚甲基（δ_C 30.8，δ_H 3.85），NOESY支持C-8a处3,5-二羟基苯基的R相对构型；该化合物此前仅作为合成衍生物已知，本研究首次从天然来源分离。2为diptoindonesin G（C₂₁H₁₂O₆），完全氧化平面菲型二聚体；3为shorealactone（C₃₄H₂₈O₁₂），白藜芦醇–抗坏血酸杂合加合物，含三环四氢呋喃内酯；4为α-viniferin（C₄₂H₃₀O₉），三聚体；5为neohopeaphenol A（C₅₆H₄₂O₁₂），对称四聚体。结论是获得结构多样、聚合度不同的寡聚白藜芦醇，并确认1的天然存在。

黄嘌呤氧化酶生物测定：体外XO实验显示，3最强（IC₅₀ = 32.53 μM），优于allopurinol（40.1 μM）；1为65.15 μM，接近allopurinol；5为106.96 μM，4为143.57 μM，2最弱（161.08 μM）。结论是XO抑制差异与寡聚体尺寸、空间位阻及催化腔地形相关，高阶三/四聚体难以深入狭窄疏水通道，低分子量且构象适配的3和1更有利；后续选择1–3进行计算，2作为结构对照，4和5因空间位阻排除。

分离化合物的药代动力学预测：SwissADME雷达和pkCSM预测显示，1和2在脂溶性、分子尺寸、极性、溶解度、柔性方面较优，但因共轭多酚骨架使不饱和度超出（Fraction C-sp³ = 0.07）；3因饱和内酯环使Fraction C-sp³ = 0.26，但分子量（628.59 g mol^?1）和极性表面积（TPSA = 195.60 ?²）超出；三者log P_o/w为2.26–4.68，人肠道吸收85.60–100%，血脑屏障渗透低。结论是1–3具外周靶向可行性；3和1代谢安全性优于2，2有CYP2C19/CYP2D6抑制和潜在肝毒性；3在吸收、代谢稳定、非肝毒性等方面平衡最佳，1也具良好类药性，2虽Lipinski/Veber无违背但代谢负担较高。

结构优化与量子性质：DFT/B3LYP/6-311++G(d,p)计算显示，2的前线轨道能隙最窄（ΔE = 3.4109 eV），化学硬度最低、亲电性最高，与其刚性共轭平面菲型骨架一致；1为中间能隙（4.2387 eV）；3能隙最宽（4.9634 eV），化学硬度最高、稳定性较高，归因于非共轭sp³抗坏血酸四氢呋喃/内酯部分限制π离域。全局反应描述符与MEP进一步支持：2在C8a羰基氧附近具最负电势（?7.513×10^?2 a.u.），1酚羟基质子附近具正电势（6.471×10^?2 a.u.），3的内酯/醚氧形成负电区并由七个羟基质子形成正电区。结论是电子性质与结构共轭程度密切相关，并影响活性位点相互作用。

分子对接分析：以DOCK6对牛乳XO（PDB ID: 1FIQ）进行柔性对接，SAL重对接RMSD为0.1916 ?，验证网格参数可靠；关键活性位点残基包括Glu233、Leu304、Ser307、Arg311、Phe345、Ser439、Phe440、Thr441、Val442、Leu445、Ala509和Ala510。Grid score顺序为3（?37.01 kcal mol^?1）<1（?34.29）
分子动力学模拟：100 ns MD显示，apo、SAL、allopurinol、1和2的蛋白骨架RMSD稳定在0.25–0.38 nm；3在前20 ns升至约0.65 nm，随后平衡于0.52–0.58 nm，体现其大体积内酯基团的诱导契合重定向。催化残基RMSF低于0.15 nm，表面环为0.3–0.5 nm；配体结合使催化腔局部刚性化。R_g总体稳定，3略高（约2.86–2.90 nm）但未展开；SASA稳定在约300–325 ?²；总氢键数保持520–600。配体RMSD显示2最稳定（0.03–0.08 nm），1为0.1–0.225 nm，3在50–55 ns发生重定向后进入平台；配体–蛋白氢键方面，1为1–4个（60 ns后3–5个），2为0–2个（主要1个），3为1–3个（峰值4个）。MM-PBSA给出ΔG_bind：2为?18.03、1为?17.17、3为?14.36 kcal mol^?1，均优于allopurinol的?9.9972 kcal mol^?1；由IC₅₀推算的ΔG_exp显示3最有利（?6.12 kcal mol^?1），优于allopurinol（?6.00），其后为1（?5.71）和2（?5.17）。结论是1–3均可稳定结合XO活性位点，3虽经历构象重定向但最终形成稳定复合物，且其水相热力学效率与体外强抑制一致。

讨论部分总结：讨论认为，XO识别受寡聚白藜芦醇骨架空间氧化模式影响；3的多氧内酯官能团以及1的茚核和酚羟基可在钼蝶呤结构域结合腔内形成有利构象；高阶三/四聚体受活性位点空间位阻限制。药代与类药性预测提示1和3较优，3类药性特征突出，2虽有良好Lipinski特征但存在CYP抑制和肝毒性风险。DFT、对接、MD和MM-PBSA从电子性质、结合姿态、动态稳定性和热力学方面支持体外活性结果。研究结论翻译：本研究确定Shorea parvifolia根皮为结构独特寡聚白藜芦醇（1–5）的来源，并报道isoampelopsin D（1）的首次天然存在，其此前仅作为合成衍生物已知。分离化合物为天然多酚衍生物，具有显著不同的XO抑制能力。Shorealactone（3）具有最强XO抑制作用（IC₅₀ = 20.45 μg mL^?1，与allopurinol的20.39 μg mL^?1相当），isoampelopsin D（1）紧随其后（29.61 μg mL^?1），而大体积三聚体α-viniferin（4）和四聚体neohopeaphenol A（5）受活性位点空间位阻限制。结构和计算分析提示，XO识别受寡聚白藜芦醇骨架空间氧化模式影响，化合物3的多氧内酯官能团以及化合物1的茚核和酚羟基在钼蝶呤结构域结合腔内形成有利构象；同时化合物1和3也表现出良好药代动力学预测，其中化合物3具有显著类药性特征。MM-PBSA计算得到化合物2、1和3的有利结合自由能，ΔG_bind分别为?18.03、?17.17和?14.36 kcal mol^?1，而allopurinol为?9.9972 kcal mol^?1。然而，化合物3达到有利热力学亲和力（ΔG_exp = ?6.12 kcal mol^?1），其后依次为化合物1（?5.71 kcal mol^?1）、allopurinol（?5.22 kcal mol^?1）和化合物2（?5.17 kcal mol^?1）。这种差异提示isoampelopsin D（1）和shorealactone（3）是进一步研究天然XO抑制剂的有前景骨架。