abdullah ahmed wrote:
Thank you for your reply,
Settings:
I am using the procedure outline by the "speptide" tutorial (i.e.
pdb2gmx, followed by editconf, genbox, grompp, and mdrun). The ".mdp"
file I've used is as follows:
; User spoel (236)
; Wed Nov 3 17:12:44 1993
; Input file
;
cpp = /usr/bin/cpp
define = -DPOSRES
constraints = none
integrator = steep
nsteps = 100
;
; Energy minimizing stuff
;
emtol = 20
emstep = 0.01
nstcgsteep = 10
nstcomm = 1
ns_type = grid
rlist = 1
rcoulomb = 1.0
rvdw = 1.0
Tcoupl = no
Pcoupl = no
gen_vel = no
coulombtype = PME
Force field: Encad all atom force field, using full sovent charges.
Local geometry is identical with the exception of the second charged
residue, which is replaced by Leucine in one structure.
Well, all you've done is a simple, restrained energy minimization, which won't
tell you much about anything. I don't see much point in using restraints during
EM anyway. Running real, unrestrained MD and obtaining suitable averages for
energies is the only way to make the conclusions I think you're after.
-Justin
------------------------------------------------------------------------
Your E-mail and More On-the-Go. Get Windows Live Hotmail Free. Sign up
now. <https://signup.live.com/signup.aspx?id=60969>
--
========================================
Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
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 gmx-users@gromacs.org
http://lists.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 gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php