《PLOS Computational Biology》：Clustering of SARS-CoV-2 membrane proteins in lipid bilayer membranes
编辑推荐：
为解决SARS-CoV-2病毒组装初期膜（M）蛋白如何自发聚集并启动病毒颗粒形成的机制难题，研究人员结合全原子分子动力学模拟、连续介质模型与原子力显微镜实验，开展了一项多尺度研究。结果显示，M蛋白间的直接相互作用是驱动其聚集的主导力量，并定量估算出其有效相互作用能和有效寡聚化能。该研究为理解冠状病毒组装和出芽的物理机制提供了定量框架，并指出了抑制病毒复制的潜在新靶点。
想象一下，在微观世界中，SARS-CoV-2病毒如同一个精密的“装配车间”，它在感染细胞内必须将各种“零件”（病毒结构蛋白）精准地组装成一个完整的、具有感染性的病毒颗粒（virion），然后才能“出厂”去感染下一个细胞。这个被称为“组装和出芽”的过程，是病毒生命周期中最关键的步骤之一。对于冠状病毒而言，这个装配线设在细胞内一个名为内质网-高尔基体中间区室（Endoplasmic Reticulum–Golgi Intermediate Compartment, ERGIC）的特殊膜结构上。在病毒的四员“大将”——刺突（S）、膜（M）、包膜（E）和核衣壳（N）蛋白中，M蛋白是病毒颗粒中含量最丰富的结构成分，被认为是组装过程的“脚手架”，负责定义病毒形状并引导其他蛋白。人们普遍认为，M蛋白能相互聚集形成簇，从而诱导膜弯曲，启动病毒颗粒的形成。然而，在这个聚集过程中，M蛋白之间的直接相互作用和由膜变形（如膜变薄）介导的间接相互作用，各自扮演了怎样的角色、孰轻孰重，一直是一个悬而未解的谜题。搞清这个基本物理机制，不仅有助于我们更深入地理解病毒的生命活动，还可能为开发针对病毒组装环节的新型抗病毒策略打开一扇窗。
为了解决这一问题，一篇发表在《PLOS Computational Biology》上的研究，巧妙地采用了多尺度研究方法，将计算机模拟、理论建模和实验观测融为一体。研究人员首先通过全原子分子动力学模拟，在原子层面观察单个M蛋白对周围脂双层膜的影响，精确量化了其引起的膜变薄效应。接着，他们利用原子力显微镜，在接近生理条件的平面支撑脂膜上，直接观测了不同密度下M蛋白的聚集行为。最后，他们构建了一个连续介质模型（Cahn-Hilliard框架），从理论上描述M蛋白密度在膜平面上的演化，并通过线性稳定性分析，找到了驱动蛋白聚集的关键参数。通过将模型预测与实验观测数据直接比较，研究人员得以定量分解出M蛋白聚集过程中的不同能量贡献。
本研究应用了几个关键技术方法：1. 全原子分子动力学模拟：使用CHARMM36m力场和GROMACS软件，模拟M蛋白“短构象”在模拟ERGIC组分的多组分脂双层中的行为，持续2微秒，以分析其诱导的膜厚度变化。2. 原子力显微镜成像与分析：在2.25微米×2.25微米的平面支撑脂双层上，对不同蛋白质量/脂质质量比的样本进行成像，通过图像处理和聚类分析，确定蛋白聚集的临界密度和簇间特征距离。3. 连续介质建模与线性稳定性分析：采用Cahn-Hilliard模型描述M蛋白密度演化，通过微扰分析和色散关系计算，确定蛋白聚集发生的临界条件（临界有效相互作用能和密度），并将理论预测与AFM实验结果进行量化对比。
M蛋白诱导的膜变薄轮廓与线张力
通过2微秒的全原子分子动力学模拟，研究人员量化了M蛋白“短构象”对周围ERGIC样脂膜的变薄效应。模拟显示，蛋白附近膜厚度减少了约0.5纳米，影响范围可达距蛋白中心12纳米处。基于此厚度轮廓和已知的膜弹性常数，计算得出由膜变薄引起的线张力约为0.10 kBT/纳米 ± 0.04 kBT/纳米。这表明膜变形确实可能作为一种膜介导的吸引力，促进邻近M蛋白的靠近。
M蛋白组装依赖于初始密度分数
原子力显微镜实验直观地展示了M蛋白在平面膜上的聚集行为。研究发现，存在一个临界蛋白密度。当蛋白面积覆盖率（密度分数）较低（约0.0016）时，M蛋白以单个体或小寡聚体形式分散存在；当密度提高到约0.161时，蛋白形成各向同性分布的、具有特征间距的紧密簇。这证实了M蛋白在无其他病毒组分的情况下，其自身相互作用足以驱动聚集，且聚集行为强烈依赖于局部蛋白浓度。
有效相互作用能定义聚集的开始
为了定量理解聚集条件，研究建立了连续介质模型并进行线性稳定性分析。分析表明，对于给定的初始蛋白密度分数，存在一个临界有效相互作用能，高于此值才会发生聚集。这个有效相互作用能包含了单个M蛋白所经历的所有最近邻相互作用（包括直接和膜介导的）。理论预测，聚集簇之间的平均距离与有效相互作用能和初始密度相关。通过数值模拟验证，在较低密度下，线性稳定性分析的预测在超越线性区后仍然成立，初始的高密度区域能存活并最终形成稳定的簇。
从AFM图像估算有效相互作用能
将理论模型与最高密度下的AFM实验数据结合，研究人员进行了量化估算。通过测量AFM图像中蛋白簇之间的平均最近邻距离（77.2 ± 1.2 纳米）和蛋白面积覆盖率（0.161 ± 0.005），利用模型反推出M蛋白的有效相互作用能范围为[5.5 kBT, 9.6 kBT]。进一步，通过减去膜变薄线张力的贡献（估算为[0.4 kBT, 1.5 kBT]），得到了主导聚集过程的有效寡聚化能（即直接M-M相互作用能）范围为[4.7 kBT, 8.9 kBT]。
高密度区域在非线性转变中存续
数值模拟表明，在线性区形成的高密度区域（极值点）在向非线性区转变后能够存续，并最终发展为清晰的蛋白簇。在较低的初始蛋白密度分数下，系统会经历一个扩散限制的生长期，随后是缓慢的奥斯特瓦尔德熟化。模拟得到的簇间平均距离与线性稳定性分析的理论预测相符，这支持了使用AFM观测到的稳定簇结构来反推相互作用参数的合理性。然而，在更高的初始密度下，簇会发生合并和粗化，行为更为复杂。
在讨论与结论部分，本研究强调了其发现的多重意义。首先，研究定量证明了在SARS-CoV-2病毒组装的早期阶段，M蛋白间的直接相互作用是驱动其自发聚集的主导力量，而膜变薄所介导的吸引力贡献相对较小。这为解决长期以来的争议提供了明确证据。其次，研究建立了一个多尺度定量框架，将原子尺度的模拟、微米尺度的实验和介观尺度的理论模型无缝衔接，能够精确估算出关键的相互作用能量参数。这些能量值与之前其他跨膜蛋白的研究结果具有可比性，且M蛋白的高阶寡聚化现象也得到了先前多项研究的支持。
研究估算出驱动组装样聚集所需的M蛋白密度分数范围为[0.118, 0.304]。根据对受感染细胞内病毒出芽区室的粗略估算，其M蛋白面积覆盖率的下限约为0.17，正好落在此预测范围内，这从生理学角度支持了本研究发现的可靠性。此外，研究指出，在真实的细胞环境中，膜曲率诱导以及其他病毒结构蛋白（如E蛋白、RNA-N复合物）的存在，可能会进一步降低聚集所需的临界密度和相互作用能阈值，从而在体内更易启动组装。
本研究聚焦于M蛋白的“短构象”，这也是AFM实验中唯一观测到的构象。一个合理的组装情景可能是：组装始于短构象M蛋白的聚集形成簇，随后在寡聚化能或其他蛋白/脂质相互作用驱动下，部分M蛋白发生构象转变，变为与高曲率区域相关的“长构象”。剩余的短构象M蛋白可能定位在出芽颈区，利用其膜变薄效应促进膜切割和病毒颗粒释放。
最后，该研究具有潜在的治疗意义。研究所确定的临界密度和临界有效相互作用能，为干扰病毒组装提供了清晰的量化靶点。例如，通过小分子药物削弱M-M结合能，使其低于临界阈值，或者降低ERGIC膜上的局部M蛋白密度，都有可能有效抑制病毒颗粒的形成，从而阻断SARS-CoV-2的复制。这为开发针对病毒组装环节的新型抗病毒策略提供了坚实的理论基础和新的思路。
