It is not clear to me how one would do this with MD. The only thing
that I can think of doing in gromacs is to create a virtual particle
that is placed at the center of the protein and then apply forces
along the vector from this virtual atom to each of the Ca atoms. You
would need to modify the gromacs source code so that the Rg was
calculated at each step and the magnitude of each force is scaled
appropriately such that you get a harmonic potential about the desired
value of Rg. It should be easier to do with MC (although that's a ways
off for gromacs unless I've missed something).
Some people appear to have done exactly what you want with MD. I
presume that they used charmm.
http://www.pnas.org/content/94/19/10161.long (and their reference 22).
Chris.
-- original message --
Hi,
Is there any way in gromacs to use radius of gyration of a polymer
as reaction coordinate for umbrella sampling ?
Thanks
Sanku
--
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