Hello,
I guess you should add your ligand.itp to aminoacids.rtp, and then use
pdb2gmx to generate a combined topology file.
After that, you can use distance or angle restraints between the protein
and the ligand.
Maybe you can try to parmed from Amber to combine your protein and
ligand topology. I never tried, but it might work.
All the best,
Qinghua
On 06/05/2018 09:50 AM, Qing Liu wrote:
Hello, everyone,
I just used Gromacs-5.1.4 to run a MD with distance and angle restraints. These restraints are
between a protein and a ligand. However, the ligand topology is included in the "ligand.itp", and
its atom numbers begin from 1. This is same as the protein atom numbers. As a result, Gromacs can not
distinguish the ligand atoms and protein atoms in the sections "distance_restraints" and
"angle_restraints". I also try to merge the topology of protein and ligand into a single file,
where the ligand atom numbers are renumbered to be consecutive from the last protein atom. But this also does
not work. Can you help me? Thank you!
--
Best wishes,
------------------------------------------------------------
Qing Liu
Fudan Univ.
E-mail: 15110700...@fudan.edu.cn
--
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.
