Replies: 5 comments
-
Hiya! The problem is most likely that your molecule has undefined stereochemistry - it's difficult to be more specific without seeing the PDB you're loading. The OpenFF Toolkit requires all molecules to have fully defined stereochemistry to reduce errors - we've made a virtue of being very strict about what we can put in a You can try inferring stereochemistry from the PDB coordinates, or you can create a |
Beta Was this translation helpful? Give feedback.
-
Thanks so much for the quick reply, I think from PDBfile, some bond information is missing, so I tried to generate ff using another way,but come with this error,is this still means the stereochemistry information is missing so that no tool can be used to generate forcefield?or it is because the ligand is name is unknown acid, so it is not able to generate forcefield?
|
Beta Was this translation helpful? Give feedback.
-
It looks like you've successfully got the molecule into the Toolkit but AmberTools doesn't like it for some reason. Could you provide the output of |
Beta Was this translation helpful? Give feedback.
-
here is the conda list output, I am refering to this link
below is the mol2 file I load in the code in text form
I have no idea about how to fix that,now . Maybe I am having too many transformations so it is hard for the toolkit to continue the job? |
Beta Was this translation helpful? Give feedback.
-
Something is wrong with how RDKit parses that file and/or how the RDKit molecule is prepared. If all you need is the molecule loaded up from SMILES anyway, the toolkit offers a more direct way of handling that in addition to many other pathways. I can assign partial charges from that SMILES pattern without issue using the most recent version of OpenFF tools and dependencies, although it takes a while since AM1 is slow.
Note also that you're using an old version (0.10.3) of the toolkit which won't be supported for much longer. Newer versions include a plethora of bugfixes useful features, such as loading PDB files directly. |
Beta Was this translation helpful? Give feedback.
-
Hello everyone,
I am creating Molecule object from a rdkit mol,with a conformer file,ligand.mol2 and a protein file 5x1a.pdb, after merging it to a complex file 5x1a_complex.pdb, it is impossible to generate Molecule object in this way.
Here is the error. I am having the openff-toolkit==0.10.7,how could I generate the Molecule object correctly.
Beta Was this translation helpful? Give feedback.
All reactions