Dear All,

A quick question as I have not really delved into code for gromacs ever, nor 
know anyone close whom has worked on it.

If I set up an MD simulation using a 4 protein complex, and 1 small peptide, 
plus waters, etc...and run the whole thing the proteins never move, only the 
amino acids within(constant temp RT and pressure 1 bar).

Two domains make one complex, and another two the other.  Basically, if I 
seperate the domains say 5, 10, 15 angstrom, etc...the amino acids will drift 
(the chains) towards each other, but the two large (global) protein units never 
move their center (I know I can make it work with Pull vectors, but why not in 
the simple system with a generated initial randomized velocities), I woundered 
why they are fixed in a normal run with minimal parameters?  Is there a reason 
(specific to developers), historical reason, or other?  As waters move around 
fine, and anything else added (salt, small molecules of 20-30 atoms, water) 
except the central molecule(s) of interest.

GrĂ¼sse

Stephan Watkins
-- 
Neu: GMX De-Mail - Einfach wie E-Mail, sicher wie ein Brief!  
Jetzt De-Mail-Adresse reservieren: http://portal.gmx.net/de/go/demail
-- 
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/Support/Mailing_Lists/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/Support/Mailing_Lists

Reply via email to