|
|
各位老师，本人正在研究跨膜受体在双分子膜上的多聚化现象。
盒子大小为72*72*12，共均匀放置了36个蛋白（2种*18个）在膜上杠，每个相距12nm呈正方形分布。能量最小化和平衡阶段都跑完了，进入生产模拟阶段大约每运行几小时就会出现lincs崩溃，从检查点重新运行后又能跑过去，就这样断断续续重跑了二三十次，累计到现在跑了大概3.5μs，检查轨迹也都一切正常。现在很希望能解决体系容易崩溃的问题，因为导师希望能跑多次重复并且将单个模拟的时间延长到体系收敛为止。我掐指一算这样跑到毕业都跑不完 。
我的具体流程是这样的：
1. 将两个蛋白的全原子模型在水中进行严格的能量最小化，并各跑了10ns的全原子模拟，取了最后一帧的构象。
2. 对这两个构象使用membrane-master的memembed功能调整orientation，并在pymol中调整蛋白在xy平面的位置，将四个蛋白（每种两个）摆放成正方形并保存。
3. 使用martinize脚本将四蛋白体系粗粒化，并用insane.py插入膜中，添加水及离子。限制backbone beads的位置并进行能量最小化和eq。
4. 对eq完成的四蛋白体系gmx editconf复制8份，变为3*3个四蛋白体系，即36个蛋白的大体系。
5. 能量最小化，并放开位置限制做了一个比较短的eq。
6. 生产模拟，跑了一段时间后出现lincs error。
很希望大家能指出问题出在哪，网上能搜到的方法基本都试过，三种控压方法都试了，tau-p也一直在调整，甚至还试了调大martini范德华和库仑力的cutoff（虽然可能这样做不对），最后的结果都大同小异，即跑一段时间后崩溃，且稳定性没有太大差异。
拜托各位大佬给提提意见！！
下面是我生产模拟的mdp，分号后的基本都是尝试过的不同的mdp设置。
- integrator = md
- tinit = 0.0
- dt = 0.02
- nsteps = 5000000000
- nstxout = 0
- nstvout = 0
- nstfout = 0
- nstlog = 100000
- nstenergy = 100000
- nstxout-compressed = 100000
- compressed-x-precision = 1000
- nstlist = 10 ;; 20
- ns_type = grid
- pbc = xyz
- rlist = 1.4 ; 1.1 ;;1.35
- coulombtype = Reaction_field
- rcoulomb_switch = 0.0
- rcoulomb = 1.2 ;; 1.1
- epsilon_r = 15
- vdw_type = cutoff
- rvdw_switch = 0.9
- rvdw = 1.2 ;; 1.1
- cutoff-scheme = verlet
- coulomb-modifier = Potential-shift
- vdw-modifier = Potential-shift
- epsilon_rf = 0
- verlet-buffer-tolerance = -1 ;;0.005
- tcoupl = v-rescale
- tc-grps = PROTEIN LIPID SOL_ION
- tau_t = 1.0 1.0 1.0
- ref_t = 323 323 323
- Pcoupl = c-rescale ;;parrinello-rahman
- Pcoupltype = semiisotropic
- tau_p = 6.0 ;;12.0
- compressibility = 3e-4 3e-4
- ref_p = 1.0 1.0
- gen_vel = no
- gen_temp = 323
- gen_seed = -1
- constraints = none
- constraint_algorithm = Lincs
- continuation = no
- lincs_order = 4
- lincs_warnangle = 30
- refcoord_scaling = com
|
|
发表于 Post on 2025-3-17 22:45:33
收藏 Add to favorites
楼主 Author