-
Notifications
You must be signed in to change notification settings - Fork 57
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.
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
Unexpected generated molecule with double cuts MMPDB #54
Comments
Hi Alice, I am not sure whether I understand your issue: The 3rd and the 5th SMILES string in your question are exactly the same - if you start with Without knowing the details, I assume you the database you create indexes a double-cut transformation that replace one nitrogen heteorcycle with a single-fluorinated form of the other. The constant parts are one *F, and the large aromatic system. If you apply this transformation to yoru starting SMILES (e.g. the 3rd SMILES), you woudl expect to get exactly the fourth SMILES. If this is not clear, please post the transformation SMILES that is applied to produce the molecule that surprises you. Bests, |
Hi Christian, I must have copy pasted the same molecules twice by mistake, sorry for the confusion! the expected molecules would be Fc1cnc2nc(CNc3cc(-c4ccccc4)c(-c4ccccc4)c4ccccc34)cn2c1. And I agree that the 4th smiles has the correct transformation, but it is attached the wrong way around! Sorry again and I hope this clarifies my question! |
Hi Alice, can you share the command you use to get the new compounds created and the full output if you run it with the --explain or --save-details option? Thanks, |
Sure :) mmpdb fragment example_dataset_double_cut.smi -o test_double_cut.fragments (where example_dataset_double_cut.smi are the SMILES: mmpdb index test_double_cut.fragments -o test_double_cut.mmpdb mmpdb transform --smiles 'c2ccc3c(c2)c(c(cc3NCc4[nH]c5c(n4)ccc(c5F)F)c6ccccc6)c7ccccc7' test_double_cut.mmpdb --explain Processing fragment Fragmentation(num_cuts=1, enumeration_label='N', variable_num_heavies=6, variable_symmetry_class='1', variable_smiles='c1ccccc1', attachment_order='0', constant_num_heavies=29, constant_symmetry_class='1', constant_smiles='c1c(-c2ccccc2)cc(NCc2nc3ccc(F)c(F)c3[nH]2)c2ccccc12', constant_with_H_smiles='Fc1ccc2nc(CNc3cc(-c4ccccc4)cc4ccccc34)[nH]c2c1F') Thanks for your help! |
Hi!
I am getting an unexpected generated molecule when using double cuts. The transformation, which in this case is the linker, is attached the wrong way around (details below). I would really appreciate your help :)
Using the two following molecules to create an MMPDB
c1ccc(cc1)c2ccc3c(c2)c(c(cc3NCc4cn5cc(cnc5n4)F)c6ccccc6)c7ccccc7 test1
c1ccc(cc1)c2ccc3c(c2)c(c(cc3NCc4[nH]c5c(n4)ccc(c5F)F)c6ccccc6)c7ccccc7 test2
to then generate new molecules from the molecule to improve
c2ccc3c(c2)c(c(cc3NCc4[nH]c5c(n4)ccc(c5F)F)c6ccccc6)c7ccccc7
the following molecule is proposed
Fc1cn2cc(CNc3cc(-c4ccccc4)c(-c4ccccc4)c4ccccc34)cnc2n1
while the expected generated molecule would be
c2ccc3c(c2)c(c(cc3NCc4[nH]c5c(n4)ccc(c5F)F)c6ccccc6)c7ccccc7.
Please let me know if any other information is needed to better understand my issue.
Thanks a lot,
Alice
The text was updated successfully, but these errors were encountered: