minnale wrote:
So I dont boughter about those sentences because I confirmed in list
archives that I can proceed further steps and moreover em.gro file is
fine, so I went for restrain, here I am getting abnormal POPC
structure
means POPC tails are tilting, water molecules structure also disturbed.
Could you please tell me where iam doing mistake?
Not really. You haven't told us much. What parameters are you using
(for the protein and the lipids)? How did you insert your protein?
Also realize that 50 ps is *extremely* short in the realm of membrane
protein MD. Even after several hundred picoseconds things may look a
bit strange. Prevailing wisdom is that between 25-30 ns of MD are
required for the lipids to reach equilibrium (although they should look
reasonable after just a few ns, but parameters like area per headgroup
may not level off for a while).
-Justin
my pr.mdp filetitle = protein in popc restrained
define = -DPOSRES_protein
constraints = all-bonds
integrator = md
dt = 0.002 ; ps !
nsteps = 10000 ; total 50 ps.
nstcomm = 1
nstxout = 50
nstvout = 1000
nstfout = 0
nstlog = 10
nstenergy = 10
nstlist = 10
ns_type = grid
rlist = 0.9
coulombtype = PME
rcoulomb = 0.9
rvdw = 1.4
pbc = xyz
; Berendsen temperature coupling is on in three groups
Tcoupl = Berendsen
tc-grps = POPC Protein SOL_CL-
tau_t = 0.1 0.1 0.1
ref_t = 310 310 310
; Anisotropic pressure coupling is now on
Pcoupl = berendsen
pcoupltype = anisotropic
tau_p = 1.0 1.0 1.0 0 0 0
compressibility = 4.5e-5 4.5e-5 4.5e-5 0 0 0
ref_p = 1.0 1.0 1.0 0 0 0
; Energy monitoring
energygrps = POPC Protein SOL_CL-
; Generate velocites is on at 300 K.
gen_vel = yes
gen_temp = 310.0
gen_seed = 173529
Thanks in advance.
------------------------------------------------------------------------
_______________________________________________
gmx-users mailing list [email protected]
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
--
========================================
Justin A. Lemkul
Graduate Research Assistant
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin
========================================
_______________________________________________
gmx-users mailing list [email protected]
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php