Dear gromacs users, I am attempting to simulate a protein with an active site containing to metals (Ni and Fe). I am trying to mantain the active site structure constraining bonds, angles and dihedrals. When I simulate the active site alone in water, it remains stable. When I insert the active site in the protein's scaffold, unfortunately, the simulation crashes because of large vibration in bond lenght and angle amplitude in the active site. I am not able to detect big clashes or any structural distortion.
I tried to play with force constants, temperature, time-step, constrains, etc.etc. but I am not able to mantain the simulation stability. I am wondering whether exists a tool to dissect the energetic contributions ( including bonded parameters ) permitting me to identify the problem giving rise the crash. 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