I have to simulate a big complex with two different proteins, GNP ligand and acetate(ACT), MG and CA ions.
I’m using gromacs software and amber force field to generate the topology for the protein.
I have faced a first issue regarding ACT ion and GNP ligand because they are not found in residue topology database. But I guess AMBER supports this molecule inherently.
For that, I have used PRODRG program to generate the GNP topology first and I have followed the same procedure as explained in Benvan Lab tutoriel.
The structure I got from this phase is good but when I tried to superpose the GNP with my reference structure it seems that the built GNP coordinates are not the same. The GNP is translated.
My first question is: is PRODRG the suitable program for generation of small molecule topologies?
Should I look at the atomic charges obtained with PRODRG for GNP before merging it with the protein?
Could any on, please, suggest me other ways to succeed this critical phase.
Grateful!
Thank you very much for your reply.
I installed ACPYPE as you suggested and the other amber tools.
Unfortunately, once I run the command *acpype –i ligand.pdb * I got this error:
I tried to write prmtop using Chimera and assign charges to the ligand but I got error: Faillure running ANTECHAMBER for ligand
Do you know about the problem?
Thank you again
Anyone to help.
I need your suggestions and comments!
The ligand of the protein am using is GNP.
I think, it's a molecule classified as non-standard residue.
In that case, how canI write its topology?