分子动力学里的生物膜模拟说透了就是拿一套力场参数把磷脂分子拆成一个个相互作用的粒子放进一个带有周期性边界的水盒子在计算机里堆出一张双层膜然后看着它逐渐稳定、演化再从中提取结构、动力学和热力学性质。这个系列写到第三十八篇我觉得最该拿来讲透的就是这件事。原因很实际它是药物分子跨膜、膜蛋白构象变化、脂质分相、纳米颗粒与膜互作等几乎所有生物物理课题的地基绕不开也没有捷径可走。如果你是从实验背景转到计算方向或者刚接触分子动力学不久这篇文章能帮你把“生物膜模拟到底怎么从零开始做”这件事讲明白。我会把模型怎么选、体系怎么搭、平衡怎么跑、数据怎么分析、常见问题怎么排查都过一遍全部基于我实际跑过的体系和踩过的坑。放心这不是上课是直接给你一套可以落到命令行的流程。生物膜模拟最大的好处是给你一双能直接“看见”原子级动态的眼睛。实验上你看到的是一个群体的平均行为模拟里你能追踪每一个脂质分子如何摆动、翻转、横向扩散甚至胆固醇分子在两层之间如何重新分布。这篇文章会把这件事的完整链路拉通。1. 生物膜模拟到底在解决什么问题1.1 生物膜不只是“一层油”很多人以为生物膜就是一层简单的油脂隔层这话没毛病但太粗糙了。真实的细胞膜是一个由磷脂、胆固醇、各类膜蛋白和糖脂共同构成的高度动态体系。磷脂分子头部亲水、尾部疏水在水环境中自发排列成双层结构亲水头部朝外接触水相疏水尾部相对聚在膜的中间区域。胆固醇分子嵌在两条尾巴中间起着调节膜流动性和机械稳定性的作用。这种结构的复杂程度在于它的物理性质不是由某一个分子决定的而是由成千上万个分子相互作用的统计结果决定的。膜厚度、单分子面积、各链段的取向有序程度、疏水核心的极性分布这些宏观参数在原子尺度上的波动恰恰是细胞执行功能的基础。比如药物分子能不能穿膜取决于膜内部的自由能势垒膜蛋白能不能正确折叠和构象变化取决于它周围脂质环境的化学组成和压力分布。传统实验手段看生物膜就像在漆黑的夜里看一座城市你能看到路灯连成的道路但看不清每个街区里的人在干什么。荧光显微镜可以标记特定脂质或蛋白但受限于光学衍射极限电子显微镜能给出高分辨静态图像但样品固定和染色过程本身就可能改变膜结构。这时候分子动力学模拟的价值就体现出来了。1.2 为什么非得用分子动力学来做分子动力学模拟的底层逻辑并不复杂你给我一套原子坐标、一个力场、一组温度和压力控制方式我就能用牛顿运动方程去积分得到整个系统下一时刻的坐标和速度。每一步积分里每一个原子都受到其他原子的作用力力的大小和方向由力场中的键伸缩、键角弯曲、二面角旋转、非键相互作用等项决定。这样重复成千上万步你就得到了一条“轨迹”包含体系中每个原子在每个时刻的位置信息。分析这条轨迹就能算出任何你关心的物理量径向分布函数、扩散系数、序参数、体系自由能变化等等。这个方法的独到之处在于当你需要回答“某一个脂质尾巴在某个时刻为什么扭成那个角度”这种问题时只要翻轨迹就能找到答案而且可以提供统计意义。当然分子动力学不是万能的。它的核心限制来自力和时间力场参数决定了模拟体系的准确性而时间步长和计算量限制了你能模拟的时空范围。全原子模拟一个包含十万个原子的膜体系最多也就能跑一两个微秒粗粒化模拟可以把时间尺度拉到几十甚至上百微秒但代价是丢失了原子级细节。所以重要的不是“能不能做”而是“怎么做才是对这个具体问题最优的选择”。1.3 一张稳定膜能回答哪些实际问题我举几个最新文献中常见的研究场景你会发现生物膜模拟的产出非常具体药物通透性阿司匹林、布洛芬等药分子从水相进入膜内部的自由能变化沿膜法线方向算个势均力场剖面就能估计它被膜吸收的倾向和穿膜速率。脂质分相胆固醇与饱和/不饱和磷脂的相互作用会导致膜上形成“筏”状微域。粗粒化模型可以直接观察到这些微域的自发形成。膜曲率和形状变化某些蛋白能局部改变膜的曲率用分子动力学可以模拟蛋白吸附后膜怎么弯曲、怎么出芽。纳米颗粒与膜的相互作用表面带电荷的纳米颗粒怎么结合到膜上是否导致膜穿孔、形成孔洞自愈所需的时间尺度。膜蛋白在特定脂质环境中的构象稳定性比如某些通道蛋白需要特定脂质结合才能保持活性构象。每一个场景的第一步都要先搭建一张稳定、均衡、无应力的双层膜。这听起来基础但恰恰是无数人反复出错的地方。下面我一个个拆开讲。2. 建模思路全原子、联合原子、粗粒化怎么选模型选择的本质是在精度和效率之间做权衡有点像手机地图。全原子模型类似街景地图每条街道、每个店铺招牌都看得清清楚楚但一张图只能覆盖很小一块区域粗粒化模型类似五环图看不到具体店铺但能一眼看完整座城市的路网结构和交通流量。选哪种取决于你要回答的问题需要哪一层分辨率。2.1 全原子模型精度高但成本也高全原子模型把每个重原子以及极性氢都显式地保留下来常见力场有CHARMM36、AMBER Lipid17、OPLS-AA等。以CHARMM36为例它经过长期针对脂质双层参数的校准在膜面积/脂质、氘序参数、核磁共振J偶合常数等实验数据上都拟合得相当好。全原子膜的优点是细节丰富你可以在轨迹里清楚地看到一条特定的脂质尾巴如何弯曲、水分子如何偶发地潜入膜内部形成水线、药分子与膜间的氢键网络如何动态变化。缺点是计算成本很高。一个典型的上百个脂质的双层体系加上水分子和离子原子总数轻松超过五万加上PME静电算法跑1微秒的轨迹在主流GPU上通常需要几十卡天的时间。所以我自己的习惯是如果研究需要精确描述氢键、电荷转移、药分子与配体间精确相互作用就用全原子模型如果只是看大尺度的相行为或统计性质全原子不一定划算。2.2 联合原子模型中间路线的取舍联合原子模型也叫UAUnited-Atom模型把碳上的非极性氢直接合并到碳原子中。这样做的好处是原子数能减少大约三分之一到一半。经典代表是GROMOS系列和Berger力场。Berger力场在历史上有很深厚的应用根基很多早期的脂质双层模拟都出自它。但联合原子模型有个尴尬的地方它缺少极性氢导致氢键网络描述不够精细和现代全原子力场相比精度有一定差距。我自己在近几年的工作中已经很少从零搭UA模型了因为主流pdb2gmx工具链对CHARMM36等全原子力场支持得更好误差也更可控。除非你的体系非常适合GROMOS力场的糖类或蛋白参数否则UA模型更多是历史遗产和教学用途。2.3 Martini粗粒化大尺度问题的利器粗粒化模型里目前最实用、社区生态最好的就是Martini力场尤其是Martini 2和Martini 3。它的核心思想是约化自由度通常把四个重原子映射为一个相互作用珠子四个水分子约化为一颗水珠子。这样系统里粒子数减少一个量级相互作用形式也从复杂的键合势变得简单因此可以稳定使用20到40飞秒的大时间步长算力效率比全原子高一到两个数量级。Martini在生物膜领域的应用非常成熟。你想观察多组分层膜体系里胆固醇的分布、脂质筏结构的形成、蛋白跨膜区的错配效应都可以用Martini在几百纳米的大体系里跑微秒甚至毫秒级的模拟。它还支持反向映射backmapping把粗粒化的最终构象映射回全原子坐标作为精细模拟的起点。Martini最需要警惕的坑是熵效应被低估。珠子之间的摩擦力太小导致动力学过程明显过快脂质横向扩散系数实验值大约是10^-8 cm²/s量级Martini模拟直接给出10^-7甚至更大所以用Martini算绝对扩散常数没有太多参考价值只能用相对趋势或配合速度校正因子。如果你想从粗粒化结果里得到精确的绝对数值建议把粗粒化轨迹反映射回全原子再跑一段做定量分析。2.4 隐式膜模型什么时候该放弃显式水隐式膜模型不显式描述脂质分子和水分子而是用平均势场来替代水相和膜化学环境。经典实现有GROMACS的MMT工具包和高介电连续溶剂模型。它的优点是计算极快适合用于膜蛋白构象采样、结构预测、氨基酸跨膜区判定这类强拓扑搜索任务。代价是丧失了一切脂质动态行为拿不到有序参数模拟不了脂质分相也看不到膜在蛋白作用下的响应变形。还有一个问题是膜厚度和疏水核心的径向函数是预设的实际尺度可能会和真实膜有偏差。所以当你需要“膜本身”作为研究对象而不是“膜作为背景环境”时我一般不推荐隐式膜模型。它属于高效但用途有限的工具适合定位在计算管线里做快速筛选。3. 从零搭一个双层膜体系关键步骤与参数这一章完全按我的实操顺序来。以最常见的全原子POPC体系为例从初始结构一路跑到平衡结束整个流程你拿去可以直接复现。3.1 准备工作工具链和力场文件我固定使用的模拟引擎是GROMACS原因很直接开源免费、粒子数量大时并行效率高、文档丰富、生态好。你还需要一个生成初始结构的工具。因为从头用代码把几十个脂质分子按正确取向摆进周期性盒子是一件没必要自己造轮子的事情所以我强烈建议直接用现成的工具。如果是全原子体系首选CHARMM-GUI的Bilayer Builder网页服务。你只需要填上脂质类型、每层脂质数目、盒子大小、水层厚度、离子浓度它就会生成整套GROMACS模拟所需的拓扑文件、坐标文件和mdp文件。这里需要注意CHARMM-GUI生成的默认mdp参数偏向保守适合平衡阶段不直接适合长时间生产模拟你需要自己根据后续目标调整。如果是Martini粗粒化体系常用的搭膜脚本是Insane.py由荷兰格罗宁根大学的研究组公开。它基于Python脚本命令简单python insane.py -l POPC:100 -l POPC:100 -x 8 -y 8 -z 8 -o system.gro -p topol.top这个命令会生成一个8纳米见方的盒子上下两层各100个POPC中间夹着水层同时自动写好topol.top文件。Insane的优势是快速、透明、可自由调整各层脂质比例非常适合多组分膜体系。3.2 搭建流程造膜、溶剂化、加盐以CHARMM-GUI生成的全原子POPC双层为例标准流程如下在CHARMM-GUI选择Bilayer Builder输入脂质类型POPC设定每层脂质数量通常128或200起步水层每侧厚度默认2.5纳米离子浓度设0.15 mol/L NaCl。提交后下载生成的GROMACS文件包解压后检查结构文件。很多初学者忽略这一步直接跑结果崩溃。至少用VMD或者PyMOL打开system.gro看一眼确认双层膜大致形状没有明显扭曲。对坐标做周期化处理。CHARMM-GUI已经帮你把脂质摆进周期性盒子但你还是要在GROMACS里用gmx pdb2gmx -f step6_input.pdb -o system.gro -ff charmm36 -water tip3重新生成拓扑因为不同力场版本对结构和拓扑的匹配要求很严格用现成的topol.top虽然省事但后续加蛋白或者改体系时会出各种兼容问题。加溶剂和离子。GROMACS的gmx solvate命令会按照你的mdp文件里的溶剂模型自动填充水分子加离子用gmx genion把体系里多余的NaCl浓度调整到0.15 mol/L。记住脂质头基本身带负电体系必须加抗衡离子使总电荷为零。平衡前对整个系统做一次能量最小化检查确认体系里没有显著的原子重叠。这里有个细节我特别强调水层厚度别省。很多新手为了减小计算量把水层压到每侧1.5纳米结果膜两侧的水分子受到周期性镜像的相互作用干扰面积/脂质直接跑偏20%。我建议全原子体系每侧水层至少2.5纳米粗粒化体系至少3.0纳米。3.3 分阶段平衡让膜“安定”下来膜体系搭建完之后不能直接放开分子动力学跑必须分阶段平衡。逻辑很简单初始结构是人为堆出来的里面很可能有局部应力和不合理的原子间距如果一上来就让原子自由运动体系可能在几百步内直接爆炸。通常的做法是先能量最小化EM常用最速下降法加共轭梯度法把结构拉到局部能量最低状态然后做NVT平衡再转NPT平衡。NPT平衡这一步尤其重要因为膜的面积和盒子尺寸必须充分调整到与力场匹配的值否则后续面积/脂质会有系统性偏差。以CHARMM-GUI生成的mdp文件为例平衡阶段它对膜内原子加了几轮位置限制力度逐渐减小。我这边的经验是不要跳过最后一步“放开全部限制”的预平衡让体系在无限制条件下至少跑100纳秒再进入生产模拟。这样得到的膜才是真正“松弛”的。命令很直观gmx grompp -f em.mdp -c system.gro -p topol.top -o em.tpr gmx mdrun -deffnm em gmx grompp -f nvt.mdp -c em.gro -r em.gro -p topol.top -o nvt.tpr gmx mdrun -deffnm nvt gmx grompp -f npt.mdp -c nvt.gro -r nvt.gro -p topol.top -o npt.tpr gmx mdrun -deffnm nptGROMACS里的-r选项用于生成位置限制参考坐标这一步是很多初学容易漏掉的地方。如果忘了写-r nvt.gro平衡时体系可能直接被拉回初始构象白跑好几小时。3.4 关键参数背后的物理依据模拟好不好一半看参数设置。这不是玄学每个参数背后都有物理原因时间步长方面全原子体系里最快的振动来自C-H键伸缩振动周期约为10飞秒量级为了保证积分稳定性时间步长必须远小于振动周期。所以CHARMM36全原子模拟默认使用2飞秒时间步长配合LINCS算法约束氢相关键长。如果用了虚拟位点方法可以把步长提到4飞秒。Martini粗粒化则因为缺少高频率振动项常用20飞秒甚至更大。温度控制方面脂质双层体系与实验温度有强耦合。POPC的主相变温度约为270 K所以在室温以上模拟它处于液态无序相Ld相比较合适DPPC的主相变温度约为314 K在300 K附近模拟DPPC双层它可能处于凝胶相面积/脂质会明显偏小。如果你在认真研究温度和相行为一定要查阅所用脂质的主相变温度别拿All-atom力场在实验无法对应的温度下硬跑。压力控制方面膜体系各向异性明显膜平面方向xy和膜法线方向z的响应不同。因此NPT平衡时一定要用各向异性或半各向异性的压力耦合方式只允许z方向独立变化。如果使用各向同性的压力耦合膜面积会被强行压缩或拉伸产生人为应力。GROMACS里用pcoupltype semiisotropic设置compressibility 4.5e-5 4.5e-5 4.5e-5 0这样的参数就可以。静电处理方面全原子体系用PMEParticle Mesh Ewald处理长程静电截断通常取1.0到1.2纳米。Martini体系里因为珠子间相互作用更平滑短程截断1.1纳米即可。如果你在Martini里也幻想能用PME省事那就错了Martini非键作用的处理方式是为它的参数专门调试过的换掉标准截断方案等于自毁校准。此外脂质双层的项部结构有极性头基和带正负电荷的基团所以离子浓度不可随意调整。常用的生理条件0.15 M NaCl是模拟的默认配置但在计算自由能时离子浓度的微小变化会造成不可忽略的系统误差所以每次对比前都必须保持相同离子环境。4. 数据分析算完之后怎么提取有效信息好不容易跑完模拟轨迹文件躺在那下一步就是怎么把轨迹变成可发表的数字。这部分我按分析对象来梳理。4.1 序参数先看力场对不对脂质碳链的序参数是检验整个体系是否正确平衡和力场是否合理的第一道关卡。实验上可以通过核磁共振测量碳氢键的氘序参数SCD每个碳链的SCD值随链段位置变化形成一个特征剖面模拟结果可以直接和实验值对比。SCD的定义是SCD (3 cos²θ - 1) / 2其中θ是C-D键与膜法线通常为z轴之间的夹角角括号表示对系综平均。GROMACS里直接用gmx order命令就能计算每条链的序参数。实际操作时注意两点一是分析前必须把轨迹的周期性处理干净用gmx trjconv -pbc whole处理二是要选对脂质类型和碳原子索引否则输出的值对不上实验。如果你的模拟体系在长时间跑完后SCD和实验值系统性偏低力场基本没问题大概率是面积/脂质偏大导致链段取向无序如果SCD严重偏离实验值检查一下力场版本、水模型和脂质拓扑是否匹配。4.2 面积/脂质最常用的平衡判据面积/脂质Area Per Lipid是脂质双层模拟里信噪比最高的指标之一。定义也很直观膜平面上的总表面积除以单层脂质分子数。在NPT模拟中xy平面盒子尺寸会随着模拟推进而波动取生产阶段平均值就能得到面积/脂质。计算公式形如A_L (Lx × Ly) / N其中N是单层脂质数目不需要除以总脂质数。原因很简单双层膜有两个单层每个单层独立占据一份表面积总面积除以单层分子数才是每个脂质分子的实际占比。数值参考上纯POPC双层在303 K、0.15 M NaCl条件下CHARMM36力场给出的面积/脂质约为0.64平方纳米。如果你的结果在0.60到0.68之间可以认为体系基本正常如果低于0.55或高于0.75平衡阶段大概率有环节出错了。面积/脂质对体系是否平衡非常敏感所以它是判断“要不要继续跑”的核心指标。我通常会在生产模拟里每10纳秒算一次面积/脂质画一条随时间变化的曲线如果趋势线还在明显爬升说明膜还在收缩或膨胀平衡还没结束。这一条判据比只看体系能量稳定可靠得多。4.3 密度分布与径向分布函数把体系沿z轴方向切成若干厚度的薄片统计每一层内不同组分的原子质量密度就得到密度分布剖面。这个分析能直观地告诉你磷脂头基的极性峰在什么位置、疏水尾巴的密度谷有多深、水分子是否渗透到了膜疏水核心。通过头基峰的位置差还能计算膜的疏水厚度和总厚度。GROMACS里用gmx density命令可以分别对POPC头基原子、尾巴原子和水分子做密度分析。输出画成图后正常的膜剖面应该是头基峰在两侧、尾巴密度在中间平缓分布、水分子基本被阻挡在膜外部只有极少量水分子短暂地潜入疏水核心。如果你看到水分子大量出现在膜中央说明膜太薄或力场参数有问题得回去查水层厚度和体系压力耦合设置。径向分布函数RDF则在分析脂质之间的空间关联时使用。比如看胆固醇与饱和脂质尾巴之间的RDF能判断胆固醇是否更倾向与特定链段堆叠看头基上的离子与磷酸基团之间的RDF可以估计离子确实结合在膜表面的哪一层。实现上一般用gmx rdf命令注意选取合适的原子组和截止距离。4.4 横向扩散系数怎么算才靠谱脂质在膜内的横向运动可以用MSD均方位移方法来描述。MSD随时间增长在二维扩散模型下MSD≈4Dt其中D是横向扩散系数。GROMACS里可以用gmx msd命令对脂质分子的质心做计算但前提是轨迹里已经去除了整体的质心漂移否则算出来的MSD会被人为抬高。实际计算时还有两个坑需要留意。第一个是长时间限制的选取MSD的线性区间通常只出现在中间时间段太短还没进入扩散区域太长会碰到周期性边界导致粒子回到起点。一般做法是对每个脂质分子在轨迹范围内最多允许移动盒子的四分之一到三分之一超出就截断。第二个是统计误差单条脂质链的MSD噪音非常大至少要计算全部脂质的平均结果再做三次重复模拟取标准差。粗粒化Martini体系里算扩散系数时要记住速度被高估的问题。Martini模拟得到的横向扩散系数通常比实验值大约4至5倍。所以如果你想发表绝对扩散常数必须用速度校正因子或者用全原子结果做交叉验证。如果只是比较不同组分膜的流动性相对变化粗粒化模拟的定性结论还是可信的。5. 常见问题与排查技巧实录最后这部分是我翻车记录里最有价值的沉淀全部来自过去几年亲测过的案例。5.1 体系崩溃和原子飞出先别急着重来模拟一开始就跑飞无非几个原因初始结构原子重叠。这是最常见的通常表现为能量最小化不收敛能量呈指数级增长。解决办法是回到CHARMM-GUI或Insane的输入文件检查脂质头基之间有没有明显靠得太近必要时减少每层脂质密度。时间步长过大。如果你把全原子模拟步长拉到4飞秒以上又没有使用虚拟位点LINCS约束算法会报错甚至直接输出NaN坐标。恢复正常的最快方法是把步长降到2飞秒或者按CHARMM-GUI默认参数走。位置限制参考坐标错误。前面说过平衡阶段用-r指定参考结构如果参考坐标搞反了系统会被强行拉向一个不存在的构象表现为体系逐渐变形而不是暴力爆炸。排查方法是在VMD里用时间轴回访确认是否还在膜状结构附近。力场与拓扑不匹配。比如把CHARMM36的topol配了个GROMOS的坐标文件那真是神仙也救不了只能在报错信息里一行行核对。记住一条原则跑崩了的模拟永远不要指望修修补补能救回来。把体系从头搭一遍通常比痛苦排查消耗的时间少得多。5.2 膜不对称和脂质翻转的坑构建双层时上下两层脂质的种类和数量如果不完全相同就叫不对称膜。真实细胞膜本来就是不对称的但模拟时要非常小心两层之间的化学势差异会导致脂质分子跨层翻转flip-flop这种翻转的天然时间尺度在几小时到几天级别天然膜里的频率极低。全原子模拟里翻转罕见但粗粒化模拟里翻转频率会显著升高。如果你在Martini模拟里发现本来应该停留在上层的高浓度胆固醇慢慢跑到了下层这可能就是翻转频率过高造成的伪像不一定是真实的物理过程。解决办法包括在初始结构里加上两层之间的组分约束或者在分析中标记每个脂质分子初始所在层按原层归属做统计。不对称膜还带来一个静电荷问题。不同脂质头基可能带有不同静电荷两层总电荷若不相等体系整体净电荷不为零会引入长程静电相互作用误差。我建议在任何模拟前用gmx grompp生成的拓扑文件确认总电荷再用gmx genion补充或去除离子。5.3 盒子尺寸与周期性边界导致的人为效应周期性边界条件下每个盒子里的粒子会和它在邻近镜像盒子里的镜像发生相互作用。最直接的人为效应就是水层太薄造成的跨盒相互作用设想你的膜每侧水层只有1.5纳米那么对面盒子里的水分子头和你的头基本团距离太近互相影响膜表面水合结构。第二个常见问题是脂质头基或长尾巴跨过盒子边界。分析轨迹时如果不做PBC解缠绕处理脂质会被不连贯地“剪断”MSD计算整个乱套。解决方法是分析前用gmx trjconv -pbc whole重新处理轨迹把穿过盒子的分子映回连续空间。第三个小坑是盒子尺寸不同时膜面积发生改变。同一个膜体系在xy平面上如果盒子太小膜的横向波动和脂质扩散都被限制这个体系就和无限大膜有本质区别。一般建议盒子边长至少是膜直径的2至3倍对128个脂质每层的体系边长约为8到10纳米还是够用的如果研究的是纳米颗粒或蛋白诱导膜变形盒子要更大。5.4 模拟时间跑多长才有统计意义很多新手跑完100纳秒就开始处理数据我却建议先问自己一个问题这个体系真正需要观测的时间尺度是多少脂质分子在膜平面内的横向扩散全原子模拟中每微秒大约只移动几个纳米一些膜曲率变化或蛋白结合事件需要几微秒到几十微秒才会出现明显信号。纯膜体系的平衡本身往往需要几百纳秒到微秒级别过了这个时间点才算“进入稳态”。因此我的建议是全原子膜平衡阶段至少跑200纳秒生产模拟至少500纳秒起步最好跑到1微秒量级粗粒化体系可以大胆跑数十微秒但要注意时间尺度膨胀问题不要把模拟时间直接解读为物理时间。统计上重复三次独立模拟几乎是所有正规期刊的底线要求因为单次模拟很难给出可信的误差估计。判断是否已经达到统计收敛我固定用三件套看面积/脂质对时间曲线是否平稳、看序参数早中晚三个时间段的平均值差是否在噪声范围内、看体系总能量漂移是否小于1%。这三关都过了心里才有底说现在这个数据可以拿去用了。最后分享一点我自己的体会。生物膜模拟这个领域真正拉开差距的地方往往不是高端分析方法而是能不能把体系搭稳、把平衡跑透、把关键参数选对。我自己搭膜这几年最深的感触是造膜不难难在别小瞧平衡和判据。很多人拿默认参数跑了一天结果发现面积/脂质飘得离谱回头查才发现力场版本和水模型对不上。所以无论从CHARMM-GUI还是Insane拿到初始结构一定要自己把拓扑、残基名、力场版本逐一核对再开始跑。如果你正开始接触生物膜模拟我的建议是先拿纯POPC双层练手把从搭建到平衡再到分析的全流程走通再逐步加胆固醇、加蛋白、加多组分复杂体系。路径跑熟了后面遇到复杂问题才有底气一步一步排查。