Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

gmx_MMPBSA error on running CHARMMff examples #83

Closed
yjwang-szbl opened this issue Oct 1, 2021 · 8 comments
Closed

gmx_MMPBSA error on running CHARMMff examples #83

yjwang-szbl opened this issue Oct 1, 2021 · 8 comments

Comments

@yjwang-szbl
Copy link

I installed Ambertools21 and gmx_MMPBSA v1.4.3 on Unbutu 18.04. The test runs on amberff examples work fine, but the tests on CHARMMff examples reported weird errors.

As detailed in error_V0001.log, the error is on reading the version stamp in COM.prmtop .
"""
At line 170 of file /home/szbl/soft/amber20_src/AmberTools/src/lib/nxtsec.F
Fortran runtime error: Bad value during floating point read
Error termination. Backtrace:
#0 0x7f92ff5afea7 in ???
#1 0x7f92ff5b0a45 in ???
#2 0x7f92ff5b11fc in ???
#3 0x7f92ff67b6ae in ???
#4 0x7f92ff67e7f4 in ???
#5 0x55d7daff15e0 in ???
#6 0x55d7daff1193 in ???
#7 0x55d7dad0ede1 in ???
#8 0x55d7dacf3cd0 in ???
#9 0x55d7dacf2d32 in ???
#10 0x55d7dacf2d8e in ???
#11 0x7f92fefacbf6 in ???
#12 0x55d7dab4ad19 in ???
#13 0xffffffffffffffff in ???
File "/home/szbl/anaconda3/bin/gmx_MMPBSA", line 8, in
sys.exit(gmxmmpbsa())
File "/home/szbl/anaconda3/lib/python3.8/site-packages/GMXMMPBSA/app.py", line 98, in gmxmmpbsa
app.run_mmpbsa()
File "/home/szbl/anaconda3/lib/python3.8/site-packages/GMXMMPBSA/main.py", line 209, in run_mmpbsa
self.calc_list.run(rank, self.stdout)
File "/home/szbl/anaconda3/lib/python3.8/site-packages/GMXMMPBSA/calculation.py", line 84, in run
calc.run(rank, stdout=stdout, stderr=stderr)
File "/home/szbl/anaconda3/lib/python3.8/site-packages/GMXMMPBSA/calculation.py", line 158, in run
raise CalcError('%s failed with prmtop %s!' % (self.program,
CalcError: /home/szbl/soft/amber20/bin/sander failed with prmtop COM.prmtop!
Exiting. All files have been retained.
"""

This error can be solved if I changed the amberformat.py (line 813) by setting VERSION_STAMP = V0001.000 to VERSION_STAMP = V00001.000 , yet another error appears as detailed in error_V00001.log,

"""
File "/home/szbl/anaconda3/bin/gmx_MMPBSA", line 8, in
sys.exit(gmxmmpbsa())
File "/home/szbl/anaconda3/lib/python3.8/site-packages/GMXMMPBSA/app.py", line 107, in gmxmmpbsa
app.parse_output_files()
File "/home/szbl/anaconda3/lib/python3.8/site-packages/GMXMMPBSA/main.py", line 1060, in parse_output_files
self.calc_types[key]['delta'] = BindClass(
File "/home/szbl/anaconda3/lib/python3.8/site-packages/GMXMMPBSA/amber_outputs.py", line 1069, in init
self.delta2()
File "/home/szbl/anaconda3/lib/python3.8/site-packages/GMXMMPBSA/amber_outputs.py", line 1294, in delta2
self.com.fill_composite_terms()
File "/home/szbl/anaconda3/lib/python3.8/site-packages/GMXMMPBSA/amber_outputs.py", line 370, in fill_composite_terms
self.data[component] = self.data[key] + self.data[component]
File "/home/szbl/anaconda3/lib/python3.8/site-packages/GMXMMPBSA/amber_outputs.py", line 59, in add
raise LengthError('length mismatch in energy vectors')
LengthError: length mismatch in energy vectors
Exiting. All files have been retained.
""""
It seems to me that the internal converstion between CHARMMff top to AMBER top is problematic. Any suggestions on solving this issues?

Thank you very much.

BTW, my posts are automatically deleted in the Google Group. That is why I have to post here.

Best,

Yingjie

error_V0001.log
error_V00001.log
@Valdes-Tresanco-MS
Copy link
Owner

Hi Yingjie.

I have never seen this error and according to your description, it seems to be a problem with the structure of the input topology, that is, an unsupported version of CHARMM.
I cannot reproduce this error. If you can send me the files (*.ndx, *.trp, 10 frames *.xtc, *.top with all dependencies [top folder and *.prm], and the groups) to try to reproduce it.

Cheers!
Mario S

@Valdes-Tresanco-MS
Copy link
Owner

Apparently, Google interprets your email as spam and after approving the publications, they disappeared ... Sorry for that ...

@yjwang-szbl
Copy link
Author

yjwang-szbl commented Oct 2, 2021 via email

@Valdes-Tresanco-MS
Copy link
Owner

Hi Yingjie
Have you made progress in discovering the reason for the problem?
I have not been able to reproduce it.
I agree with you that it may be due to the charmmff in the GROMACS directory since our examples work with such content.
Have you tried with a system prepared in CHARMM-GUI?

@yjwang-szbl
Copy link
Author

Hi Mario,

I have tried both the given CHARMM examples (CHARMMff from CHARMM-GUI) inside gmx_MMPBSA_examples and my own systems (local CHARMMff directory). The errors reamin the same as previous posts. I also updated the ParmEd module to the latest version, but found no help. So the problem might be on the internal generation of AMBER topology. Would it be possible that you provide me a zip file containing successful outputs and all the intermediate files for a given example, say Protein_ligand_CHARMMff, so that I can have a detailed comparison with my local output to locate the source of errors? You can send the file to my email. Thank you very much.

@Valdes-Tresanco-MS
Copy link
Owner

I attach the example solved with CHARMM ff.
Please keep me aware of the problem.

Cheers!
Mario S.

Protein_ligand_CHARMMff.zip

@yjwang-szbl
Copy link
Author

yjwang-szbl commented Oct 9, 2021 via email

@Valdes-Tresanco-MS
Copy link
Owner

Excellent.
I was going to recommend reinstalling Ambertools, but it got ahead of me.

I close this issue now. If you have another error, please open a new issue

Cheers!
Mario S.

@laodeaman-ai laodeaman-ai mentioned this issue Jul 12, 2023
Closed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Assignees
No one assigned
Labels
None yet
Projects
None yet
Milestone
No milestone
Development

No branches or pull requests
2 participants
@Valdes-Tresanco-MS @yjwang-szbl