综述：量化小分子与脂质膜的关联：方法、模型与局限性

《Membranes》：Quantifying Small-Molecule Association with Lipid Membranes: Methods, Models, and Limitations

　　

小分子与脂质膜的关联在药物递送、药代动力学、毒性和膜生物物理学中起着核心作用，并且由于大多数药物靶点是膜相关蛋白，因此对于药物药效学也具有根本重要性。然而，由于实验系统的多样性、膜环境的复杂性以及各方法固有的局限性，准确测定溶质-膜关联亲和力仍然具有挑战性。本综述全面概述了当前用于量化小分子与脂质膜关联的常用实验和计算方法。讨论了标准实验技术，包括基于光谱的方法、量热法、电泳测量和表面敏感方法，以及从连续介质模型到全原子分子动力学模拟的成熟计算策略。特别强调了数据分析所需的形式体系，包括分配模型和热力学框架，以及每种方法背后的假设。批判性地审视了有效性界限、不确定性来源以及常见的实验和解释陷阱。通过提供统一和比较视角，本工作建立了研究溶质-膜相互作用的定量结构化框架，指导新研究人员选择适当方法并严格分析实验和计算结果。此外，它使得跨文献报道的亲和力参数能够一致且定量地合理化，支持精选数据集和预测关系的开发，从而为设计新的更有效药物提供参考。

1. 1.
   引言
   药物与生物膜的关联是决定药物生物利用度和药代动力学的关键步骤之一。通常使用药物在水与辛醇之间的分配来估计这一性质，但辛醇无法捕获生物膜的显著特征。脂质双分子层是一种结构化的、极薄的、强各向异性介质，从与水接触的极性头部基团区域到双分子层中心的非极性区域仅在几纳米内变化。由于每个单分子层的厚度与大多数生物活性分子的尺寸相当，它们跨越了溶剂化性质截然不同的膜区域，采取最大化与膜有利相互作用的构象、位置和取向。认识到脂质双分子层是界面相，膜分配本质上不同于分配到本体相。将现有关于药物/膜相互作用的数据统一化，并建立公开可用数据库对于提高从分子结构预测膜亲和力的能力至关重要。
2. 2.
   量化膜亲和力的参数
   量化小分子对脂质膜亲和力的参数有多种，参考参数为溶质摩尔浓度之比定义的分配系数K_mw。其他常见参数包括基于重量摩尔浓度、摩尔分数或质量比的参数，这些参数在不同形式下相互可转换。处理脂质膜为与水不混溶的相导致分配系数形式体系，适用于低浓度小分子。对于某些溶质和方法，关联最好用结合形式体系描述，考虑膜中的结合位点以及结合导致的结合容量下降。平衡结合常数K_b与分配系数的关系由方程给出。此外，将整个脂质囊泡视为结合剂的形式体系用于从结合和解离速率常数获得平衡亲和力。早期工作中还使用了定义为与膜表面接触的水介质厚度的参数K。

2.1 静电相互作用和离子强度对观察到的溶质-膜亲和力的影响

在分配形式体系中，假设溶质不影响膜性质，亲和力与浓度无关。但对于带电溶质，区分静电相互作用贡献和疏水相互作用贡献很重要。观察到的分配系数K_obs与不存在静电效应时的分配系数K₀通过膜表面电势ψ₀和溶质电荷z相关。膜表面电势由膜电荷密度σ决定，而σ又依赖于溶质在膜中的量，形成迭代求解过程。对于小表面电势和主要是单价离子的电解质溶液，可使用简化方程直接计算表面电势。在体内条件下，由于膜浓度高、药物总浓度低以及生理离子强度高，带电溶质的局部浓度通常很低，观察到的亲和力与未带电膜相近，但仍需校正以获得内在亲和力参数。

2.2 相互作用动力学对表观溶质-膜亲和力的影响

