[gmx-users] Membrane protein tutorial

Pawan Kumar pawan.chinari at gmail.com
Thu Apr 9 06:38:37 CEST 2009


Hello sir,

Thanks for such a explanatory tutorial.
I am stuck with one error now.
I have protein, popc, sol and cl- in my system.
I have merged the sol and cl- using make_ndx.
But when I run grompp I get the error like this "popc not defined in the
index file".
The grompp command : grompp -f pr.mdp -c box_em.pdb -p topol.top -o
box_pr.tpr
Any suggestions please.
*
pr.mdp file :*
title               =  protein in popc bilayer
cpp                 =  /usr/bin/cpp
define              =  -DPOSRES -DPOSRES_LIPID
constraints         =  all-bonds
constraint-algorithm=  Lincs
integrator          =  md
dt                  =  0.002
nsteps              =  5000
nstcomm             =  1
nstxout             =  50
nstvout             =  1000
nstfout             =  0
nstlog              =  10
nstenergy           =  10
nstlist             =  10
ns_type             =  grid
rlist               =  1
coulombtype         =  PME
rcoulomb            =  1
vdw-type            =  Cut-off
rvdw                =  1
; Berendsen temperature coupling is on in two groups
tcoupl              = berendsen
tc_grps             = Protein POPC SOL_CL-
tau_t               = 0.1 0.1 0.1
ref_t               = 300 300 300
; Energy monitoring
energygrps        =  Protein POPC SOL_CL-
; Pressure coupling is on
;Pcoupl              =  berendsen
tau_p               =  2.0 2.0
compressibility     =  4.5e-5 4.5e-5
ref_p               =  1.0 1.0
Pcoupl_type         =  semiisotropic
; Generate velocites is on at 300 K.
gen_vel             =  yes
gen_temp            =  300.0
gen_seed            =  173529

Thanking you,
Pawan
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20090409/6c972da3/attachment.html>


More information about the gromacs.org_gmx-users mailing list