You are here

Job 「№15623」has failed

4 posts / 0 new
Last post
Job 「№15623」has failed
#1

Hi

Iam a newby with Rosie, I sent a PDB and sdf file with many molecules for Rosie docking, and I got the following input error message:

Job 「№15623」has failed with the message:

molfile_to_params step failed, please double check your input files!

I searched the documentation and for what I know the input PDB and the SDF where in correct format, but of course I maybe wrong...

Could you please explain me what exactly is wrong with the input and how to correct it?

ADDITIONALLY
I think I may have found the problem, the hydrogens where missing in the SDF file, I added them with babel, it will work?

Thanks a lot,

Fabian Glaser PHD
Technion, Haifa
Israel

Category: 
Post Situation: 
Tue, 2015-07-07 00:17
fglaser

Unfortunately, that particular job has been deleted, so I can't see the details about why things failed. The fact that the molfile_to_params step failed indicates that it's an issue with the ligand input, though.

I'll mention that *by far* the most common error when doing ligand docking is using the wrong format for the ligand input. While the protein should be in PDB format, the ligand needs to be in SDF/molfile format. See http://ligand-expo.rcsb.org/reports/H/HEM/HEM_ideal.sdf for an example of what a valid SDF file looks like. -- Since you're using babel to convert things, I'm guessing that you probably have that down.

You're right that the ligand should have hydrogens added (for the appropriate pH: deprotonate your carboxylates and protonate your amines) - Using OpenBabel with the "-p" option (or Avogadro with the "Protonate for pH" option) is what I normally do, too. I don't think missing hydrogens would result in an error, though. You'd just get suboptimal docking results.

Other possible errors: Your ligand should be completely connected by bonds. If you have something like a chelated metal ion, you'll need to add a dummy bond (or an extra single bond) from the ligand to the metal ion.

You also should have three dimensional coordinates on your ligand. SDF files from molecular drawing programs sometimes have just 2D coordinates. Those won't do very well in docking, and may confuse the molfile_to_params step.

Tue, 2015-07-07 08:03
rmoretti

Hi thanks a lot for the detalied answer.

I tried to add the hydrogens with babel but still get an error, see for example run 15649. But I am puzzled since the hdrogens are there see for example this two molecules from the beginning of the file I am using the approved_hyd.sdf from the drugbank. You can see the full file on the run 15649. I guess there is still something wrong... but what?

Alternatively do you have a built in approved relatively small (thousands) of molecules to dock built in in rosie that I could use for docking?

Best regards,

Fabian
Technion
Haifa, Israel

[NO NAME]
OpenBabel07071510262D
Created with ChemWriter - http://chemwriter.com
175180 0 0 1 0 0 0 0 0999 V2000
12.8548 -2.6382 0.0000 O 0 0 0 0 0 0 0 0 0 0 0 0
13.9726 -2.5226 0.0000 O 0 0 0 0 0 0 0 0 0 0 0 0
10.1766 -3.9327 0.0000 O 0 0 0 0 0 0 0 0 0 0 0 0
11.2019 -0.7961 0.0000 O 0 0 0 0 0 0 0 0 0 0 0 0
8.7800 -1.3064 0.0000 O 0 0 0 0 0 0 0 0 0 0 0 0
16.8589 -3.2421 0.0000 O 0 0 0 0 0 0 0 0 0 0 0 0
10.3562 1.2163 0.0000 O 0 0 0 0 0 0 0 0 0 0 0 0
3.2702 4.8341 0.0000 O 0 0 0 0 0 0 0 0 0 0 0 0
2.3500 8.2734 0.0000 O 0 0 0 0 0 0 0 0 0 0 0 0
3.8213 4.2201 0.0000 O 0 0 0 0 0 0 0 0 0 0 0 0
5.2178 1.5938 0.0000 O 0 0 0 0 0 0 0 0 0 0 0 0
7.8960 2.8883 0.0000 O 0 0 0 0 0 0 0 0 0 0 0 0
7.1271 0.5358 0.0000 O 0 0 0 0 0 0 0 0 0 0 0 0
12.4834 3.1249 0.0000 O 0 0 0 0 0 0 0 0 0 0 0 0
13.1495 -4.0364 0.0000 N 0 0 0 0 0 0 0 0 0 0 0 0
11.2402 -2.9784 0.0000 N 0 0 0 0 0 0 0 0 0 0 0 0
15.2089 -3.2393 0.0000 N 0 0 0 0 0 0 0 0 0 0 0 0
10.1383 -1.7503 0.0000 N 0 0 0 0 0 0 0 0 0 0 0 0
15.6226 -2.5255 0.0000 N 0 0 0 0 0 0 0 0 0 0 0 0
12.3806 -6.3890 0.0000 N 0 0 0 0 0 0 0 0 0 0 0 0
9.2926 0.2619 0.0000 N 0 0 0 0 0 0 0 0 0 0 0 0
3.1485 7.0391 0.0000 N 0 0 0 0 0 0 0 0 0 0 0 0
4.3338 5.7884 0.0000 N 0 0 0 0 0 0 0 0 0 0 0 0
16.8613 -1.8132 0.0000 N 0 0 0 0 0 0 0 0 0 0 0 0
8.1907 1.4900 0.0000 N 0 0 0 0 0 0 0 0 0 0 0 0
5.1795 3.7761 0.0000 N 0 0 0 0 0 0 0 0 0 0 0 0
6.2814 2.5480 0.0000 N 0 0 0 0 0 0 0 0 0 0 0 0
11.3170 -7.3433 0.0000 N 0 0 0 0 0 0 0 0 0 0 0 0
12.6753 -7.7873 0.0000 N 0 0 0 0 0 0 0 0 0 0 0 0
3.8596 0.7098 0.0000 N 0 0 0 0 0 0 0 0 0 0 0 0
6.7335 6.3823 0.0000 N 0 0 0 0 0 0 0 0 0 0 0 0
5.9507 7.4636 0.0000 N 0 0 0 0 0 0 0 0 0 0 0 0
13.9701 -3.9516 0.0000 C 0 0 1 0 0 0 0 0 0 0 0 0
14.3043 -4.7058 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
13.6903 -5.2569 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
12.9766 -4.8431 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
12.5984 -3.4224 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
11.7912 -3.5925 0.0000 C 0 0 2 0 0 0 0 0 0 0 0 0
14.3838 -3.2378 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
11.5349 -4.3767 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
12.0859 -4.9907 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
9.8819 -2.5345 0.0000 C 0 0 2 0 0 0 0 0 0 0 0 0
10.4330 -3.1486 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
9.0747 -2.7046 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
8.8184 -3.4889 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
11.8295 -5.7750 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
8.0112 -3.6590 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
9.3694 -4.1029 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
9.5873 -1.1363 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
9.8436 -0.3521 0.0000 C 0 0 1 0 0 0 0 0 0 0 0 0
10.6509 -0.1820 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
16.4476 -2.5269 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
2.9756 6.2324 0.0000 C 0 0 1 0 0 0 0 0 0 0 0 0
12.0091 -0.6259 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
8.9979 1.6602 0.0000 C 0 0 2 0 0 0 0 0 0 0 0 0
9.5489 1.0461 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
12.1242 -7.1732 0.0000 C 0 3 0 0 0 0 0 0 0 0 0 0
2.1549 6.1475 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
4.8848 5.1744 0.0000 C 0 0 1 0 0 0 0 0 0 0 0 0
4.9232 2.9920 0.0000 C 0 0 1 0 0 0 0 0 0 0 0 0
1.8207 6.9019 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
9.2542 2.4444 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
4.1159 2.8219 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
3.5266 5.6183 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
2.4348 7.4528 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
5.6921 5.3445 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
6.8324 1.9340 0.0000 C 0 0 2 0 0 0 0 0 0 0 0 0
4.6285 4.3902 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
12.1793 -1.4332 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
12.8164 -0.4557 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
11.8389 0.1814 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
3.8596 2.0376 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
5.4742 2.3779 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
7.6397 2.1041 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
10.0615 2.6145 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
3.0790 1.7862 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
5.9484 6.1287 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
4.3411 1.3737 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
3.0790 0.9612 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
6.5761 1.1498 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
2.3645 2.1987 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
10.3178 3.3987 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
10.6125 2.0005 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
2.3645 0.5487 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
5.4646 6.7970 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
1.6500 1.7862 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
1.6500 0.9612 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
11.1251 3.5688 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
11.4198 2.1706 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
6.7349 7.2073 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
11.6761 2.9548 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
6.8258 -0.3862 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
13.1091 2.4275 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
11.5460 -2.0426 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
15.6997 -4.0928 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
11.1016 -1.5473 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
15.1318 -1.6721 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
8.3293 0.0588 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
4.0486 7.4380 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
4.6396 6.7242 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
17.8458 -1.8149 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
16.3705 -0.9597 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
7.8849 0.5542 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
6.1428 3.9791 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
6.5872 3.4838 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
10.6595 -6.6105 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
11.0111 -8.2791 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
12.3695 -8.7231 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
13.6386 -7.5842 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
4.1654 -0.2260 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
7.5290 5.8023 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
13.3889 -3.0532 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
15.3505 -4.9300 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
14.5796 -5.7398 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
13.8003 -6.3212 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
12.7919 -5.8381 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
11.9983 -5.2766 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
12.9187 -5.9115 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
12.1236 -2.5754 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
10.4879 -4.5974 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
11.2025 -5.3938 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
13.1329 -4.7701 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
12.8005 -5.7871 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
10.9289 -2.3139 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
8.3601 -1.9082 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
8.0277 -2.9253 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
8.4860 -4.5060 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
10.7825 -5.9957 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
11.1149 -4.9786 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
7.2966 -2.8626 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
7.6788 -4.6761 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
6.9642 -3.8796 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
9.0371 -5.1200 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
10.4164 -3.8822 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
10.0840 -4.8993 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
8.7966 -0.5728 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
10.9833 0.8350 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
9.9363 0.6144 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
4.0096 6.5076 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
8.6655 0.6431 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
1.6183 5.2218 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
1.0865 6.2054 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
5.2172 6.1915 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
4.2086 3.7884 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
0.7745 7.1263 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
1.5457 7.9359 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
8.5396 3.2408 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
9.5866 3.4614 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
3.4013 3.6183 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
3.0689 2.6012 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
6.4067 4.5481 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
6.7391 5.5651 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
7.1648 2.9511 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
11.3830 -2.1479 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
13.1964 -1.7655 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
12.4000 -2.4802 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
13.5311 -1.2520 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
13.1487 0.5614 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
13.8634 -0.2350 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
12.6352 0.8961 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
10.8218 0.5137 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
11.6182 1.2284 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
5.3731 1.3737 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
5.5291 0.9291 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
5.8615 1.9462 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
2.3645 3.2307 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
9.6286 4.1668 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
10.2919 1.0196 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
2.3645 -0.4833 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
4.4326 6.7988 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
0.7563 2.3022 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
0.7563 0.4452 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
11.4457 4.5497 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
12.1090 1.4025 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
7.5708 7.8125 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0
1 37 2 0 0 0 0
2 39 2 0 0 0 0
3 43 2 0 0 0 0
4 51 1 0 0 0 0
4 54 1 0 0 0 0
5 49 2 0 0 0 0
6 52 2 0 0 0 0
7 56 2 0 0 0 0
8 64 2 0 0 0 0
9 65 2 0 0 0 0
10 68 2 0 0 0 0
11 73 2 0 0 0 0
12 74 2 0 0 0 0
13 80 1 0 0 0 0
13 92 1 0 0 0 0
14 91 1 0 0 0 0
14 93 1 0 0 0 0
15 33 1 0 0 0 0
15 36 1 0 0 0 0
15 37 1 0 0 0 0
16 43 1 0 0 0 0
16 94 1 0 0 0 0
17 19 1 0 0 0 0
17 39 1 0 0 0 0
17 95 1 0 0 0 0
18 49 1 0 0 0 0
18 96 1 0 0 0 0
19 52 1 0 0 0 0
19 97 1 0 0 0 0
20 46 1 0 0 0 0
20 57 2 0 0 0 0
21 56 1 0 0 0 0
21 98 1 0 0 0 0
22 53 1 0 0 0 0
22 65 1 0 0 0 0
22 99 1 0 0 0 0
23 64 1 0 0 0 0
23100 1 0 0 0 0
24 52 1 0 0 0 0
24101 1 0 0 0 0
24102 1 0 0 0 0
25 74 1 0 0 0 0
25103 1 0 0 0 0
26 68 1 0 0 0 0
26104 1 0 0 0 0
27 73 1 0 0 0 0
27105 1 0 0 0 0
28 57 1 0 0 0 0
28106 1 0 0 0 0
28107 1 0 0 0 0
29 57 1 0 0 0 0
29108 1 0 0 0 0
29109 1 0 0 0 0
30 78 1 0 0 0 0
30 79 1 0 0 0 0
30110 1 0 0 0 0
31 77 1 0 0 0 0
31 90 1 0 0 0 0
31111 1 0 0 0 0
32 85 1 0 0 0 0
32 90 2 0 0 0 0
33 34 1 0 0 0 0
33 39 1 0 0 0 0
33112 1 6 0 0 0
34 35 1 0 0 0 0
34113 1 0 0 0 0
34114 1 0 0 0 0
35 36 1 0 0 0 0
35115 1 0 0 0 0
35116 1 0 0 0 0
36117 1 0 0 0 0
36118 1 0 0 0 0
37 38 1 0 0 0 0
38 16 1 0 0 0 0
38 40 1 0 0 0 0
38119 1 6 0 0 0
40 41 1 0 0 0 0
40120 1 0 0 0 0
40121 1 0 0 0 0
41 46 1 0 0 0 0
41122 1 0 0 0 0
41123 1 0 0 0 0
42 18 1 0 0 0 0
42 43 1 0 0 0 0
42 44 1 0 0 0 0
42124 1 6 0 0 0
44 45 1 0 0 0 0
44125 1 0 0 0 0
44126 1 0 0 0 0
45 47 1 0 0 0 0
45 48 1 0 0 0 0
45127 1 0 0 0 0
46128 1 0 0 0 0
46129 1 0 0 0 0
47130 1 0 0 0 0
47131 1 0 0 0 0
47132 1 0 0 0 0
48133 1 0 0 0 0
48134 1 0 0 0 0
48135 1 0 0 0 0
49 50 1 0 0 0 0
50 21 1 0 0 0 0
50 51 1 0 0 0 0
50136 1 6 0 0 0
51137 1 0 0 0 0
51138 1 0 0 0 0
53 58 1 0 0 0 0
53 64 1 0 0 0 0
53139 1 6 0 0 0
54 69 1 0 0 0 0
54 70 1 0 0 0 0
54 71 1 0 0 0 0
55 25 1 0 0 0 0
55 56 1 0 0 0 0
55 62 1 0 0 0 0
55140 1 1 0 0 0
58 61 1 0 0 0 0
58141 1 0 0 0 0
58142 1 0 0 0 0
59 23 1 0 0 0 0
59 66 1 0 0 0 0
59 68 1 0 0 0 0
59143 1 1 0 0 0
60 26 1 0 0 0 0
60 63 1 0 0 0 0
60 73 1 0 0 0 0
60144 1 6 0 0 0
61 65 1 0 0 0 0
61145 1 0 0 0 0
61146 1 0 0 0 0
62 75 1 0 0 0 0
62147 1 0 0 0 0
62148 1 0 0 0 0
63 72 1 0 0 0 0
63149 1 0 0 0 0
63150 1 0 0 0 0
66 77 1 0 0 0 0
66151 1 0 0 0 0
66152 1 0 0 0 0
67 27 1 0 0 0 0
67 74 1 0 0 0 0
67 80 1 0 0 0 0
67153 1 6 0 0 0
69154 1 0 0 0 0
69155 1 0 0 0 0
69156 1 0 0 0 0
70157 1 0 0 0 0
70158 1 0 0 0 0
70159 1 0 0 0 0
71160 1 0 0 0 0
71161 1 0 0 0 0
71162 1 0 0 0 0
72 76 1 0 0 0 0
72 78 2 0 0 0 0
75 82 2 0 0 0 0
75 83 1 0 0 0 0
76 79 1 0 0 0 0
76 81 2 0 0 0 0
77 85 2 0 0 0 0
78163 1 0 0 0 0
79 84 2 0 0 0 0
80164 1 0 0 0 0
80165 1 0 0 0 0
81 86 1 0 0 0 0
81166 1 0 0 0 0
82 88 1 0 0 0 0
82167 1 0 0 0 0
83 89 2 0 0 0 0
83168 1 0 0 0 0
84 87 1 0 0 0 0
84169 1 0 0 0 0
85170 1 0 0 0 0
86 87 2 0 0 0 0
86171 1 0 0 0 0
87172 1 0 0 0 0
88 91 2 0 0 0 0
88173 1 0 0 0 0
89 91 1 0 0 0 0
89174 1 0 0 0 0
90175 1 0 0 0 0
M CHG 1 57 1
M END
> <DATABASE_ID>
DB00014

