How to get the chemical information for each rule environment with a given rule? about mmpdb HOT 5 CLOSED

Dear Jen,

If I understand you correctly, you would like to see something like a SMILES or SMARTS for a given environment. This is currently not stored in the DB, so there is no direct way to access it.

A workaround could be that for a given rule_environment, you extract one sample pair and display this, together with the radius (as integer). If the from_smiles and to_smiles are known, the environment for a given radius can then be inferred from looking at the pair.

Please note that the environment is grouped via the Morgan fingerprint, rooted at the attachment atom. The Morgan Fingerprint for each atom amongst others features the number of neighboring atoms. Therefore, a radius 1 Fingerprint in reality is more a kind of radius 1.5 fingerprint. For this reason, if you wanted to depict the Morgan Fingerprint environment as a SMILES, this would get be rather complicated.

Does this answer your question?

from mmpdb.

Hi Christian,

Thanks for the reply. Yes, I would like to see the SMILES or SMARTS for a given environment. Good to know there is a workaround. I think it should work.

I will give a try to extract one example pair/or all pairs (if possible) from a given rule_environment. Would you mind give an example of how to do that? Is there any existing functions in the package I can use?

Thanks,
Jen

from mmpdb.

Hi Jen,

there is no existing function that does exactly what you ask. But if you are coding already anyway, and you know the rule_environment_id for which you want to see a pair, you could try something like this:

sql = """
SELECT pair.id,
lhs_compound.clean_smiles, rhs_compound.clean_smiles,
lhs_compound.public_id, rhs_compound.public_id
FROM pair,
compound as lhs_compound, compound as rhs_compound
WHERE pair.rule_environment_id = ?
AND pair.compound1_id = lhs_compound.id
AND pair.compound2_id = rhs_compound.id
LIMIT 1
"""
c = cursor.execute(sql, (environment_id, ))

Disclaimer: The concept of the statement should work, but I have not tested it and my python has gotten a little rusty lately, you may need to fix some of the commas etc. Please let me know if this works.

from mmpdb.

Hi Christians,

Thank you so much for the details. I will code the query and try to grab the pair information. I will keep you posted how it works.

Best, Jen

from mmpdb.

Hi Christian, It works as expected!

Thanks, Jen.

from mmpdb.