gmx构建初始模拟体系

  1. 把蛋白A、蛋白B、DNA-C三个组分,按照B-A-C这种夹心的空间结构放置到一起,使用方法:
gmx editconf -f proteinB.gro -o proteinB-box.gro -box 15 15 40 -center 7 7 7
gmx editconf -f proteinA.gro -o proteinA-box.gro -box 15 15 36 -center 7 7 19
gmx editconf -f DNA.gro -o DNA-box.gro -box 15 15 40 -center 7 7 31
cat proteinB-box.gro proteinA-box.gro DNA-box.gro >> complex.gro	# cat之后要修改complex.gro中的内容
#添加溶剂
gmx solvate -cp complex.gro  -p topol.top -o solv.gro
# 可以用vmd观看三者的相对位置是否是符合我们的要求

  1. 构建多离子体系
gmx insert-molecules -ci A_GMX.gro -box 8 8 8 -nmol 20 -o A-box.gro
gmx solvate -cp A-box.gro -p topol.top -o solv.gro
gmx grompp -f ions.mdp -c solv.gro -p topol.top -o ions.tpr -maxwarn 2
gmx genion -s ions.tpr -p topol.top -o solv-CA.gro -pname CA -np 50
gmx grompp -f ions.mdp -c solv-CA.gro  -p topol.top -o ions.tpr -maxwarn 2
gmx genion -s ions.tpr -p topol.top -o solv-CA-MG.gro -pname MG -np 20
gmx grompp -f ions.mdp -c solv-CA-MG.gro   -p topol.top -o ions.tpr -maxwarn 2
gmx genion -s ions.tpr -p topol.top -o solv-CA-MG-Cl.gro -nname CL -nn 100

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值