> <DATABASE_NAME>
drugbank

> <SMILES>
CC(C)C[C@H](NC(=O)[C@@H](COC(C)(C)C)NC(=O)[C@H](CC1=CC=C(O)C=C1)NC(=O)[C@H](CO)NC(=O)[C@H](CC1=CNC2=CC=CC=C12)NC(=O)[C@H](CC1=CN=CN1)NC(=O)[C@@H]1CCC(=O)N1)C(=O)N[C@@H](CCCN=C(N)N)C(=O)N1CCC[C@H]1C(=O)NNC(N)=O

> <INCHI_IDENTIFIER>
InChI=1S/C59H84N18O14/c1-31(2)22-40(49(82)68-39(12-8-20-64-57(60)61)56(89)77-21-9-13-46(77)55(88)75-76-58(62)90)69-54(87)45(29-91-59(3,4)5)74-50(83)41(23-32-14-16-35(79)17-15-32)70-53(86)44(28-78)73-51(84)42(24-33-26-65-37-11-7-6-10-36(33)37)71-52(85)43(25-34-27-63-30-66-34)72-48(81)38-18-19-47(80)67-38/h6-7,10-11,14-17,26-27,30-31,38-46,65,78-79H,8-9,12-13,18-25,28-29H2,1-5H3,(H,63,66)(H,67,80)(H,68,82)(H,69,87)(H,70,86)(H,71,85)(H,72,81)(H,73,84)(H,74,83)(H,75,88)(H4,60,61,64)(H3,62,76,90)/t38-,39-,40-,41-,42-,43-,44-,45+,46-/m0/s1

