shirtsgroup / InterMol

Conversion tool for molecular simulations
MIT License
187 stars 54 forks source link

Desmond conversion #344

Open tjgdamen opened 6 years ago

tjgdamen commented 6 years ago

When I try to convert a Desmond .cms file to any other format I get the following error:

Traceback (most recent call last): File "convert.py", line 811, in main() File "convert.py", line 226, in main system, prefix = _load_desmond(cms_in=cms_in) File "convert.py", line 589, in _load_desmond system = des.load(cms_file=cms_in) File "C:\ProgramData\PyMOL\lib\site-packages\intermol-0.1.0.dev0-py2.7.egg\intermol\desmond\desmond_parser.py", line 39, in load return parser.read() File "C:\ProgramData\PyMOL\lib\site-packages\intermol-0.1.0.dev0-py2.7.egg\intermol\desmond\desmond_parser.py", line 1141, in read self.load_ffio_block(molname, self.ffio_blockpos[i], self.fmct_blockpos[i+1]-1) File "C:\ProgramData\PyMOL\lib\site-packages\intermol-0.1.0.dev0-py2.7.egg\intermol\desmond\desmond_parser.py", line 903, in load_ffio_block self.sysDirectivetype File "C:\ProgramData\PyMOL\lib\site-packages\intermol-0.1.0.dev0-py2.7.egg\intermol\desmond\desmond_parser.py", line 576, in parse_bonds new_bond = self.create_forcetype(self.desmond_bonds[key], atoms, params) KeyError: 'INVALID'

I have read several topics with seem to have similar issues, but never found a solution for this error.

Thanks in advance for helping me out.

mrshirts commented 6 years ago

I can try to take a look, if you could send an example file that fails (the smaller, the better).

tjgdamen commented 6 years ago

Dear Dr. Michael Shirts,

Thanks a lot for your quick reply. Because of a National Holiday in The Netherlands last Friday I could only come back to you now; hopefully you are able to help me out.

Please find attached the .cms file created by Desmond and zipped with 7-zip to reduce its size. As an example I did a 1 ns MD simulation of pdb 5FNU. When I prepared the system, I removed all crystal waters before solvation in SPC waters. I used the Desmond version as is implemented in Maestro of Schrodinger Suite 2018-1.

When I try to convert the .cms file to gromacs (or any other format), I got the following error:

[cid:image001.png@01D3E08E.B4547C90]

BTW, I am running the convert.py script in the most recent Anaconda (python).

Thanks in advance for taking a deeper look at it and hopefully you can help me out!

Very kind regards,

Eddy

Eddy Damen, Ph.D. Group Leader Comp. Chem. & Drug Design


Mercachem Kerkenbos 1013, 6546 BB Nijmegen P.O. Box 6747, 6503 GE Nijmegen The Netherlands phone: +31(0)24-3723307 fax: +31(0)24-3723305 e-mail: eddy.damen@mercachem.com


LEGAL NOTICE Unless expressly stated otherwise, this message is confidential and may be privileged. It is intended for the addressee(s) only. Access to this e-mail by anyone else is unauthorised. If you are not an addressee, any disclosure or copying of the contents of this e-mail or any action taken (or not taken) in reliance on it is unauthorised and may be unlawful. If you are not an addressee, please inform the sender immediately. All outgoing messages have been virus checked by sender. Sender will not be liable for and/or in connection with all and any possible consequences and/or damage resulting from this e-mail message, such as damage caused by viruses. Mercachem is a registered trade name of Mercachem Holding B.V.. Mercachem is made up of four subsidiaries: Mercachem b.v., Mercatorial b.v., Mercachem Process Research b.v. and Mercaleads b.v. All statutory seats are in Nijmegen. The respective Trade Register numbers for these corporations are as follows: 10147906, 09127273, 09117608, 09155544 and 55343104 registered in Arnhem. A fifth subsidiary: Mercachem Prague s.r.o. has its statutory seat in Prague, CZ. Its registration number with the commercial court is 032 40 878.


From: Michael Shirts [mailto:notifications@github.com] Sent: donderdag 26 april 2018 22:31 To: shirtsgroup/InterMol InterMol@noreply.github.com Cc: Eddy Damen Eddy.Damen@Mercachem.Com; Author author@noreply.github.com Subject: Re: [shirtsgroup/InterMol] Desmond conversion (#344)

I can try to take a look, if you could send an example file that fails (the smaller, the better).

— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHubhttps://github.com/shirtsgroup/InterMol/issues/344#issuecomment-384779259, or mute the threadhttps://github.com/notifications/unsubscribe-auth/Ak9dW6BfUiOCg8pS1VapWnud6eZf22ZPks5tsi6SgaJpZM4TlBKl.

shansun7994 commented 6 years ago

I ran into the same problem. Do we have a solution now?

mehdim390 commented 6 years ago

Dear all,

Any solution for this issue with Maestro (Schrodinger Suite 2018-1) ?

mrshirts commented 6 years ago

I apologize for the lack of response. We unfortunately don't have the manpower in the lab to bugfix this semester.

Thuravi commented 2 years ago

Dear Mr. Shirts,

Is there any updates on this? I am still facing the same error. It would be really helpful if you could reply regarding this issue. Thanks in advance.

Best, Santhosh