Curvature feedback between anisotropic membrane proteins and flexible membranes drives protein sorting and shape remodelling

技术摘要：各向异性膜蛋白与柔性膜之间的曲率反馈

问题陈述
虽然诱导曲率的蛋白质已知可以塑造细胞膜，但控制蛋白质组织与膜形态之间反馈循环的物理原理仍未得到完全理解。具体而言，需要一个统一的理论框架来解释各向异性蛋白质（如弯曲棒状 BAR 结构域蛋白）如何感知曲率、如何在异质表面进行自我筛选，以及如何主动重塑膜形状。关键的未解决问题包括：总蛋白质浓度如何影响在复杂几何表面上的空间分布；向列有序（nematic ordering）和拓扑缺陷（TDs）如何共同调节定位；以及在何种条件下，集体结合会驱动诸如颈部稳定（neck stabilization）或管状化（tubulation）等形态转变。

研究方法
作者开发了一个介观理论框架，对部分覆盖着弯曲各向异性棒状蛋白质的闭合柔性膜进行建模，这些蛋白质表现出面内向列有序性。该系统被视为一个曲率耦合的向列型表面。

• 自由能泛函： 通过自洽最小化总自由能 ($F_{tot}$) 来确定蛋白质分布、取向场和膜形状。该泛函包含：
  • 各向同性弯曲 ($f_b$)： 由脂双层中的 Helfrich 自发曲率模型描述。
  • 取向凝聚 ($f_c$)： 通过 Landau-de Gennes 项进行建模，以描述在相变温度以下的向列有序。
  • 固有弹性 ($f_i$)： 惩罚向列序参数的空间变化（直接相互作用）。
  • 偏离弯曲 ($f_e$)： 一个关键项，将蛋白质的固有曲率 ($C_p$) 与膜的局部主曲率 ($C_1, C_2$) 相耦合。该项捕捉了曲率失配的能量代价。
• 约束条件： 总相对表面蛋白质浓度 ($\phi_{tot}$) 是固定的，而局部分布 ($\phi$)、向列序矢量场 ($\mathbf{n}$) 和膜形状是变量。
• 模拟方法： 数值模拟利用在离散轴对称网格（$101 \times 101$ 个点）上的蒙特卡洛方法进行。膜形状由通过傅里叶级数展开的轮廓曲线定义。最小化循环在确定平衡向列配置、计算局部浓度分布以及更新膜形状之间进行迭代，直到收敛。

主要结果

1. 曲率筛选机制：

   • 蛋白质定位受定量曲率匹配而非仅仅是曲率符号的支配。蛋白质聚集在局部主曲率与它们的固有曲率 ($C_p$) 最匹配的地方。
   • 拓扑缺陷 (TDs)： 向列场中的拓扑缺陷充当了能量障碍。即使局部曲率在其他方面是有利的，蛋白质也会在缺陷核心处减少，因为缺陷的能量代价超过了曲率匹配带来的收益。
   • 反直觉的筛选： 具有弱负固有曲率 ($C_p = -2$) 的蛋白质优先占据弱正曲率球形区域 ($C \approx 1.05$)，而非高度负曲率的颈部区域 ($C \approx -24$)，因为前者的量级失配较小。

2. 膜重塑与颈部稳定：

   • 弯曲蛋白质可以在相对较低的浓度下主动稳定薄膜颈部。
   • 颈部半径控制： 稳定颈部的半径随总蛋白质浓度 ($\phi_{tot}$) 的增加和固有曲率量级 ($|C_p|$) 的增加而减小。
   • 竞争关系： 这种行为反映了各向同性膜弯曲能（惩罚高曲率）与蛋白质诱导的偏离弯曲能（倾向于曲率匹配）之间的竞争。
   • 解析推导： 作者推导出了稳定颈部所需的临界浓度 ($\phi_c$) 的解析表达式（式 19），该表达式与数值模拟结果一致。

3. 形态转变：

   • 管状突起： 在中等浓度下，蛋白质会聚集到与其固有曲率半径相匹配的窄管状突起中，而剩余的膜则向其偏好的各向同性弯曲形状弛豫。
   • 赤道环： 在较低浓度下，系统转变为具有赤道环的形状。这种几何结构代表了一种能量折衷：它减少了膜的整体各向同性弯曲能（该能量不利于形成长管），同时仍能在曲率失配有利的区域容纳弯曲蛋白质。
   • 取向依赖性： 正曲率蛋白质通过垂直于颈轴线方向进行取向来稳定颈部，而负曲率蛋白质则平行于轴线取向。

意义与主张
本文声称提供了一个关于各向异性膜及附着的弯曲蛋白质系统中曲率感知与曲率生成的统一物理描述。研究阐明了膜重塑是一个由以下因素相互作用驱动的集体现象：

• 定量曲率匹配。
• 膜弹性。
• 向列取向有序。
• 非均匀蛋白质分布。
• 全局几何约束。

作者断言，这些见解推进了对软凝聚态物质中曲率介导组织的根本理解，并为由各向异性夹杂物驱动的颈部稳定、膜管状化和赤道环形成等生物过程提供了理论基础。该框架不仅适用于生物膜，也适用于被弯曲向列分子覆盖的液晶壳层。