> <INCHI_KEY>
BLCLNMBMMGCOAS-URPVMXJPSA-N

> <JCHEM_ACCEPTOR_COUNT>
18

> <JCHEM_AVERAGE_POLARIZABILITY>
131.22118283534343

> <JCHEM_BIOAVAILABILITY>
0

> <JCHEM_DONOR_COUNT>
17

> <JCHEM_EXACT_MASS>
1268.641439486

> <JCHEM_FORMAL_CHARGE>
0

> <JCHEM_FORMULA>
C59H84N18O14

> <JCHEM_GHOSE_FILTER>
0

> <JCHEM_IUPAC>
(2S)-1-[(2S)-2-[(2S)-2-[(2R)-3-(tert-butoxy)-2-[(2S)-2-[(2S)-3-hydroxy-2-[(2S)-2-[(2S)-3-(1H-imidazol-5-yl)-2-{[(2S)-5-oxopyrrolidin-2-yl]formamido}propanamido]-3-(1H-indol-3-yl)propanamido]propanamido]-3-(4-hydroxyphenyl)propanamido]propanamido]-4-methylpentanamido]-5-[(diaminomethylidene)amino]pentanoyl]-N-(carbamoylamino)pyrrolidine-2-carboxamide

> <ALOGPS_LOGP>
0.30

