脂质膜面积压缩模量的分子分解：POPC双分子层中烃链与界面水合对面积弹性的层级贡献

《Biomolecules》：Molecular Decomposition of the Area Compressibility Modulus of Lipid Membranes

　　

面积压缩模量（area compressibility modulus）是脂质膜力学性质的关键描述符，但决定该弹性的分子相互作用无法仅从总模量中显现。本研究将粗粒化Martini 3力场描述的POPC双分子层分子动力学模拟与维里应力（virial stress）分析相结合，将面积压缩模量分解为分子基团及其相互作用的贡献。研究人员将部分侧向张力（partial lateral tension）评估为膜面积应变的函数，并利用其在平衡态处的导数获得对面积压缩模量的相应贡献。烃链提供主要的正贡献，约占总模量的70%，而水–脂质相互作用贡献约30%，且主要通过水–头基相互作用实现。在双分子层内部，头基与头基–链贡献量级相近但符号相反，几乎相互抵消。直接的单层间相互作用仅占总模量的约2%，表明两单层的弹性响应近乎可加。弹性贡献的层级与平衡部分张力的层级显著不同，表明维持静态应力平衡的相互作用未必是支配膜刚度的相互作用。该分解为脂质双分子层面积弹性提供了微观解释。

研究背景方面，生物膜界定细胞区室、调控分子输运并为众多生化过程提供动态平台，其功能不仅取决于组成与分子组织，还取决于脂质双分子层在形变过程中承受与再分配机械应力的能力。侧向张力（lateral tension，γ）与面积压缩模量（area compressibility modulus，k_A）是描述该力学行为的两个核心量：前者表征给定面积下膜的力学状态，后者衡量双分子层对平面内膨胀或压缩的微分阻力。面积弹性远不止于均匀拉伸，膜弯曲必然在双分子层厚度方向产生相反的面积应变，因而与单层的侧向压缩和扩张内在耦合，相关形变能量代价参与膜介导相互作用、孔形成、脂质畴边界力学以及膜分裂与融合。在连续介质层面这些效应通常由少量弹性模量描述，但在分子层面刚度的来源并不透明：究竟脂质分子的哪些部分或哪些脂质–溶剂相互作用承载了膜面积变化的力学响应并不清楚。分子动力学模拟通过解析单个分子相互作用产生的应力为这种微观描述提供了直接途径。已有应力分析的关键结果表明，无张力双分子层未必局部无应力，力学平衡可源于膜不同区域与分子组分所关联的大幅正负贡献的抵消；然而这些分解刻画的是静态力学状态，本身并不能确定哪些分子子系统主导弹性响应。这一区分直接来自面积压缩模量作为侧向张力对相对面积应变导数的定义：某个分子子系统在平衡态可承载很大的部分张力，但其张力随形变变化很弱；反之，小的平衡贡献可能具有陡峭的应变依赖从而对模量贡献很大。因此膜刚度的分子起源编码于各分张力对应变的响应而非其幅值之中。基于此，研究人员对Martini 3力场描述的POPC双分子层开展了本研究，旨在将应变相关的侧向张力分解并微分，从而把部分面积压缩模量归属到特定结构基团与相互作用，分离链堆积、水合与单层间耦合的作用。研究得出，烃链基团提供面积弹性的主导贡献，水–脂质相互作用提供可观的次级贡献，而直接单层间耦合尽管承载很大的平衡部分张力却对面积模量贡献很弱；通过在同一框架内解析平衡应力与微分刚度，该分析为膜结构与水合如何决定脂质双分子层抵抗面积变化的分子图景提供了依据。其重要意义在于明确了承载静态应力平衡的相互作用与支配微分刚度的相互作用并不一致，部分应力对应变的依赖性蕴含了平衡应力分解所无法提供的力学信息。

关键技术方法方面（不超过250字）：研究采用粗粒化Martini 3力场对含256个POPC脂质（每单层128个）的双分子层进行分子动力学模拟，水含量为每脂质11个标准水珠子（相当于每脂质44个水分子，共11264个水分子），体系经CHARMM-GUI构建，使用GROMACS 2024.5模拟，在?100至150 bar nm的一系列施加侧向张力下各平衡100 ns并进行3 μs生产运行；采用GROMACS-LS基于Irving–Kirkwood–Noll定义计算维里应力张量，并通过从轨迹中移除不属于目标基团的珠子、以单网格单元配置计算应力，实现水、脂质双分子层、上下单层及其相互作用、头基与烃链及其相互作用、水–头基与水–链相互作用的多层级分解；部分侧向张力按含瞬时膜面积平均的关系式计算，部分模量由部分张力–应变关系的二次多项式最小二乘拟合在零应变处的线性系数确定，并以分块重采样的自举法（1000次）估计不确定性。

