溶剂效应对Rhodocoranes F-J分子性质、理化参数、非线性光学(NLO)行为及分子对接影响的计算研究

《Next Materials》：Insights into solvent effects on molecular properties, physicochemical parameters, NLO behaviour, and Molecular docking of Rhodocoranes F-J: A computational study

　　

摘要：本研究专门针对Rhodocoranes（F-J）——一类已知具抗真菌与抗菌活性的倍半萜类化合物——进行计算机(in-silico)研究。该类化合物从掌状射纹菌(Rhodotus palmatus)蘑菇中分离得到，此前未被鉴定且显示药物应用潜力。研究重点之一为考察不同溶剂对电子光谱的影响。含时密度泛函理论(TD-DFT)研究表明，除少数情况外，最高波长处电子跃迁主要来自HOMO→LUMO，对多数化合物贡献达98%。溶剂极性可稳定分子——Rh-F的HOMO能级由庚烷(heptane)中的?6.09 eV移至水中变为?6.18 eV。Rh-G显示最高化学软度(softness, σ=0.565)及最低化学硬度(hardness, η=1.767 eV)。非线性光学(NLO)分析表明，各化合物的一阶超极化率(βtot)为标准尿素(urea)的6–9倍，其中Rh-G最高达3.69×10?30esu。分子对接重点评估化合物与19种蛋白的结合亲和力，以筛选黑色素瘤、乳腺癌及宫颈癌潜在治疗药物——Rh-F对皮肤癌蛋白(?8.4 kcal/mol)、乳腺癌蛋白(?8.7 kcal/mol)及细菌蛋白(?8.3 kcal/mol)具显著结合亲和力；Rh-G对宫颈癌蛋白结合能最优(?8.7 kcal/mol)。此外，借助ADME及ProTox-II进行毒性预测以完善安全性评价，Rh-H、I、J的LD50=1000 mg/kg，属第4类毒性(class 4 toxicity)，提示较低毒性水平。

论文解读：溶剂效应对Rhodocoranes F-J分子性质、理化参数、NLO行为及分子对接影响的计算研究

Rhodocoranes（F-J）是从掌状射纹菌(Rhodotus palmatus) submerged cultures中分离得到的一类新型倍半萜(sesquiterpenoid)及降倍半萜(norsesquiterpenoid)类化合物（Rhodocoranes F-L），前期研究仅报道其轻微抗真菌活性及全合成路线，尚未对其开展系统的量子化学计算、非线性光学(nonlinear optical, NLO)性质预测、类药物性质评价及多靶点抗肿瘤/抗菌分子对接筛选。鉴于天然产物开发中早期in-silico筛选可大幅降低试错成本，研究人员选取文献已报道结构的Rhodocoranes F-J为对象，综合溶剂效应下的电子结构分析、前线分子轨道(frontier molecular orbital, FMO)参数、NLO响应、ADME-Tox及针对皮肤癌（黑色素瘤）、乳腺癌、宫颈癌及相关致病菌靶蛋白的分子对接，全面评估其作为光子学材料与候选药物的潜力。该论文发表于《Next Materials》。

主要关键技术方法： 研究人员以文献中Rhodocoranes F-J结构为初始模型，用ChemDraw构建并经Gaussian 09W依次以半经验、Hartree-Fock及B3LYP/6–311++G(d,p)基组优化几何；采用含时密度泛函理论(time-dependent DFT, TD-DFT)结合Conductor-like Polarizable Continuum Model(CPCM)模拟庚烷、甲苯、乙醇、甲醇、水五种溶剂中的紫外–可见吸收；用优化后结构计算偶极矩、平均极化率(α)及一阶超极化率(β_tot)；依Kohn-Sham轨道能计算电离能(ionization potential, IP)、电子亲和能(electron affinity, EA)、化学硬度(η)、化学软度(σ)、电负性(χ)、亲电指数(ω)；用SWISSADME预测类药物性及Lipinski类药五规则(Rule of Five, Ro5)，用ProTox-II/ProTox-3.0预测急性毒性类别(LD₅₀)及肝毒性、致突变性等端点；从RCSB PDB下载19个癌症/细菌靶蛋白，经UCSF Chimera 1.18处理（去水、加氢）后用AutoDock Vina进行分子对接，并以CB-DOCK与Achilles Blind Docking Server交叉验证，选最低结合自由能复合物用iMODS服务器做正则模分析(normal mode analysis, NMA)评估复合物稳定性。

3.1. Solvent Effects（溶剂效应—UV–Vis与前线轨道）

研究人员用TD-DFT/CPCM计算五种溶剂中Rhodocoranes F-J的第一激发态吸收波长(λ_max)、激发能与振子强度(oscillator strength, f)。结果表明，除Rh-G含HOMO?1→LUMO混合跃迁外，其余化合物最长波长吸收均源自HOMO→LUMO（贡献98%）；Rh-H/I/J在甲苯中出现高色效应(hyperchromic shift)且振子强度最大、激发能最低。随溶剂介电常数(ε)增大，HOMO能级下降（如Rh-F由庚烷?6.09 eV至水?6.18 eV），表明极性溶剂稳定HOMO；LUMO能级基本不受溶剂ε影响；HOMO–LUMO能隙(ΔE_g)随溶剂极性增大略减小，Rh-G降幅最明显，说明高极性溶剂中分子反应活性增强。

3.2. Frontier molecular orbital analysis（前线分子轨道分析）

研究人员依气相B3LYP/6–311++G(d,p)计算ΔE_g=E_LUMO?E_HOMO，得ΔE_g顺序为Rh-F(4.500 eV)＞Rh-I(4.288 eV)＞Rh-H(4.285 eV)＞Rh-J(4.261 eV)＞Rh-G(3.534 eV)。Rh-F具最大动力学稳定性与最高化学硬度，Rh-G具最小能隙与最高化学软度(σ)及亲电性(electrophilicity, ω)，窄能隙利于电子跃迁，预示更高化学反应性与NLO响应。

3.3. Global and local reactivity parameters（全局与局部反应性参数）

研究人员据IP与EA计算全局描述符：Rh-F硬度η=2.250 eV最高，Rh-G η=1.767 eV最低、σ=0.565最高、ω=44.23最高；Rh-F/I/J亲电指数较低。分子静电势(molecular electrostatic potential, MEP)图谱显示Rh-F的氧原子区域（红色/黄色）为亲电位点，蓝色/绿色区为亲核中心，与蛋白结合位点互补。

3.4. Non-linear optics（非线性光学性质）

研究人员在气相下计算静态偶极矩(μ)、平均线性极化率(α)及一阶超极化率(β_tot)。尿素参比β_tot=0.3728×10^?30esu，Rhodocoranes F-J的β_tot为其6–9倍，其中Rh-G达3.69×10^?30esu（约9.9倍），归因于非中心对称结构、共轭π电子离域及邻位羟基–吸电子乙酰基促进分子内电荷转移(intramolecular charge transfer, ICT)，β_xxx分量主导总超极化率，表明电荷转移主要沿x轴。

3.5. ADME（吸收、分布、代谢、排泄预测）

研究人员用SWISSADME评估类药物性。所有Rhodocoranes F-J满足Lipinski Ro5、Ghose、Egan及Veber规则，生物利用度评分(bioavailability score)=0.55；雷达图除sp³碳分数略超限外其余参数落于类药物允许区(pink region)。共识log P_o/w介于2.48–3.29，提示适中亲脂性；ESOL与Ali模型预测水溶性良好，SILICOS-IT模型仅Rh-F/G显示中等水溶性。

3.6. Toxicity investigation（毒性预测）

研究人员用ProTox-3.0预测急性毒性：Rh-H/I/J LD₅₀=1000 mg/kg，属第4类毒性(class 4 toxicity, 500–5000 mg/kg)，提示低急性毒性；Rh-G被划为class 1（高毒），其余毒性端点（ hepatotoxicity、mutagenicity、immunotoxicity等）文中仅述及筛选保留class 4及以上者进一步考量。

3.7. Normal mode analysis of the docked complex（对接复合物的正则模分析）

研究人员对四个最优对接复合物（Rh-F–5ZC9、Rh-F–1G50、Rh-G–4EOJ、Rh-F–3V99）做NMA。各复合物特征值(eigenvalue)极低（10^?6~10^?8量级），方差(variance)较高（30%–55%），协方差矩阵呈明暗区对应残基相互作用，弹性网络模型示靶蛋白刚性低、易变形但形成稳定配体–受体复合物，证实Rh-F与皮肤癌蛋白5ZC9、乳腺癌1G50、细菌3V99，及Rh-G与宫颈癌4EOJ间形成稳定结合。

3.8. Molecular docking studies（分子对接研究）

研究人员将Rh-F～J对接19个癌症/细菌靶蛋白，以AutoDock Vina结合同期CB-DOCK与Achilles Blind Docking验证。筛选出最优结合：Rh-F与黑色素瘤相关蛋白5ZC9结合能?8.4 kcal/mol（氢键GLN115、ASP337，疏水PHE366/PHE52/VAL119/ARG365）；与雌激素受体相关乳腺癌蛋白1G50 ?8.7 kcal/mol（氢键ARG1477/GLN1470/ASN1304，疏水HIS1474/LEU1370）；与革兰氏阴性菌靶3V99 ?8.3 kcal/mol（氢键ARG370/ARG457/TYR，疏水ALA453/LEU244/VAL243/ILE454/PHE450）。Rh-G与宫颈癌相关蛋白4EOJ ?8.7 kcal/mol（氢键ASN132/ASP145，疏水ALA31/LEU134/PHE80/VAL64/VAL18/LYS33）。

结论(Conclusion)翻译：

TD-DFT计算表明Rh-H、I、J因扩展共轭在甲苯中HOMO–LUMO能隙缩小，与实验趋势一致。NLO研究显示Rh-G为一阶超极化率β_tot=3.69×10^?30esu的高性能NLO材料，约尿素9倍，适用于光子学器件。类药物参数确认所有分子符合Lipinski规则，生物利用度评分0.55且具较好水溶性。分子对接与MD模拟指出Rh-F为多靶点先导化合物，对皮肤癌(?8.4 kcal/mol)与乳腺癌(?8.7 kcal/mol)蛋白及细菌靶蛋白(?8.3 kcal/mol)具强结合；Rh-G对宫颈癌蛋白结合能?8.7 kcal/mol最优。Rh-F、H、I、J毒性分类为第4类(LD₅₀=1000 mg/kg)，安全性可接受。综上，Rhodocorane-F被确立为极具前景的候选药物分子，值得进一步药学实验验证。