> <JCHEM_LOGP>
-5.17557331843203

> <ALOGPS_LOGS>
-4.65

> <JCHEM_MASS>
1269.4105

> <JCHEM_MDDR_LIKE_RULE>
1

> <JCHEM_NUMBER_OF_RINGS>
6

> <JCHEM_PHYSIOLOGICAL_CHARGE>
2

> <JCHEM_PKA>
9.818594492171766

> <JCHEM_PKA_STRONGEST_ACIDIC>
9.267832115311407

> <JCHEM_PKA_STRONGEST_BASIC>
10.816089522841512

> <JCHEM_POLAR_SURFACE_AREA>
495.8899999999997

> <JCHEM_REFRACTIVITY>
325.83880000000016

> <JCHEM_ROTATABLE_BOND_COUNT>

Sun, 2015-07-12 01:57
fglaser

One of your issues with job 15649 is that your input SDF file contains more than one molecule. Each ROSIE ligand docking run is expecting to get an SDF file with a single chemical species in it. (The multiple structures in the SDF file should be different conformers of the same chemical compound.) If you want to dock multiple ligands, you need to submit multiple jobs. If you're attempting to do high-throughput screening, you're probably better off trying to use Rosetta locally - Rosie isn't set up for HTS docking.

The more proximal cause of your error is that some of your molecules (Porfimer, Ergoloid mesylate, etc.) have '*' for element types for some of their atoms. Rosetta needs to know what the atom types are for the ligand, so having things like '*' isn't going to work. -- Also, not every element is supported. You should be fine with standard CHNOPS based compounds and things like halogens, but once you get into more esoteric metal complexes you may have issues. (The core problem is that Rosetta doesn't know how to evaluate interactions with these complexes.)

Rosetta has a small number of small molecules it knows about, but those are normally biologically based. If you're looking for drug-like molecules, you'll need to supply your own set.

Also, note that 15649 still has an issues with 2D structures. The input structures for ligand docking should already have 3D coordinates added. (The --gen3D option of OpenBabel should do a passable job of generating them.) In fact, you normally want multiple such 3D conformations, varying based on the internal degrees of freedom.

Tue, 2015-07-14 14:35
rmoretti