Dear all,

I have used the PULL code in gmx 4 to unwind part of my protein. 10 runs were 
performed (all pulling the pulled group from 1.3 nm to 4.3 nm away from the 
reference). The output seems ok from the pullx.xvg and the trajectories. 
However, when using g_wham to obtain the PMF there seems to be a problem.

The options for pulling in the mdp file:
; PULLING STUFF
pull            = umbrella
pull_geometry   = distance
pull_dim        = n y n
pull_start      = yes
pull_init1      = 0.0
pull_vec1       = 0.0 1.0 0.0
pull_rate1      = 0.005 
pull_ngroups    = 1
pull_group0     = Glu33
pull_group1     = Glu81
pull_k1         = 5000
pull_nstxout    = 500
pull_nstfout    = 500

Using the command (with all the files where they should be)
 g_wham -ix -it -o -hist 
results in the following:
@    title "Umbrella potential"
@    xaxis  label "z"
@    yaxis  label "E (kJ mol\S-1\N)"
@TYPE xy
1.324382e+00    0.000000e+00
1.338367e+00    0.000000e+00
1.352352e+00    0.000000e+00
1.366337e+00    0.000000e+00
....
2.401216e+00    -1.638628e+03
2.415201e+00    -1.730001e+03
2.429185e+00    0.000000e+00
2.443170e+00    -inf
2.457155e+00    -inf
2.471140e+00    -inf
2.485125e+00    -inf
2.499110e+00    -inf
2.513095e+00    -inf
....etc

Using the command (with all the files where they should be)
 g_wham -if -it -o -hist 
results in a more reasonable plot. However, I don't understand what is on the 
x-axis. I would expect the reaction coordinate (so distance between pulled and 
reference groups) there, but this doesn't seem to be the case. 

I a post of a very similar problem in the mailing list archive of April 2009. 
It seems, however, that nobody has answered it. I wondered if anyone has a 
solution to it now?

Thanks,
Marieke
_______________________________________________
gmx-users mailing list    gmx-users@gromacs.org
http://www.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

Reply via email to