[gmx-users] Free energy in many conformation of molecular
Mark.Abraham at anu.edu.au
Thu Oct 12 11:42:24 CEST 2006
> I want to find a best (or stable) conformation for one molecular in
> First, I built many conformations for my molecular. Then, I run MD to get
> local minima. I want to know the free energy difference for all conformation
> in molecular. If one conformation had lowest free energy difference, I
> it is better than others.
> Now I want to know the free energy difference between conformation N and
> conformation A.(N = A, B, C …, etc). Compare the results.
> I don’t know how to do this work.
Your approach is upside-down. To evaluate free-energy differences with
explicit solvation you need to evaluate/estimate the density of states
over energy space. You can't do this with MD for some preconceived
configurations because you won't re-sample the same configurations in
accessible simulation times. See the paper below for a step in the right
F A Hamprecht, C Peter, X Daura, W Thiel and W F van Gunsteren (2001).
A strategy for analysis of (molecular) equilibrium simulations:
Configuration space density estimation, clustering, and visualization.
In: Journal of Chemical Physics, 114(5):2079-2089.
More information about the gromacs.org_gmx-users