Why not do a steered MD or umbrella sampling, where you start with the
ligand in the binding pocket (in its correct conformation) and gradually
pull it out.  If done correctly, you should get a nice PMF.

> Hello!
> I need to calculate the free energy of complex formation between protein A
> and ligand B. I am about to simulate in solution the molecular dynamics of
> A, B and AB complex separately and calculate the free energy of each
> system using g_lie. Then, I suppose, deltaG = G(AB) - (G(A) + G(B)).
> Is it a correct algorithm? Or maybe it's better to use "slow-growth"
> method (whith lambda = 1 for the absence of ligand in the binding pocket)
> to calculate the difference between G(AB) and G(A)?
> Thank you!
>
>
>
> _______________________________________________
> 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 [EMAIL PROTECTED]
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php
>

_______________________________________________
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 [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Reply via email to