Dear all,

I am calculating the mutation free energies using TI method for a small
globulin protein which contains about 40 residues.
I meet convergence problems in the FE calculations:
(1) when simulations start with different velocities for the some lambda
points, the variance of dGdlamba is big and comparable with ddG (several
kJ/mol). That makes the calculation not well reproducible.
(2) the error estimates using block averaging method is also large. The
error is comparable calculated ddG ( ~4 kJ/mol ). Thus the results seems not
very reliable.

By reading, I understand this is due to inefficient sampling for mutations
in proteins.
I attach one of my dgdl plots for LJ mutation(W->A) step (charges for
mutated atoms are removed in previous step). I try to increase the dgdl data
collection time to 2ns, which is relative long compared to time used in
protein mutation in the references. But the large variance in lambda=0.7,
0.8, 0.9 seems not get improved. The error is also not well reduced.
The follwing are the dG results and error in kJ/mol.
dG(500ps)          =  1.01       error = 3.92
dG(500ps_run1) =  -0.40       error = 3.83
dG(1ns)             =   2.21       error = 3.67
dG(2ns)             =   2.44       error = 3.22
dG(2ns_run1)     =  -1.02       error = 2.48



Can anybody give some suggestions on how to improve the convergence and
reduce the statistical error?

I appreciate anyone with free energy calculation experience can give some
advice.



Thanks

Qiang

<<attachment: dgdl.jpg>>

_______________________________________________
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