Hello everyone, I am trying to minimize the potential energy of a metalloenzyme containing Ni and Fe atoms. What is the best way to constrain (fix?) the position of the active site atoms during the geometry optimization? I have tried introducing bonds with relatively high force constants and alternatively, tried introducing a [constraints] section, but the atoms are still not staying put. Or should one use extra position restraints?
Any tips are welcome! Thank you. -- 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 gmx-users-requ...@gromacs.org.