GROMACS version: 2020.1-1
GROMACS modification: No
Here post your question
Dear GROMACS experts
Thanks for the tutorial and question&answers in the forum, I succeed to simulate peptide in DPPC membrane
Now I tried to run 2 peptides in DPPC membrane.
However, when I tried to use ‘grompp’ commend,
; gmx grompp -f minim_inflategro.mdp -c system_inflated.gro -p topol.top -r system_inflated.gro -o system_inflated_em.tpr
an error message appeared.
;
ERROR 1 [file strong_posre.itp, line 212]:
Atom index (208) in position_restraints out of bounds (1-207).
This probably means that you have inserted topology section
“position_restraints”
in a part belonging to a different molecule than you intended to.
In that case move the “position_restraints” section to the right molecule.
According to the manual, this error appeared when placing position restraint files for multiple molecules out of order.
https://manual.gromacs.org/current/user-guide/run-time-errors.html
So I modified topol.top files from
; Include forcefield parameters
#include “gromos53a6_lipid.ff/forcefield.itp”
; Include chain topologies
#include “topol_Protein_chain_A.itp”
#include “topol_Protein_chain_B.itp”
; Strong position restraints for InflateGRO
#ifdef STRONG_POSRES
#include “strong_posre.itp”
#endif
; Include DPPC chain topology
#include “dppc.itp”
; Include water topology
#include “gromos53a6_lipid.ff/spc.itp”
#ifdef POSRES_WATER
; Position restraint for each water oxygen
[ position_restraints ]
; i funct fcx fcy fcz
1 1 1000 1000 1000
#endif
; Include topology for ions
#include “gromos53a6.ff/ions.itp”
to
; Include forcefield parameters
#include “gromos53a6_lipid.ff/forcefield.itp”
; Include chain topologies
#include “topol_Protein_chain_A.itp”
#ifdef POSRES
#include “posre_Protein_chain_A.itp”
#endif
#include “topol_Protein_chain_B.itp”
#ifdef POSRES
#include “posre_Protein_chain_B.itp”
#endif
; Strong position restraints for InflateGRO
#ifdef STRONG_POSRES
#include “strong_posre.itp”
#endif
; Include DPPC chain topology
#include “dppc.itp”
; Include water topology
#include “gromos53a6_lipid.ff/spc.itp”
#ifdef POSRES_WATER
; Position restraint for each water oxygen
[ position_restraints ]
; i funct fcx fcy fcz
1 1 1000 1000 1000
#endif
; Include topology for ions
#include “gromos53a6.ff/ions.itp”
But same error message still is appearing, and I cannot find what I’m missing.
Could you spare me some advice?
Thanks for spare time and attention.