Dear all,
I've been having some problems using the 'les split' functionality in ptraj. I've been running a number of LES simulations and using ptraj to split the copies both into separate trajectories and into a series of restarts.
I noticed that where my LES region includes the very first (N-terminal) residue in the complex, ptraj does not split the trajectory properly, and instead only spits out the co-ordinates of the first LES copy. If the first residue is excluded from the LES region, the 'split' works fine and co-ordinates of all copies are printed out separately.
I tried narrowing this down further and it seems that the problem lies with including the very first atom (i.e. Atom ID 1 in the PDB) in the LES region. Picking the sidechain only of the N-terminal residue results in the correct splitting of trajectories. Equally, picking the N-teminal residue of the ligand in the complex (i.e. with atom IDs > 1) also results in correct splitting.
When including the first atom, ptraj doesn't quit with any errors, it simply processes the trajectory as normal but only the co-ordinates of the first copy are printed. I've tried looking through the ptraj.c code (and others) to see why this might be the case but am by no means an expert and can't see anything obvious...
I'm running Amber 10/Ambertools 1.3 on Mac OSX 10.5.8, gcc compiled but have also tested this on our Linux-based HPC cluster (Amber 10/Ambertools 1.3, icc compiled) and seen exactly the same results. Ptraj passes all test cases on both systems.
Any advice?
Thanks for your help,
Richard
--
Richard Bradshaw
PhD Student, Institute of Chemical Biology
Room 537, Chemistry C1
Imperial College London
South Kensington Campus
SW7 2AZ
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Oct 11 2010 - 05:30:04 PDT
