Dear Gmx Users, I am fighting with the position restraints. My system cosnsists of 6 ligands - all of them have the same topology - ligand.itp. I want to:
1. Run NPT restraining all of ligands (easy: -DPOSRES_L, ifdef POSRE_L, posre ligand1.itp 2. Restrain 5 of them and pull one of them away. Here I came across some difficulties. I used genrestr -f 6ligands.gro -n index.ndx -o posre1.itp and by index files I create an index of the ligands I wan to restrain. Each my ligands has 46 atoms. When I process to grompp: Atom index (256) in position_restraints out of bounds (1-46). This probably means that you have inserted topology section "position_restraints Shall I rename Ligands as different residue names? If my posre includes 1-46 atom all ligands are restraints? How to restrain just 5 out of 6 of them? Best, Steven
-- 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