I haven't done it myself, but both of these methods appear as if they
could be implemented out of the box in gromacs.
On the use of orientational restraints and symmetry corrections in
alchemical free energy calculations. Mobley DL, Chodera JD, Dill KA.
http://www.ncbi.nlm.nih.gov/pubmed/16965052
Absolute Binding Free Energy Calculations Using Molecular Dynamics
Simulations with Restraining Potentials. Jiyao Wang, Yuqing Deng, and
BenoƮt Roux
http://www.ncbi.nlm.nih.gov/pmc/articles/PMC1578458/
Chris
-- original message --
I'm applying COM pull on a two-domain protein by fixing the distance between
the centers of mass of each domain with an umbrella potential. However, I
find that the domains can rotate so that the inner surfaces of each domain
are not parallel any more. How can I apply a rotational restraint so that
each domain doesn't rotate during simulation? Thanks in advance.
--
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