24 January 2020 2 1K Report

I have a molecule with MG ions in the pdb file. When I am creating topology file using pdb2gmx, it does not recognize the MG ions (with gromos 54a7 ff). So I created pdb with only protein atoms and ran pdb2gmx to get a topology file (protein.itp) and protein.gro file. I used editconf for the only MG ions pdb and created MG.gro. I copied the contents from MG.gro into protein.gro and named it complex.gro . This file I used to solvate and there was no problem, but when i add ions using

gmx grompp -f ions.mdp -c solv.gro -p topol.top -o ions.tpr I get following error

Fatal error:

number of coordinates in coordinate file (solv.gro, 216418)

does not match topology (topol.top, 216496)

Please someone help me to create a topology for my molecule with Mg ions.

Thanking you.

Similar questions and discussions