《The Journal of Physical Chemistry B》：Free-Energy and Position-Dependent Diffusion Profiles Shape Passive Permeation of a Semisynthetic Naphthoquinone across a Model Lipid Bilayer
编辑推荐：
被动膜渗透既取决于溶质分配，也取决于坐标依赖性迁移率，但二者各自及其空间耦合的贡献很少在统一的阻力框架内被量化。研究人员采用四组分模型脂质双层，研究了一种在临床前模型中具有低微摩尔细胞毒性和肿瘤生长抑制活性的半合成萘醌CNFD（6b,7-dihydro-5H-
被动膜渗透既取决于溶质分配，也取决于坐标依赖性迁移率，但二者各自及其空间耦合的贡献很少在统一的阻力框架内被量化。研究人员采用四组分模型脂质双层，研究了一种在临床前模型中具有低微摩尔细胞毒性和肿瘤生长抑制活性的半合成萘醌CNFD（6b,7-dihydro-5H-cyclopenta[b]naphtho[2,1-d]furan-5,6(9aH)-dione）。研究人员结合了成熟的模拟退火增强伞形采样（simulated tempering-enhanced umbrella sampling, ST-eUS）、加权直方图分析法（weighted histogram analysis method, WHAM）、Hummer位置自相关形式体系以及非均匀溶解度-扩散（inhomogeneous solubility-diffusion, ISD）模型。三个独立重复产生了有利的界面分配（?4.66 kJ mol?1）和正的膜中心自由能（+12.36 kJ mol?1）。受控剖面替换表明，PMF异质性和位置依赖性扩散分别将CNFD阻力放大了10.54和14.89倍，而其残余空间耦合因子为1.134。PMF最大值与扩散最小值相距2.14 nm。主要内在有效渗透系数（Peff）为9.09 × 10?2 cm s?1。将相同分析应用于β-拉帕醌（β-lapachone）和白花丹醌（plumbagin），识别出对比鲜明的阻力机制，其中白花丹醌结合了有利的中心PMF和低得多的中心阻力局域化。这些结果表明，空间分离的自由能景观和迁移率景观如何决定内在脂质相渗透。
半合成萘醌被动渗透的物理机制：基于ST-eUS/WHAM、Hummer自相关与ISD模型的定量分解
**研究背景与问题**
被动膜渗透是决定药物生物利用度、分布和疗效以及环境外源性物质毒理学特征的关键环节。传统的脂水分配规则虽然提供了经验性指导，但对结构复杂分子的渗透系数（P
eff）进行定量预测仍然困难，因为热力学分配和动力学屏障之间存在复杂耦合。萘醌类化合物在药物化学中具有重要地位，但其膜渗透的分子机制尚未被充分理解。实验方法如Caco-2实验和PAMPA只能给出整体渗透估计，无法提供原子级分辨率的渗透途径和能量信息。分子动力学（MD）模拟可以显式模拟溶质-脂质-水相互作用，其中非均匀溶解度-扩散（ISD）模型将P
eff分解为平均力势（PMF）和位置依赖性扩散系数D
⊥(z)的积分。然而，PMF异质性、扩散异质性以及二者空间排列对膜阻力的独立贡献很少在同一框架内被定量比较。针对此问题，研究人员选择半合成萘醌CNFD，在四组分模型脂质双层中结合多种增强采样和扩散分析手段，精确量化这些贡献。该研究发表在《The Journal of Physical Chemistry B》。
**主要技术方法**
研究人员利用密度泛函理论（DFT）在B3LYP/6-311++G(d,p)水平优化CNFD几何结构，并使用CGenFF参数化。采用CHARMM-GUI构建由POPC、POPE、胆固醇和PSM组成的四组分双层膜系统（共51,966个原子）。通过恒速拉伸分子动力学（SMD）从水相拉向膜中心，生成46个伞形采样窗口。模拟退火增强伞形采样（ST-eUS）采用9个温度态（310.15–374.15 K），对每个窗口进行三副本、每副本20 ns的生产模拟。PMF重建采用加权直方图分析法（WHAM），位置依赖性扩散系数D
⊥(z)通过Hummer位置自相关形式由约束轨迹估算（数值积分作为主要估计，拟合PACF作为敏感性分析），最终代入ISD模型计算内在脂质相P
eff。CNFD、β-拉帕醌（LAP）和白花丹醌（PLU）均采用相同的46窗口、三副本工作流。
**研究结果**
**3.1 界面分配与中心自由能垒**
通过ST-eUS/WHAM重建CNFD在半个对称双层中的PMF，以体相水为参考。结果显示界面处存在?4.66 kJ mol
?1的极小值，而膜中心为+12.36 kJ mol
?1的正值。这表明CNFD在界面相对水相热力学有利，但进入双层中心是不利的。
**3.2 位置依赖性扩散的独特最小值**
Hummer位置自相关分析138条约束轨迹（46窗口×3副本）得到强位置依赖的D
⊥(z)。数值自相关积分显示异质性跨度（D
max/D
min）为57.9倍，其中最小值D
min = 2.79 × 10
?7 cm
2 s
?1位于z = 2.20 nm，最大值D
max = 1.62 × 10
?5 cm
2 s
?1位于z = 4.41 nm。拟合PACF敏感性分析显示约39.8倍跨度。PMF最大值在膜中心，而D
⊥(z)最小值更靠近头基-核心转变区，二者相距2.14 nm。
**3.3 ISD积分的内在渗透率**
对半双层阻力积分并对对称性加倍后取倒数，数值自相关积分给出CNFD内在P
eff = 9.09 × 10
?2 cm s
?1，拟合PACF敏感性分析为1.14 × 10
?1 cm s
?1。阻力曲线显示，操作性核心区（z ≤ 0.5 nm）虽然只占坐标区间10.3%，却容纳了81.4%的积分阻力。
**3.4 热力学与动力学阻力的反事实分离**
通过受控剖面替换（保持PMF/扩散之一不变，另一设为平坦或常数），相对于零PMF、常数扩散的平坦参考，PMF-only放大因子AF = 10.54，扩散-only因子AD = 14.89，残余空间耦合因子Γ = 1.134；对应RT ln(A)分别为6.07、6.97、0.33 kJ mol
?1。相同分解应用于比较物时发现不同阻力机制：LAP更受热力学主导（AF = 21.19，AD = 4.85，Γ = 1.453），PLU则PMF因子低于1（AF = 0.845），扩散因子为9.07，Γ = 0.800，表明其有利热力学部分抵消动力学阻力。中心区阻力占比：CNFD 81.4%，LAP 94.0%，PLU仅6.5%，PLU中位阻力位于z = 2.35 nm。
**3.5 PMF收敛与重复变异性**
三个独立ST-eUS副本的全剖面两两PMF RMSD分别为1.260、1.388、1.364 kJ mol
?1，均值1.34 ± 0.07 kJ mol
?1；限定z ≤ 3.5 nm区间为1.51 ± 0.08 kJ mol
?1，表明副本间一致性良好。
**讨论总结**
研究人员从物理机制上解释了CNFD的阻力来源：PMF异质性和扩散异质性各自将阻力放大一个数量级，而空间耦合贡献很小，证实中心阻力并非二者极值共定位的瓶颈。溶剂化分析显示，CNFD从界面移向中心时水合数下降、脂质接触增加，但组分富集分析未发现可重复的脂质特异性关联。方法论上，Hummer位置自相关在约束窗口内只能给出局域平衡D
⊥(z)，不能等同于无约束长时爱因斯坦扩散系数；绝对P
eff值依赖CHARMM/TIP3P力场模型。三种萘醌的比较揭示不同阻力景观：CNFD热力学和动力学贡献相当且耦合弱，LAP热力学主导，PLU具有有利中心PMF但扩散惩罚偏近界面。已发表的LAP和PLU转运数据为模拟提供外部背景，但实验测量涉及水边界层、细胞结构、转运蛋白等因素，不能与内在脂质相P
eff直接互换。在被动扩散为重要途径的前提下，模拟预测CNFD行为更接近LAP而非PLU，该预测需匹配的转运实验检验。
**结论翻译**
这项工作采用成熟的ST-eUS/WHAM、Hummer位置自相关和ISD分析，解析了CNFD在模型脂质双层中阻力的物理起源。受控剖面替换表明，PMF异质性和位置依赖性扩散分别将阻力放大10.54和14.89倍，而残余空间耦合因子仅为1.134。由于PMF最大值与扩散最小值相距2.14 nm，中心阻力并非共定位的热力学-动力学瓶颈。主要内在脂质相P
eff为9.09 × 10
?2 cm s
?1，该值依赖于当前CHARMM/TIP3P模型和数值自相关估计器。相同的分解区分了三种萘醌的阻力机制：CNFD具有显著且相当的热力学与动力学放大，LAP更受热力学主导，PLU则结合有利自由能与偏离膜中心的扩散惩罚。PLU紧凑、平面、脱水且接触丰富的中心构型与其有利中心PMF一致，但现有数据未指定唯一因果相互作用。已发表的LAP和PLU转运测量使条件性实验含义明确：若被动脂质双层渗透是重要共同途径，则预测CNFD的表现更像LAP而非PLU。该预测需匹配的CNFD/LAP/PLU转运实验加以检验。