膜体积指溶质在平衡时间内可接近的部分。溶质在水相和膜之间的平衡涉及多个步骤，速率限制步骤可能是插入、解吸或易位。对于相对极性溶质，易位是速率限制步骤；对于非常亲脂性溶质，解吸附是限制步骤。使用多层脂质体时，由于可接近脂质分数小且平衡依赖解吸速率，可能导致表观分配系数严重低估，因此推荐使用单层脂质体。

2.3 溶质对膜性质及其自身膜亲和力的影响

除了静电和动力学效应，溶质可能引起膜性质的广泛变化，导致浓度依赖性亲和力。例如多价离子与带电膜相互作用、胆固醇等刚性溶质影响膜动力学和侧向组织、两亲性去污剂样药物稳定膜孔等。这些效应在体内可能不显著，但在体外表征时可能影响亲和力测定。

1. 3.
   获取溶质-膜亲和力的方法
   3.1 实验方法
   表征溶质与脂质膜关联的方法多样，可分为直接定量、间接通过溶质性质变化、间接通过膜或系统性质变化。

3.1.1 直接定量至少一个相中的溶质

物理分离水相和膜富集溶液后直接定量溶质。常用方法包括离心、平衡透析、过滤和尺寸排阻色谱。计算分配系数时需校正膜富集溶液中的水相体积。最准确程序是定量两相中的溶质量，仅测水相溶质可能导致高估。

3.1.2 间接通过溶质性质变化定量

当溶质性质随膜关联改变时，可间接定量。主方程假设溶质在水和膜中的信号可加。有效性依赖于可接近膜分数、信号常数以及在高脂质浓度下无膜扰动。常用方法包括紫外-可见光谱、荧光光谱、荧光相关光谱等。荧光各向异性必须按分数荧光强度加权，否则导致错误估计。最佳报告性质是固定波长荧光强度。

3.1.3 间接通过系统整体和/或膜性质变化定量

等温滴定量热法（ITC）测量溶质从水相转移到膜的热变，可提供亲和力、焓变和动力学信息。简单分配模型或独立结合位点模型可用于分析。ITC还可表征电离变化和动力学。其他方法包括基于膜静电性质变化（如电泳迁移率）、脂质体尺寸变化（如动态光散射、热泳）、脂质相变温度变化（如差示扫描量热法）以及固相吸附膜相互作用（如表面等离子体共振、石英晶体微天平）。这些方法通常需要高局部溶质浓度，可能不反映未扰动膜的亲和力。

3.2 计算方法

3.2.1 分子动力学模拟

分子动力学（MD）模拟广泛用于评估溶质-膜相互作用。局部分配系数从过剩自由能曲线计算。总体分配系数可通过积分自由能曲线获得，需考虑积分限和预指数因子以匹配实验定义。无约束模拟可直接从概率密度计算分配系数，但受限于采样时间。增强采样方法如伞形采样需注意膜变形问题。力场选择对结果准确性影响显著。

3.2.2 隐式溶剂模型

隐式溶剂模型如异介电广义Born（HDGB）和COSMOmic提供了计算自由能曲线的有效途径。HDGB将膜建模为异介电环境，加入静电和非极性贡献，并扩展至动态膜变形。COSMOmic基于量子化学和统计热力学，将双分子层分层计算位置依赖的化学势，可快速预测分配系数，适合高通量筛选。

3.2.3 纳入溶质电离变化的考虑

溶质与膜关联时常伴随电离状态变化。可通过独立模拟不同电离形式的自由能曲线，结合热力学循环获得各物种分配系数和膜内电离常数。另一种方法是恒定pH分子动力学（CpHMD），在模拟中动态更新质子化状态，已用于研究药物透膜机制。正确确定位置依赖的电离对准确预测膜分配和渗透至关重要。

1. 4.
   结论
   本文系统分析了量化小分子与模型脂质膜关联亲和力的多种实验和计算方法，明确了各种参数的形式体系和转换关系。通过统一视角，为建立可靠的溶质-膜相互作用数据库和预测模型奠定了基础，有助于药物设计和膜生物物理研究。