Dear gromacs user, I noted a strange behavior in topology generation workflow.
I have a protein whose active site contains a Ni atom linked to 4 Cys. I inserted the 4 bonds by specbonds.dat. When I run pdb2gmx, I have inspected the Linking messages and they fit with the expected topology, I mean pdb2gmx correctly detects the four bonds nickel-sulphur. So I run grompp and the topol.tpr has been correctly generated. Then I compared topol.top and topol.tpr bonds and I noticed that in topol.top there a lot of fake-bonds, i.e Ni-hydrogen, Ni-Carbonium. I dumped topol.tpr through gmxdump, and I noticed that these bonds are no more present and only "right" bonds are present. So, I am wondering why the fake bonds are inserted in topol.top and how grompp detectes these bonds. Francesco
-- gmx-users mailing list gmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. Can't post? Read http://www.gromacs.org/Support/Mailing_Lists