GROMACS version:2022.5
GROMACS modification: Yes (plumed 2.9.0)
Dear gromacs users,
I am trying to pull a calcium ion from a metal-protein complex to find the unbinding energy using charmm36 ff. The simulation went smooth and the pulling force looks reasonable (mdp parameters are very similar to Justin’s umbrella sampling tutorial). However, PMF at the early stage of the pulling is very noisy. I am missing something here to understand and correct this, unfortunately still unable to catch it. Any help to understand this is appreciated.
How did you calculate the PMF? I don’t understand what it the ‘Umbrella Potential’ you are plotting.
That’s just the default title printed by gmx wham
. Perhaps a bit misleading, but what’s shown is the free energy surface.
Yes, Justin is correct, it is a free energy curve. Sorry for the misleading figure title.
I don’t know the exact settings you have used, but my assumption would be that sampling the unbinding of a calcium ion might require tighter placed sampling windows with higher pull force constant than what was used in the tutorial you based your simulations on. You would probably also need to simulate longer in each of the sampling simulations.
Have you checked the histogram overlap from gmx wham
?