On 4/30/12 12:26 PM, Hovakim Grabski wrote:
Dear Gromacs users,

I want to set up an umbrella sampling simulation of an Hypericin dimer.
I have created position restraint filesfor hypericin with genrestr and then
use it with the topology file, so I want to know if it is just restraining one
of the Hypericin Molecules or both of them?

That depends on how your topology is set up. Position restraints are applied per moleculetype, so if you have defined one moleculetype for this molecule and simply specified two of those molecules in the .top, it is applied globally.

If the two molecules are restrained , how can I restrain just one of the
Hypericin molecules?

In this case you would need two distinct [moleculetype] directives in the topology, though they would be identical. Only one of them would specify a [position_restraints] directive.

For running umbrella sampling of two molecules in water, there is no need for position restraints at all. The only biasing potential you need is the one maintaining a distance between the two molecules.

-Justin

--
========================================

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/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