[gmx-users] Com Pulling

[chris.neale at utoronto.ca]

If I understand correctly, you want to apply a biasing force and you  
also want to have an "unbiased" simulation. These two things are not  
possible in combination. There are a variety of techniques available  
that will give you particular information. For example, you can  
extract the unbiased PMF from umbrella sampling (but not the unbiased  
dynamics); you can also study some things with reaction path sampling  
or slow AFM pulling, but for us to help you you will need to describe  
what you want to do.

Chris.

-- original message --

Dear All,

I was setting up things to use COM pulling for a multiprotein complex.  
  My question after going through the gromacs documentation(online) is:

Is there a way to limit either the time or distance two groups are  
pulled towards each other(reguarding only the applied force), as I  
could not find anything.  Basically, if I run something in vacuo, the  
groups (proteins) do not move globally even when constraints, etc are  
set to none.  To measure any Gibbs , delta H etc...I need the groups  
too move towards each other and form the interactions I am measuring,  
but do not want the pull force to interfear with overall calculations.

My thoughts were to do a very short run, and then either remove the  
pull force, etc from the .mdp file, and use a checkpoint to re-start  
it, but do not know if this would work or would maintain the movement  
however arbitrary in the general correct direction.

that or use such a small (.2 kcal/mol nm2) that it does not matter very much.

Sincerely,

Stephan Watkins
-- 
Sicherer, schneller und einfacher. Die aktuellen Internet-Browser -
jetzt kostenlos herunterladen! http://portal.gmx.net/de/go/chbrowser