[gmx-users] Membrane protein tutorial

Dallas B. Warren Dallas.Warren at pharm.monash.edu.au
Thu Apr 9 06:44:25 CEST 2009


Well, as the error says, popc isn't defined in the index file.  Seems
grompp is looking for a popc entry in the index file, and can't find it.
So, I would hazard to say, you need to define it.
 
Catch ya,

Dr. Dallas Warren
Department of Pharmaceutical Biology and Pharmacology
Pharmacy and Pharmaceutical Sciences, Monash University
381 Royal Parade, Parkville VIC 3010
dallas.warren at pharm.monash.edu.au
+61 3 9903 9167
---------------------------------
When the only tool you own is a hammer, every problem begins to resemble
a nail. 
 


________________________________

	From: gmx-users-bounces at gromacs.org
[mailto:gmx-users-bounces at gromacs.org] On Behalf Of Pawan Kumar
	Sent: Thursday, 9 April 2009 2:39 PM
	To: jalemkul at vt.edu; Discussion list for GROMACS users
	Subject: Re: [gmx-users] Membrane protein tutorial
	
	
	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/791f15f9/attachment.html>


More information about the gromacs.org_gmx-users mailing list