gromacs小案例

//////
//*********************4.5.3////
/////////////
pdb2gmx -f protein.pdb -o protein.gro -p protein.top -ignh
/////////////////
grompp -v -f minim.mdp -c protein-water.gro -p protein.top -o protein-water.tpr////File input/output error:minim.mdp
/////////////////////////
genion -s protein-water.tpr -o protein-solvated.gro -conc 0.15 -neutral -pname NA+ -nname CL-
///////////////////////
grompp -v -f minim.mdp -c protein-solvated.gro -p protein.top -o protein-EM-solvated.tpr
///////////////////////
mdrun -v -deffnm protein-EM-solvated -c protein-EM-solvated.gro
////////////////////
grompp -v -f pr.mdp -c protein-EM-solvated.gro -p protein.top -o protein-PR.tpr
//////////////////////
mdrun -v -deffnm protein-PR
////////////////////////////////
g_energy -f protein-PR.edr -o thermodynamics-PR.xvg
////////////////////////////
xmgrace -nxy thermodynamics-PR.xvg
////////////////////////////////
grompp -v -f nvt.mdp -c protein-PR.gro -p protein.top -o protein-NVT.tpr
////////////////////////////////
mdrun -v -deffnm protein-NVT
//////////////////////////////////
g_energy -f protein-NVT.edr -o thermodynamics-NVT.xvg
///////////////////////////////////
grompp -v -f npt.mdp -c protein-NVT.gro -p protein.top -o protein-NPT.tpr
////////////////////////////////////
mdrun -v -deffnm protein-NPT
////////////////////////////////
g_energy -f protein-NPT.edr -o thermodynamics-NPT.xvg
////////////////////////////////
g_energy -f protein-NPT.edr -o density-NPT.xvg
////////////////////////////////
grompp -v -f md.mdp -c protein-NPT.gro -p protein.top -o topol.tpr
//////////////////////////////////////
mdrun -v
/////////////////////////////////
gmxdump -s topol.tpr | less
//////////////////////////////

评论/留言