Laura,
if your organic ligand left its binding pocket during the simulation,
this is no "PBC effect", but the result of your simulation. (If this is
unexpected or not in accordance with experimental findings, maybe your
setup or your ligand parameters are not right.)
Perhaps you could use cpptraj to unwrap or autoimage your trajectory.
If VMD displays "broken" bonds in a structure, it may just be a
visualization artifact, because the distance between the two atoms is
larger than the program-specific built-in cutoff for a bond. This does
not mean, that the bond broke and maybe formed again.
Concerning a MMPBSA analysis:
This post-production analysis makes only sense for a valid trajectory,
i.e. when the ligand stays intact and bound to the protein.
Maybe this helps.
Best,
Anselm
Bioinformatik | NHR.FAU
Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU
Germany
Am 02.11.2022 um 16:19 schrieb laura zhang via AMBER:
> Dear all,
>
> I did AMBER simulations on a protein binding with a small organic molecule. Checking the trajectory, I can see that the small organic compounds had PBC effect that it departured from the protein. However, when I used pbc unwrap in vmd program to unwrap the trajectory, the small compound side chains broke here and there. I am wondering how I can unwrap the amber trajectory correctly.
> Another question is, under such kind of condition, when I did mmgbsa binding energy calculation on the protein binding with the small organic molecule, will the result be reliable? what should I do in order to get reliable binding energy result? Thank you very much.
>
> Laura
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Nov 02 2022 - 09:00:03 PDT