[gmx-users] restraining to center of geometry?

David spoel at xray.bmc.uu.se
Tue Jan 27 17:14:00 CET 2004

On Sun, 2004-01-25 at 04:00, Michael Shirts wrote:
> Hello, all-
> So, supposing I had a protein, and wanted to restrain a ligand by attaching
> it with a spring to some point that wasn't actually an atom -- say it was
> the center of geometry of a number of atoms around the binding pocket, so that
> it's a point that redefines itself at each step of the simulation.
> 1) Is this possible to do using either the position restraint or distance
> restraint tools in GROMACS?
> 2) If not, how might one modify the code to put this in?
If this is dynamic, like the centre of mass of a group, you might add
something using a user group. Compute the center of mass of your user
group and call the position restraint routine. The good thing (TM) is
that you don't have to program any I/O.

> Thanks,
> Michael Shirts
> Stanford University.
> _______________________________________________
> gmx-users mailing list
> gmx-users at gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-users
> Please don't post (un)subscribe requests to the list. Use the 
> www interface or send it to gmx-users-request at gromacs.org.
David van der Spoel, PhD, Assist. Prof., Molecular Biophysics group,
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596,  	75124 Uppsala, Sweden
phone:	46 18 471 4205		fax: 46 18 511 755
spoel at xray.bmc.uu.se	spoel at gromacs.org   http://xray.bmc.uu.se/~spoel

More information about the gromacs.org_gmx-users mailing list