Error running gmx_MMPBSA as "failed when querying step5_1.tpr"

164 views
Skip to first unread message

Nitesh Mani

unread,
Oct 8, 2021, 5:05:52 AM10/8/21
to gmx_MMPBSA
Dear Mario,

I am running binding energy calculations of a protein-peptide system.The MD simulation system was prepared using CHARMM-GUI and it is showing the error:

[ERROR  ] MMPBSA_Error /usr/local/gromacs/bin/gmx editconf failed when querying step5_1.tpr.
           Check the gmx_MMPBSA.log file to report the problem.
  File "/home/ablab/anaconda3/envs/MMPBSA/bin/gmx_MMPBSA", line 8, in <module>
    sys.exit(gmxmmpbsa())
  File "/home/ablab/anaconda3/envs/MMPBSA/lib/python3.9/site-packages/GMXMMPBSA/app.py", line 95, in gmxmmpbsa
    app.make_prmtops()
  File "/home/ablab/anaconda3/envs/MMPBSA/lib/python3.9/site-packages/GMXMMPBSA/main.py", line 576, in make_prmtops
    self.FILES.mutant_receptor_prmtop, self.FILES.mutant_ligand_prmtop) = maketop.buildTopology()
  File "/home/ablab/anaconda3/envs/MMPBSA/lib/python3.9/site-packages/GMXMMPBSA/make_top.py", line 115, in buildTopology
    self.gmx2pdb()
  File "/home/ablab/anaconda3/envs/MMPBSA/lib/python3.9/site-packages/GMXMMPBSA/make_top.py", line 166, in gmx2pdb
    GMXMMPBSA_ERROR('%s failed when querying %s' % (' '.join(self.editconf), self.FILES.complex_tpr))
  File "/home/ablab/anaconda3/envs/MMPBSA/lib/python3.9/site-packages/GMXMMPBSA/exceptions.py", line 169, in __init__
    raise exc(msg + '. Check the gmx_MMPBSA.log file to report the problem.')
MMPBSA_Error: /usr/local/gromacs/bin/gmx editconf failed when querying step5_1.tpr. Check the gmx_MMPBSA.log file to report the problem.
Exiting. All files have been retained.

Logfile is attached for your reference

Kindly help in this regard.

Best Regards,
Nitesh

gmx_MMPBSA.log

Mario Ernesto

unread,
Oct 8, 2021, 11:51:50 AM10/8/21
to Nitesh Mani, gmx_MMPBSA
Hi and welcome! In the command line, you specify the smae group as receptor and ligand groups '-cg 0 0' and that's why gmx_MMPBSA fails on identify successfully the receptor and the ligand groups... Make sure to create a separated group for the receptor and one for the ligand in the index file and then use those in the command line...

hope this helps!

Mario E

CONFIDENTIALITY NOTICE: The contents of this email message and any attachments are intended solely for the addressee(s) and may contain confidential and/or privileged information and may be legally protected from disclosure. If you are not the intended recipient of this message or their agent, or if this message has been addressed to you in error, please immediately alert the sender by reply email and then delete this message and any attachments. If you are not the intended recipient, you are hereby notified that any use, dissemination, copying, or storage of this message or its attachments is strictly prohibited.

--
You received this message because you are subscribed to the Google Groups "gmx_MMPBSA" group.
To unsubscribe from this group and stop receiving emails from it, send an email to gmx_mmpbsa+...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/gmx_mmpbsa/0f2415b7-fe1d-4cd1-8e5d-86d8058ddfcan%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
Reply all
Reply to author
Forward
0 new messages