研究结果方面，第一，平衡部分张力：总侧向张力在参考态为零（?0.006 ± 0.011 k_BT nm^?2），但各分子贡献并不小，力学平衡源于大幅正负部分张力的抵消，且该抵消可在若干体系划分下独立验证；脂质子系统内上下单层分别承载3.203 ± 0.010与3.199 ± 0.011 k_BT nm^?2的正部分张力，单层–单层相互作用贡献为?9.045 ± 0.003 k_BT nm^?2，使双分子层总张力为负；以头基、烃链及头基–链表示时三者分别为11.586 ± 0.010、?3.281 ± 0.011与?10.948 ± 0.006 k_BT nm^?2，求和与双分子层值一致；水单独贡献6.841 ± 0.011 k_BT nm^?2，水–头基与水–链相互作用分别贡献?2.424 ± 0.010与?1.779 ± 0.002 k_BT nm^?2，溶剂相关合计约+2.639 k_BT nm^?2，补偿双分子层值约?2.643 k_BT nm^?2；平衡部分张力幅值分布很不均匀，最大单基团贡献为正的头基张力，头基–链与单层–单层相互作用承载相近量级的负张力，水亦贡献可观正值，该排序与相应弹性导数的层级差异很大。第二，部分面积压缩模量的提取：按定义将各贡献的γ_ij(ε)以二次多项式拟合，k_A,ij取ε = 0处导数即线性项系数；以烃链贡献为例，其部分张力在模拟应变范围内平滑变化，拟合在原点处的导数给出k_A,ij = 50.93 ± 0.19 k_BT nm^?2，其余部分张力以同样方式分析，故所有模量均来自同一定义而非平衡张力幅值。第三，面积弹性的分子基团层级：烃链基团提供主导正项50.93 ± 0.19 k_BT nm^?2，占总模量71.85 ± 0.34 k_BT nm^?2的70.9%，该贡献本身几乎等于双分子层模量50.48 ± 0.26 k_BT nm^?2；其余双分子层内部项大而强烈补偿，头基贡献为?15.01 ± 0.10 k_BT nm^?2，头基–链相互作用贡献+14.56 ± 0.12 k_BT nm^?2，二者之和接近零，与烃链项相加得50.48 k_BT nm^?2；水–脂质相互作用构成第二大正贡献但在两类脂质基团间高度不对称，水–头基项为22.08 ± 0.12 k_BT nm^?2，水–链项略负为?1.42 ± 0.04 k_BT nm^?2，二者直接评估的合计为20.67 ± 0.12 k_BT nm^?2，占总模量的28.8%；水自身直接贡献很小，为0.70 ± 0.16 k_BT nm^?2，与上述合计相加约21.37 k_BT nm^?2，占29.7%；因此两个最大正弹性贡献分别关联烃链与水–头基相互作用，其他若干项主要作为补偿修正。第四，单层分解与单层间耦合：上下单层分别贡献24.62 ± 0.18与24.36 ± 0.21 k_BT nm^?2，仅相差0.26 k_BT nm^?2，其和48.98 k_BT nm^?2已占双分子层模量约97%，剩余1.50 ± 0.06 k_BT nm^?2（约占总模量2%、双分子层模量约3%）归属单层–单层相互作用，三者相加与双分子层值一致，表明两单层对均匀面积形变近乎可加响应，直接单层间耦合仅提供小修正；这与平衡张力分解形成鲜明对比，同一单层–单层相互作用承载?9.045 ± 0.003 k_BT nm^?2的大负平衡张力，对静态应力平衡重要但对微分面积刚度较弱。第五，平衡应力与微分刚度的对比：烃链项具有中等的负平衡张力?3.281 ± 0.011 k_BT nm^?2却提供主导弹性贡献50.93 ± 0.19 k_BT nm^?2；水–头基相互作用平衡张力仅?2.424 ± 0.010 k_BT nm^?2而对k_A贡献达22.08 ± 0.12 k_BT nm^?2；反之，若干最大平衡应力对面积模量贡献很小甚至为负，单层–单层相互作用平衡时承载?9.045 ± 0.003 k_BT nm^?2但对k_A仅贡献1.50 ± 0.06 k_BT nm^?2，水承载6.841 ± 0.011 k_BT nm^?2的大正平衡张力而直接弹性贡献仅0.70 ± 0.16 k_BT nm^?2，最显著的是头基项具有最大的正平衡部分张力11.586 ± 0.010 k_BT nm^?2却对应?15.01 ± 0.10 k_BT nm^?2的负弹性贡献，头基–链相互作用则相反，承载?10.948 ± 0.006 k_BT nm^?2的大负平衡张力而弹性贡献为正的14.56 ± 0.12 k_BT nm^?2；可见ε = 0处γ_ij的量级与符号可与k_A,ij无关，正负部分张力以相反符号参与平衡张力平衡，而正负部分模量分别表示相应部分张力随面积应变增大而增大或减小，正部分模量增大总模量、负贡献减小总模量；对面积弹性而言相关量是参考态处局部斜率dγ_ij/dε而非γ_ij本身；由相关关系式可见，对烃链第二项贡献47.6 k_BT nm^?2显著超过第一项3.3 k_BT nm^?2，而对单层–单层相互作用第二项?7.5 k_BT nm^?2被第一项+9.0 k_BT nm^?2大幅补偿而得小模量，说明两项相对大小因贡献而异，零应变处部分张力本身无法确定部分模量的大小或排序。综上，分解确立了POPC面积压缩响应的四个定量特征：烃链主导正弹性响应；水–脂质相互作用提供可观的次级贡献且集中于水–头基相互作用；两单层近乎可加响应且直接单层间耦合弱；弹性贡献层级与平衡部分张力层级根本不同，分子基团、溶剂/双分子层与单层划分提供了同一力学响应的替代可加表示。

