On Tue, 2014-05-20 at 13:21 +0000, Aronica, Pietro wrote:
> I have run the same calculation with a different but similar system (it has a staple between two residues) and I see the same error as before: when it generates the trajectory for the mutated complex, the newly created alanine has unrealistic dynamics. Furthermore, this happens in the same residue (146) while not happening in the ones that didn't have an error even in systems without the staple (143-145 and 147). I have not checked the other residues that give this error, but I expect them to show the same behaviour.
> I have then run the same calculation on an older version of AmberTools. I'm not sure which one it is, but the command
> MMPBSA --version
> returns 2.7. This calculation runs fine, generating a completely reasonable trajectory for the mutated complex and giving no error whatsoever.
> Clearly, the error is within the new version of the protocol, which, in some consistent cases, fails to give proper results.
> I hope this can be fixed somehow.
I've run your example and I have reproduced the problem. The alanine
scanning code has not changed appreciably since it was written ~5 years
ago. I still haven't been able to track down the root cause, but when
it is fixed it is likely that only AmberTools 14 will receive the
update.
Thanks,
Jason
--
Jason M. Swails
BioMaPS,
Rutgers University
Postdoctoral Researcher
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue May 20 2014 - 07:00:03 PDT
