《Protein Science》：Understanding cell-penetrating peptide mechanisms using computational electrophysiology simulations
编辑推荐：
本研究聚焦于细胞穿膜肽（CPPs）的跨膜机制这一未完全阐明的关键科学问题，研究人员引入计算电生理学（CompEL）这一创新计算工具，通过离子失衡诱导膜应力并促使孔道形成，结合双脂双层分子动力学模拟，系统研究了Arg9、MAP、TP10、TP2四种代表性CPPs的跨膜行为。研究发现CPPs可利用预先形成的孔道进行跨膜转运，而疏水性非CPP肽（Leu9）则主要参与孔道稳定，在多个肽共存的体系中还观察到了肽间的协同作用。该研究不仅为理解CPPs的跨膜分子机制提供了新视角，也为下一代药物递送系统的理性设计提供了有价值的见解。
细胞就像一座戒备森严的城堡，其外围的细胞膜构成了坚固的城墙，严格控制着物质的进出。然而，有一类特殊的“信使”——细胞穿膜肽（Cell-penetrating peptides, CPPs），却能够在不破坏城墙结构、不引起细胞毒性的情况下，自由穿梭于城墙内外。它们甚至能携带药物、DNA、纳米颗粒等“重要货物”进入细胞内部，因此在药物递送、基因治疗等领域展现出巨大的应用潜力。尽管CPPs的发现和应用已有时日，但一个根本性的科学谜题仍未完全解开：它们究竟是如何穿越疏水的细胞膜“城墙”的？其内在的分子机制尚不明确，这严重制约了我们对这类分子的理性设计和优化。
为了回答这个关键问题，一篇发表在《Protein Science》上的研究论文，独辟蹊径地引入了一种名为“计算电生理学”（Computational Electrophysiology, CompEL）的分子动力学模拟新策略，深入探究了CPPs与细胞膜的相互作用机制。传统分子动力学模拟虽然能提供原子级细节，但观察完整的CPPs跨膜过程（通常需要秒到分钟）计算成本过高。而CompEL方法通过构建一个双脂双层系统，并在其内外水室人为制造离子浓度差（?Q），从而在膜两侧诱导产生跨膜电位。这种膜应力会促使脂双层自发形成水孔，为CPPs的跨膜提供了“通道”。研究人员利用这种方法，模拟了四种典型的CPPs（带正电的Arg9、两亲性的MAP和TP10、疏水性的TP2）以及一个作为阴性对照的疏水非穿膜肽（nonaleucine, Leu9）的行为。
研究采用的主要关键技术方法包括：1. 计算电生理学模拟：通过建立双脂双层POPC膜系统，设定不同离子失衡值（?Q 0, 8, 12, 16, 24）以诱导跨膜电位和孔道形成，模拟时长为250 ns，并进行了三次重复。2. 分子动力学模拟：使用GROMACS软件和CHARMM36m力场，在350 K温度下运行模拟。肽段结构通过ColabFold（基于AlphaFold）建模。3. 多肽系统模拟：在选定?Q 16条件下，构建包含8个肽的体系（肽脂比1:32），模拟时间延长至500 ns，以研究肽浓度和协同效应的影响。4. 数据分析：利用MDAnalysis、PyLipID等工具分析肽的跨膜行为、脂质占有率、二级结构、氢键、孔道半径及脂质有序参数等。
研究结果
?Q 16条件下单肽系统的跨膜行为
研究人员首先在不同?Q下测试了单肽行为，发现?Q 16时观察到的行为最为多样，便于区分不同肽。他们将肽-膜相互作用分为三类：吸附于内膜、插入孔道但不完全穿越、以及完全跨膜易位。结果表明：
- •
TP10 显示出最高的跨膜易位倾向。
- •
Arg9 和 MAP 主要以吸附为主，但也具备跨膜能力。
- •
TP2 和阴性对照 Leu9 则主要表现为插入并稳定孔道，而无法完成跨膜。
- •
脂质占有率分析显示，带正电的Arg9偏好与脂质极性头部（如磷酸基团）相互作用，而疏水的Leu9则主要与脂质疏水尾部作用。TP2中的色氨酸（W12）残基对稳定其在膜疏水核心至关重要。
- •
二级结构分析表明，Arg9和MAP在模拟中基本无规则，TP10和TP2结构较稳定，而Leu9在插入膜后显著转变为α-螺旋构象。
?Q 16条件下八肽系统的协同效应
考虑到CPPs在生理条件下可能以多聚体形式作用，研究进一步模拟了包含8个肽的体系。结果发现：
- •
肽浓度增加影响了跨膜行为：Arg9和MAP的跨膜事件有所增加；TP10的跨膜反而减少，倾向于形成聚集体稳定在膜界面；TP2在多个肽共存时展现出跨膜能力，表明其存在协同效应；Leu9依然无法跨膜，但能高效插入并稳定孔道。
- •
氢键分析揭示了不同的作用模式：Arg9主要形成分子内氢键，以单体形式通过脂质介导的孔道易位；而TP10、TP2和Leu9则表现出更高的分子间氢键比例，表明其跨膜或插入过程存在肽-肽协同作用。
- •
孔道性质发生变化：在单肽体系中，孔道主要由脂质极性头部稳定，呈现“环形孔”特征；而在八肽体系中，TP10、TP2和Leu9等肽在稳定孔道中起主导作用，脂质参与减少，孔道呈现“桶状孔”特征。
- •
观察到了脂质翻转现象，且仅发生在有肽的体系中，表明CPPs的插入或易位过程可以扰动脂双层结构，促进脂质运动。
跨膜与插入的分子机制
通过对模拟轨迹的细致分析，研究描绘了典型的跨膜（以Arg9为例）和插入（以Leu9为例）过程：
- •
跨膜过程：1. 孔道起始：跨膜电位诱导水缺陷形成。2. 孔道构建与成熟：脂质头部迁移并稳定水孔。3. 肽易位：CPPs被吸引至孔道并穿越。4. 孔道解体：肽完成跨膜后，孔道关闭。
- •
插入过程：疏水肽如Leu9倾向于在孔道形成后插入膜疏水核心，并通过肽间相互作用稳定孔道，阻止其关闭，但自身无法完成跨膜。
研究结论与意义
本研究系统论证了CompEL作为一种强大、快速且可重复的计算工具，在CPPs机理研究中的有效性和独特性。主要结论包括：1. CompEL诱导的跨膜电位可有效促进脂双层孔道形成，CPPs能够利用这些预先存在的孔道进行跨膜转运，而非自身创造孔道。2. 不同理化性质的CPPs展现出迥异的跨膜行为与分子机制：带正电的Arg9主要通过静电作用与脂质头部结合，以单体形式易位；两亲性肽如MAP和TP10具备与脂质头部和尾部相互作用的能力，其行为受浓度和协同效应影响；疏水性肽TP2在达到一定浓度时可借助协同作用实现跨膜；而高度疏水的Leu9则主要起稳定孔道作用，缺乏跨膜能力。3. 肽浓度显著影响其跨膜效率和作用模式，揭示了从单体易位到协同聚集的多种机制“图谱”。
该研究的发现具有多重重要意义。在基础科学层面，它以前所未有的原子细节揭示了CPPs与生物膜相互作用的动态过程，特别是明确了跨膜电位、预成孔道、肽的理化性质及其协同效应在CPPs内化过程中的关键作用，加深了对这一复杂生物物理过程的理解。在方法论上，本研究成功将原本主要用于研究膜蛋白（如离子通道）的CompEL技术拓展至CPPs研究领域，证明了其作为一种增强采样手段，在研究膜活性肽跨膜这类稀有事件中的强大能力，为计算生物学工具包增添了新利器。最终，在应用转化层面，这些分子机制的阐明为基于CPPs的下一代药物递送系统的“理性设计”提供了直接指导。例如，通过合理设计肽的序列（调整电荷、疏水性、两亲性），可以精准调控其与膜的相互作用模式、跨膜效率以及对膜完整性的影响，从而开发出更高效、更安全、更具靶向性的治疗性递送载体，助力癌症、神经退行性疾病、免疫疗法等多种重大疾病的治疗。
