On Sat, Aug 29, 2015, Suchetana Gupta wrote:
> My protein is a homodimer of some ~100 amino acids. I wish to calculate its
> energy required for dimerisation. The process that I am following is:
> 1. Take the dimer. Run 40ns MD. Calculate MMGBSA free energy.
> 2. Take a monomer. Run 40ns MD. Calculate MMGBSA free energy.
> Dimerisation energy=1-(2*2)
This sounds correct to me. Make sure that you are including an estimate of
the configurational entropy for both calculations.
>
> Now doing this, I am getting less dimerisation energy for a wild type
> protein than the drug resistant mutant one.
Do you know if this is right or wrong (i.e. compared to experiment)? If the
drug binds in the dimer interface, then your result would be consistent with
the mutant having drug resistance (since the binding pocket is less exposed
in the mutant than in the WT protein.)
> Moreover, the drug resistant variety (both monomer as well as dimer) is
> showing stabler energy value than WT counterpart.
You cannot compare (say) the free energy of the WT monomer with that of the
mutant monomer. They are not isomers of each other, and hence molecular
mechanics gives no information about their free energy difference.
A more accurate calculation (at least in principle) would be to use TI to
mutate the WT into the mutant sequence. Do this twice: once in the monomer,
and once in the dimer, and construct a thermodynamics cycle to generate
relative binding free energies.
...dac
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sat Aug 29 2015 - 05:00:05 PDT
