On Sun, Oct 04, 2020, Debarati DasGupta wrote:
>
>While applying a restraint of 10kcal/mol on my C3N ligand, why does the
>rms fluctuations look like this...( Plot of 12 lambdas C3N RMSD over
>the production phase of TI runs)???
Is this perhaps the problem that was discussed a while ago? Are you
restraining the *ligand* to stay in its original location, but allowing
the protein to move? In that case, the restraint on the ligand
essentially has no effect at all: letting the protein move away from the
ligand is equivalent to letting the ligand move away from the protein.
If you then run your cpptraj script, which superimposes the protein
frames onto a reference, the ligand will appear to be moving all over
the place.
Visualization of one of the trajctories is a must here (if you've not
already done that.) You can see what the ligand is doing. If it's
moving out of the pocket, there's nothing that cpptraj commands can do
after the fact.
If you want restraints to keep the ligand in the pocket, then need to be
relative restraints connecting the ligand and the protein, not
"absolute" ones that just keep the ligand (only) in one place.
Apologies if I'm off base here or mis-remembering things!
....dac
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sun Oct 04 2020 - 18:00:02 PDT