讨论部分总结：本研究通过考察各分子贡献对侧向张力的应变依赖性，扩展了以往对平衡部分张力的分析，区分了承载静态侧向应力与贡献面积压缩模量这两种力学角色，并因若干分子贡献在平衡部分张力与应变导数比较中排序与符号均改变而使该区分尤为清晰。总模量为71.85 ± 0.34 k_BT nm^?2，在300 K下约相当于298 mN m^?1，略高于若干早期流体POPC膜的实验估计（约200–215 mN m^?1），但更接近Drabik等报道的249 ± 35与255 ± 31 mN m^?1，实验值在研究中差异明显，本值略高于中心实验估计但量级相当。烃链贡献50.93 ± 0.19 k_BT nm^?2占总模量70.9%，为主要正项，表明对均匀侧向面积形变的主要阻力关联酰基尾区域本身；这种近乎相等源于其余双分子层内部贡献大体相互补偿，头基单独贡献大且为负（?15.01 ± 0.10 k_BT nm^?2），头基–链相互作用贡献+14.56 ± 0.12 k_BT nm^?2，二者近乎抵消；负头基贡献应解释为参考态处头基部分张力随面积的负斜率而非孤立头基子系统的不稳定，部分模量是可加分解的组分，各自不必为正。水合提供第二大贡献，水–头基22.08 ± 0.12 k_BT nm^?2，水–链仅?1.42 ± 0.04 k_BT nm^?2，合计20.67 ± 0.12 k_BT nm^?2占28.8%，而水自身直接项仅0.70 ± 0.16 k_BT nm^?2；溶剂力学角色主要通过其与膜的相互作用体现，且在本划分中主要定域于头基–水界面，因头基含甘油部分以及胆碱和磷酸珠子，该贡献代表更广的水合界面区域而非仅带电珠子。单层分解揭示类似区分：上下单层贡献之和48.98 k_BT nm^?2已占双分子层模量约97%，直接单层–单层相互作用仅1.50 ± 0.06 k_BT nm^?2，尽管同一相互作用承载?9.045 ± 0.003 k_BT nm^{?2的大平衡部分张力；两单层对均匀面积压缩与拉伸近乎可加，尽管其相互作用对平衡双分子层的内张力平衡至关重要。POPC中较小的单层间贡献不应假定适用于分子组织不同的膜，更强的单层间相互作用、更大的链互穿或链有序性变化可能改变单层间部分张力的应变依赖性从而改变其对面积压缩模量的贡献，链堆积变化也可能影响单层内贡献，因此更强耦合或更有序膜中单层内与单层间相互作用的相对重要性有待进一步研究。约占总模量2%的单层–单层直接小贡献可置于单层–双分子层等效性问题背景下考察，例如将双分子层内单层与气–水界面Langmuir单层比较，弹性响应等效与单层贡献可加性将意味着双分子层面积压缩模量为相应Langmuir单层的两倍，但本研究中双分子层内两单层近乎可加的弹性响应并不意味着双分子层与相应Langmuir单层面积压缩模量呈2:1之比，因为两个体系界面环境与平衡分子构型不同，尤其单层间相互作用即使直接贡献小也可能改变烃链堆积进而改变单层内弹性响应。本研究计算的侧向张力指水合双分子层体系，应与Langmuir单层的表面压区分，后者定义为裸气–水界面与脂质覆盖气–水界面表面张力之差；在与完全水合双分子层比较时，约30 mN m?1的Langmuir单层表面压被认为对应相似的分子堆积条件，取裸气–水界面表面张力约72 mN m?1，则对应脂质覆盖界面张力约42 mN m?1；原则上本研究的维里应力形式也可用于分解脂质覆盖气–水界面的张力，但这些贡献必须针对Langmuir单层体系评估而不能由本双分子层计算直接推断；由本双分子层部分张力可将单层内贡献与水和水–脂质贡献的一半组合，构造出在所含相互作用意义上类似理想化水–真空界面单层张力的量，其在双分子层系综而非孤立单层平衡态中评估，约为18.7 mN m?1，该构造值与脂质覆盖气–水界面估计张力（约42 mN m?1）的差异并不表示不一致，因为界面环境与平衡分子构型不同。对本文考虑的对称双分子层，可将上下单层内贡献的平均值与水和水–脂质贡献的一半组合构造平均有效单层贡献，约为35.17 k_BT nm?2，等于总模量与直接单层间贡献之差的一半，故双分子层总模量约为该有效贡献的两倍；如上所述该关系涉及在双分子层系综中评估的贡献，并不确立Langmuir单层的相应关系。更一般地，平衡部分张力不是弹性刚度的替代指标：头基项提供最大正平衡张力却对k_A为负贡献，头基–链相互作用呈相反符号模式，烃链平衡张力相对中等却主导模量，水承载可观平衡应力而直接对刚度贡献很小；弹性相关量不是参考态部分应力幅值而是该应力在膜受应变时变化的多快，这在解释分子应力分解时很重要，维持力学平衡的相互作用未必是支配形变能量代价的相互作用。该分解还为考察膜组成或环境如何再分配弹性提供自然框架，例如改变酰链堆积（即烃链的空间排列与构象）的扰动可能通过改变链相互作用对面积形变的响应而影响主导的烃链贡献，水合、头基化学或离子条件的变化应主要通过水–头基及头基相关贡献体现，这些分子组织变化可用傅里叶变换红外（FTIR）光谱表征，CH₂伸缩振动提供烃链构象与堆积信息，磷酸与羰基伸缩振动提供水合信息。胆固醇等链有序化添加剂为将该分析扩展到多组分膜提供了自然测试案例，胆固醇可作为独立分子基团处理，其相互作用纳入分解并在各组成下评估部分模量，分解虽保持可加，但贡献的大小与相对层级可能非线性依赖胆固醇浓度，量化这些变化需显式模拟而非从本POPC双分子层外推。面积压缩与拉伸因弯曲在单层厚度方向产生相反面积应变而与膜弯曲内在耦合，但弯曲与均匀面积形变在若干重要方面不同，尤其涉及空间非均匀应变场并可能伴随额外分子重排，包括脂质堆积与水合变化及可能形成的疏水缺陷，因此本研究所得面积压缩模量的分子分解不能直接移用于弯曲弹性，不同分子基团与相互作用的相对贡献在两种形变模式间可能不同，弯曲模量的定量分解需专门针对膜弯曲另行研究。本分析也限于对均匀面积形变的平衡响应，未刻画时间相关的剪切与黏弹性质，单层间相互作用对面积压缩模量的小贡献并不意味着弱单层间摩擦，因为二者刻画不同力学响应，烃链间分子重排与相互作用可能在剪切弛豫与单层相对滑移中起显著作用，量化剪切黏度与单层间摩擦需单独的动力学研究。本研究限于Martini 3力场描述的单组分POPC双分子层，所报道部分面积压缩模量的数值与相对大小不应视为脂质膜的普适性质，不同脂质组成与力场可能导致弹性响应在分子基团及其相互作用间分布不同；此外各贡献依赖所选脂质珠子划分，本分析将甘油珠子归入头基而烃链基团仅含酰基尾珠子，但可加关系确保部分贡献重构相应双分子层或总响应，所提出的分解框架本身不限于POPC或Martini 3，可应用于其他膜体系。研究结论部分翻译如下：通过将应变相关的侧向张力分解为分子基团与相互作用贡献，本研究为POPC膜的面积压缩性提供了微观解释。烃链区域提供主导正贡献，约占总面积压缩模量的70%，而水–脂质相互作用贡献约30%，主要通过水–头基相互作用实现。在双分子层内部，负头基贡献被正头基–链项近乎补偿，使烃链成为双分子层刚度的主要来源。两单层近乎可加响应，直接单层间相互作用仅占总面积压缩模量的约2%。与平衡部分张力的比较进一步表明，维持静态应力平衡的分子相互作用未必是支配微分刚度的相互作用，因此部分应力的应变依赖性蕴含了仅凭平衡应力分解无法获得的力学信息。尽管定量结果特定于本文考虑的POPC双分子层与Martini 3力场，所提出的框架可扩展到其他脂质组成与膜环境，以确定分子堆积、水合与界面相互作用的变化如何再分配膜弹性。}