On Tue, Nov 03, 2020, SATYAJIT KHATUA wrote:
>
>Thanks a lot. If you don't mind, I am asking a different scenario with two
>protein chains bound to each other through a metal center zinc & I would
>like to calculate the protein-protein interaction energy using MMPBSA.py,
>then how should I proceed?
To me, this seems like a really difficult problem. First, estimating
the energetics of zinc-ligand interactions is difficult via force
fields. You could look at the literature for trying to calibrate Mg --
ligand interactions and think about following that path.
Second, what is the reaction you are looking at? Is is 2P + Zn ->
complex (where P is protein)? Or is it P-Zn + P -> complex, where the
Zinc is pre-coordinated to one of the two proteins? The first, a
tri-molecular reaction, is outside the scope of MMPBSA.
Third, in forming the complex, presumably a water ligand on the zinc is
being replaced by coordination to some atom of a protein. Estimating this
interaction energy by continuum dielectric theories is very sensitive to
the paramters you choose, so much so that people rarely try this for
divalent ions.
That said, you might learn things by studying this, but you should
expect to put a lot of effort into the project. Of course, that is how
simulation science advances....
...dac
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue Nov 03 2020 - 06:00:03 PST