Open yinglongmiao opened 6 years ago
There is no file called test_gromacs.py
. Check out the testing directory for more info: https://github.com/shirtsgroup/InterMol/tree/master/intermol/tests
P.S.: It's generally not a good idea to sudo pip install
since this could execute arbitrary code buried in the setup.py
Thanks for your reply, Christoph.
I have re-installed the package without sudo, but with the —user option. Also it runs fine with py.test -v.
Now I’m wondering whether we can convert maestro files as used by the DE Shaw’s group into the Amber .prmtop and .crd format with InterMol. What command(s) do we need for that?
All the best, Yinglong
Yinglong Miao, Ph.D. Assistant Professor Center for Computational Biology and Department of Molecular Biosciences University of Kansas 2030 Becker Drive, Suite 200D MRB Lawrence, KS 66047 miao@ku.edumailto:miao@ku.edu; 785-864-8298 (office) http://miao.compbio.ku.edu
On Oct 9, 2017, at 10:44 AM, Christoph Klein notifications@github.com<mailto:notifications@github.com> wrote:
There is no file called test_gromacs.py. Check out the testing directory for more info: https://github.com/shirtsgroup/InterMol/tree/master/intermol/tests
P.S.: It's generally not a good idea to sudo pip install since this could execute arbitrary code buried in the setup.py
— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHubhttps://github.com/shirtsgroup/InterMol/issues/340#issuecomment-335197506, or mute the threadhttps://github.com/notifications/unsubscribe-auth/AfG_nyp73P5ZFKDeBqgzSFhXFqB-H4Uzks5sqj94gaJpZM4Px0Jd.
We have basic support for the .cms
file format. Check out the readme for basic instructions and let us know if you have any specific questions!
Just found an error as follows:
python convert.py --des_in tests/desmond/stress_tests/2BP4_amber99_bulk/2BP4_amber99_bulk.cms --amber INFO 2017-10-10 17:05:31 Beginning InterMol conversion
INFO 2017-10-10 17:05:31 InterMol is research software. If you make use of InterMol in scientific publications please cite the following reference: Shirts, M.R., Klein, C., Swails, J.M. et al. J Comput Aided Mol Des (2016). doi:10.1007/s10822-016-9977-1
ERROR 2017-10-10 17:05:32 nrexcl 0 not currently supported Traceback (most recent call last): File "convert.py", line 663, in _save_amber top = pmd.load_file(oname + '.top', xyz=oname + '.gro') File "/usr/local/lib/python2.7/dist-packages/ParmEd-2.7.3-py2.7.egg/parmed/formats/registry.py", line 214, in load_file return cls(filename, *args, **kwargs) File "/usr/local/lib/python2.7/dist-packages/ParmEd-2.7.3-py2.7.egg/parmed/gromacs/gromacstop.py", line 247, in init self.read(fname, defines, parametrize) File "/usr/local/lib/python2.7/dist-packages/ParmEd-2.7.3-py2.7.egg/parmed/gromacs/gromacstop.py", line 435, in read GromacsWarning) GromacsWarning: nrexcl 0 not currently supported INFO 2017-10-10 17:05:32 Finished!
From: Christoph Klein notifications@github.com Sent: Tuesday, October 10, 2017 11:07 PM To: shirtsgroup/InterMol Cc: Miao, Yinglong; Author Subject: Re: [shirtsgroup/InterMol] Issues with installing InterMol (#340)
We have basic support for the .cms file format. Check out the readmehttps://github.com/shirtsgroup/InterMol/blob/master/README.md for basic instructions and let us know if you have any specific questions!
- You are receiving this because you authored the thread. Reply to this email directly, view it on GitHubhttps://github.com/shirtsgroup/InterMol/issues/340#issuecomment-335676535, or mute the threadhttps://github.com/notifications/unsubscribe-auth/AfG_n5Dz5GDcI3KIGI7fMbIDhU6OuFPOks5srD8GgaJpZM4Px0Jd.
Can you upload an example file that it fails on? Much easier to debug that way!
Hi Michael,
The file was “InterMol/intermol/tests/desmond/stress_tests/2BP4_amber99_bulk/2BP4_amber99_bulk.cms”, a test Desmond file I found within the InterMol package.
Eventually, I would like to convert the maestro files used by the DE Shaw’s group for their Anton simulations (available here: http://pubs.acs.org/doi/suppl/10.1021/acs.jctc.7b00172) to the Amber format if it’s possible ...
Thanks! Yinglong
On Oct 12, 2017, at 11:57 PM, Michael Shirts notifications@github.com<mailto:notifications@github.com> wrote:
Can you upload an example file that it fails on? Much easier to debug that way!
— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHubhttps://github.com/shirtsgroup/InterMol/issues/340#issuecomment-336350451, or mute the threadhttps://github.com/notifications/unsubscribe-auth/AfG_n1IZroh5iRvX_6U1Nd7mfsVo5WBVks5sru25gaJpZM4Px0Jd.
I'll look into converting that over the weekend.
That particular file was converted from amber to desmond and back if I recall. It would be useful if you tried to do a conversion of a file you were actually interested in which is Desmond native. It might work.
OF course, we want to fix the bug, but making sure it works for your usage is the first priority.
Then it would be wonderful if we can make InterMol work for the maestro files you can download at: http://pubs.acs.org/doi/suppl/10.1021/acs.jctc.7b00172
Right now, I run into the following error:
python convert.py --des_in dmso.mae --amber INFO 2017-10-12 22:24:03 Beginning InterMol conversion
INFO 2017-10-12 22:24:03 InterMol is research software. If you make use of InterMol in scientific publications please cite the following reference: Shirts, M.R., Klein, C., Swails, J.M. et al. J Comput Aided Mol Des (2016). doi:10.1007/s10822-016-9977-1
Traceback (most recent call last):
File "convert.py", line 811, in
Thanks, Yinglong
On Oct 13, 2017, at 10:21 PM, Michael Shirts notifications@github.com<mailto:notifications@github.com> wrote:
I'll look into converting that over the weekend.
That particular file was converted from amber to desmond and back if I recall. It would be useful if you tried to do a conversion of a file you were actually interested in which is Desmond native. It might work.
OF course, we want to fix the bug, but making sure it works for your usage is the first priority.
— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHubhttps://github.com/shirtsgroup/InterMol/issues/340#issuecomment-336606132, or mute the threadhttps://github.com/notifications/unsubscribe-auth/AfG_nzDDTb5yZerOY1gN1PVryJWHX-iwks5ssCi1gaJpZM4Px0Jd.
OK, I've downloaded that one. Lots of grant writing this weekend, so if I can't get it to work with a bit of a tweak, it might take a while.
Sure, thanks a lot!
On Oct 13, 2017, at 10:41 PM, Michael Shirts notifications@github.com<mailto:notifications@github.com> wrote:
OK, I've downloaded that one. Lots of grant writing this weekend, so if I can't get it to work with a bit of a tweak, it might take a while.
— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHubhttps://github.com/shirtsgroup/InterMol/issues/340#issuecomment-336606985, or mute the threadhttps://github.com/notifications/unsubscribe-auth/AfG_n-WnH2_5WPNIIDvjtLeUHQF9HVdVks5ssC1egaJpZM4Px0Jd.
Michael, hope all went well with your grant writing!
Any good news on converting the file?
Thanks, Yinglong
On Oct 13, 2017, at 10:41 PM, Michael Shirts notifications@github.com<mailto:notifications@github.com> wrote:
OK, I've downloaded that one. Lots of grant writing this weekend, so if I can't get it to work with a bit of a tweak, it might take a while.
— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHubhttps://github.com/shirtsgroup/InterMol/issues/340#issuecomment-336606985, or mute the threadhttps://github.com/notifications/unsubscribe-auth/AfG_n-WnH2_5WPNIIDvjtLeUHQF9HVdVks5ssC1egaJpZM4Px0Jd.
Not done yet, hopefully tomorrow. Apologies.
Thanks for the update, Michael. I have been waiting for converting these files so that we can start running the simulations here. Hope it’s not too much effort for the coding … ;-)
On Oct 26, 2017, at 2:41 AM, Michael Shirts notifications@github.com<mailto:notifications@github.com> wrote:
Not done yet, hopefully tomorrow. Apologies.
— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHubhttps://github.com/shirtsgroup/InterMol/issues/340#issuecomment-339579813, or mute the threadhttps://github.com/notifications/unsubscribe-auth/AfG_nzBxLEEHOJGK5iii6-cA9cxhpf4vks5swDexgaJpZM4Px0Jd.
Hi, Yinglong-
I think I have a workaround for now. The problem is actually ParmEd, because the number of atoms to exclude is not properly being sent to the GROMACS file.
The workaround is:
1) convert desmond to GROMACS. This is actually what InterMol currently does anyway.
2) manually fix the GROMACS files. Wherever you see an entry of form
[ moleculetype ] molecule_name 0
Then change the 0 to 3 (which it should be for AMBER or GAFF).
3) Convert GROMACS to AMBER.
Let me know if that works. I'll look at the underlying code reason, but it might take a bit longer.
Hi Michael,
I tried your suggestion on the dmso.mae, but still get the error for converting to gromacs. This step worked for example .cms file included in the intermol/tests, but conversion to amber showed "TypeError: AmberParm does not support all of the parameters defined in the input Structure", and no output for Amber files.
I have attached the program log below. Any further suggestions for fixing dmso.mae and step 3 - conversion to Amber will be much appreciated!
Thanks,
Yinglong
python ~/programs/InterMol/intermol/convert.py --des_in dmso.mae --gromacs INFO 2017-11-10 13:41:08 Beginning InterMol conversion
INFO 2017-11-10 13:41:08 InterMol is research software. If you make use of InterMol in scientific publications please cite the following reference: Shirts, M.R., Klein, C., Swails, J.M. et al. J Comput Aided Mol Des (2016). doi:10.1007/s10822-016-9977-1
Traceback (most recent call last):
File "/home/y147m353/programs/InterMol/intermol/convert.py", line 811, in
! ~/programs/InterMol/intermol/convert.py --des_in ~/programs/InterMol/intermol/tests/desmond/stress_tests/2BP4_amber99_bulk/2BP4_amber99_bulk.cms --gromacs INFO 2017-11-10 13:27:57 Finished!
! vim 2BP4_amber99_bulk_converted.* [ moleculetype ] Molecule_1 3
! ~/programs/InterMol/intermol/convert.py --gro_in 2BP4_amber99_bulk_converted.gro 2BP4_amber99_bulk_converted.top --amber ERROR 2017-11-10 13:31:24 AmberParm does not support all of the parameters defined in the input Structure Traceback (most recent call last): File "/home/y147m353/programs/InterMol/intermol/convert.py", line 664, in _save_amber top.save(oname + '.prmtop', overwrite=True) File "/usr/local/lib/python2.7/dist-packages/ParmEd-2.7.3-py2.7.egg/parmed/structure.py", line 1533, in save s = amber.AmberParm.from_structure(self) File "/usr/local/lib/python2.7/dist-packages/ParmEd-2.7.3-py2.7.egg/parmed/amber/_amberparm.py", line 326, in from_structure raise TypeError('AmberParm does not support all of the ' TypeError: AmberParm does not support all of the parameters defined in the input Structure INFO 2017-11-10 13:31:24 Finished!
From: Michael Shirts notifications@github.com Sent: Thursday, November 9, 2017 11:27 PM To: shirtsgroup/InterMol Cc: Miao, Yinglong; Author Subject: Re: [shirtsgroup/InterMol] Issues with installing InterMol (#340)
Hi, yinglong-
I think I have a workaround for now. The problem is actually ParmEd, because the number of atoms to exclude is not properly being sent to the GROMACS file.
The workaround is:
[ moleculetype ] molecule_name 0
Then change the 0 to 3 (which it should be for AMBER or GAFF).
Let me know if that works. I'll look at the code, but it might take a bit longer.
- You are receiving this because you authored the thread. Reply to this email directly, view it on GitHubhttps://github.com/shirtsgroup/InterMol/issues/340#issuecomment-343378429, or mute the threadhttps://github.com/notifications/unsubscribe-auth/AfG_n06obtf4_5G5-L3y1TMxTv6LiWrAks5s097GgaJpZM4Px0Jd.
For the parmed problem, this is an issue with the ParmEd libraries that InterMol uses to convert from GROMACS to AMBER - it simply can't take all representations of dihedrals. However, I know that they were looking to resolve this issue. Your best bet is to check with ParmEd and see if the have resolved this issue. I will check the dmso.mae problem.
Hi Michael,
That sounds good. Please let me know when you find a solution to the dmso.cms problem. I will also check with the latest version of ParmEd regarding the conversion from GROMACS to AMBER …
Thanks, Yinglong
Yinglong Miao, Ph.D. Assistant Professor Center for Computational Biology and Department of Molecular Biosciences University of Kansas http://miao.compbio.ku.edu
On Nov 13, 2017, at 12:33 AM, Michael Shirts notifications@github.com<mailto:notifications@github.com> wrote:
For the parmed problem, this is an issue with the ParmEd libraries that InterMol uses to convert from GROMACS to AMBER - it simply can't take all representations of dihedrals. However, I know that they were looking to resolve this issue. Your best bet is to check with ParmEd and see if the have resolved this issue. I will check the dmso.cms problem.
— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHubhttps://github.com/shirtsgroup/InterMol/issues/340#issuecomment-343826709, or mute the threadhttps://github.com/notifications/unsubscribe-auth/AfG_n3XY62aGcZxU4p87-5lSo23i7aFjks5s1-LDgaJpZM4Px0Jd.
I have confirmed the error: something to do with the way the blocks are set up in this particular MAE file. It will take a bit of time to fix correctly.
On Mon, Nov 13, 2017 at 8:02 AM, yinglongmiao notifications@github.com wrote:
Hi Michael,
That sounds good. Please let me know when you find a solution to the dmso.cms problem. I will also check with the latest version of ParmEd regarding the conversion from GROMACS to AMBER …
Thanks, Yinglong
Yinglong Miao, Ph.D. Assistant Professor Center for Computational Biology and Department of Molecular Biosciences University of Kansas http://miao.compbio.ku.edu
On Nov 13, 2017, at 12:33 AM, Michael Shirts <notifications@github.com< mailto:notifications@github.com>> wrote:
For the parmed problem, this is an issue with the ParmEd libraries that InterMol uses to convert from GROMACS to AMBER - it simply can't take all representations of dihedrals. However, I know that they were looking to resolve this issue. Your best bet is to check with ParmEd and see if the have resolved this issue. I will check the dmso.cms problem.
— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHubhttps://github.com/ shirtsgroup/InterMol/issues/340#issuecomment-343826709, or mute the threadhttps://github.com/notifications/unsubscribe- auth/AfG_n3XY62aGcZxU4p87-5lSo23i7aFjks5s1-LDgaJpZM4Px0Jd.
— You are receiving this because you commented. Reply to this email directly, view it on GitHub https://github.com/shirtsgroup/InterMol/issues/340#issuecomment-343945853, or mute the thread https://github.com/notifications/unsubscribe-auth/AEE31OpAa_fuMO4U7L2DdwREacQt2CnIks5s2FoZgaJpZM4Px0Jd .
Can we perhaps check with Albert Pan to see how they came up this different version of MAE file, and is there a way to generate a file that’s compatible with InterMol? Thanks!
On Nov 13, 2017, at 9:50 AM, Michael Shirts notifications@github.com<mailto:notifications@github.com> wrote:
I have confirmed the error: something to do with the way the blocks are set up in this particular MAE file. It will take a bit of time to fix correctly.
On Mon, Nov 13, 2017 at 8:02 AM, yinglongmiao notifications@github.com<mailto:notifications@github.com> wrote:
Hi Michael,
That sounds good. Please let me know when you find a solution to the dmso.cms problem. I will also check with the latest version of ParmEd regarding the conversion from GROMACS to AMBER …
Thanks, Yinglong
Yinglong Miao, Ph.D. Assistant Professor Center for Computational Biology and Department of Molecular Biosciences University of Kansas http://miao.compbio.ku.edu
On Nov 13, 2017, at 12:33 AM, Michael Shirts notifications@github.com<mailto:notifications@github.com< mailto:notifications@github.com>> wrote:
For the parmed problem, this is an issue with the ParmEd libraries that InterMol uses to convert from GROMACS to AMBER - it simply can't take all representations of dihedrals. However, I know that they were looking to resolve this issue. Your best bet is to check with ParmEd and see if the have resolved this issue. I will check the dmso.cms problem.
— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHubhttps://github.com/ shirtsgroup/InterMol/issues/340#issuecomment-343826709, or mute the threadhttps://github.com/notifications/unsubscribe- auth/AfG_n3XY62aGcZxU4p87-5lSo23i7aFjks5s1-LDgaJpZM4Px0Jd.
— You are receiving this because you commented. Reply to this email directly, view it on GitHub https://github.com/shirtsgroup/InterMol/issues/340#issuecomment-343945853, or mute the thread https://github.com/notifications/unsubscribe-auth/AEE31OpAa_fuMO4U7L2DdwREacQt2CnIks5s2FoZgaJpZM4Px0Jd .
— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHubhttps://github.com/shirtsgroup/InterMol/issues/340#issuecomment-343961752, or mute the threadhttps://github.com/notifications/unsubscribe-auth/AfG_nzTG2FYm6HtL9ClMFFbtHIzt9Ec9ks5s2GVWgaJpZM4Px0Jd.
Hello,
I tried to install InterMol, but it kept giving me the following errors:
====
Similar error was found for the convert.py. Would you have any suggestions?
Thanks, Yinglong