I have been working on post modifying my trajectories using trjconv. My protein has two zinc ions, and I have used PBC (explicit waters). After a week of attempting to get my simulations to behave due to boundary effects, I have finally got my protein to not fall apart, and even managed to get the waters to stay within the box (visually using VMD). However, the problem now is my zinc ion ligands. To simplify my problem with the zincs, I removed the waters (using groups in my index file). I then created a group in the index called "protein_ZN" which contain protein and two zinc atoms..this did work..Then I did the following:
trjconv -s pr.trp -f pr13.trr -pbc nojump -n index.ndx -o protein13.trr: The output group was "protein_zn". then: trjconv -s pr.trp -f protein13.trr -n index.ndx -fit rot+trans -o protein13A.trr When I visualized the traj. in VMD it showed that this worked, but only for ONE zinc ligand ion...the other left the box with ligand to sulfur bond on 30-40 angstroms. Can someone tell me what is going on here...??? Thanks, Steve -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to [email protected].
