openmolecules.org

 
Home » DataWarrior » Functionality » Calculating similarity to one query compound (molecular similarity)
Re: Calculating similarity to one query compound [message #1902 is a reply to message #1889] Mon, 22 May 2023 21:55 Go to previous message
thomas is currently offline  thomas
Messages: 715
Registered: June 2014
Senior Member
SMARTS and IDCodes use slightly different concepts. Therefore, a one-to-one conversion is not possible. Not all atom/bond query features, which are available in SMARTS, do exist in OpenChemLib (i.e. DataWarrior) molecule fragments and vice versa. DataWarrior does not support recursive substructures and SMARTS don't support 'exclude groups'. Most query features are translated by DataWarrior close enough, but not necessarily in an exact way. For instance, atoms that are considered aromatic in the Daylight (SMARTS) world are not necessarily considered aromatic in OpenChemLib. And OpenChemLib distinguished aromatic and delocalized atoms. The latter is not existing in SMARTS. During the recent year DataWarrior's SMARTS parser and creator were improved to closer translate features, but it will never be 100%, because the underlying concepts don't match. Another reason is that SMARTS don't support enhanced stereo recognition.

The conversion is done in DataWarrior this way: If you paste a table that contains SMARTS (or open a text file) then DataWarrior should recognize that and a create chemical structure (reaction) column from it, which is column containing IDCodes and is tagged to contain chemical information. To create a SMILES column (or SMARTS, if a structure column contains query features) call "Chemistry->From Chemical Structure->Add SMILES Code".

Your question referred to reactions rather than molecules and I realize that there is no "Chemistry->From Chemical Reaction->Add SMILES Code". You can, however, use "Copy From Reaction->Reaction As->Reaction SMILES". I could easily add the functionality "Chemistry->From Chemical Reaction->Add SMILES Code" with the limitations discussed above, if that would help?

[Updated on: Mon, 22 May 2023 21:59]

Report message to a moderator

 
Read Message
Read Message
Read Message
Read Message
Previous Topic: Add reaction from structure
Next Topic: 'Search ChEMBL Database' function
Goto Forum:
  


Current Time: Fri Nov 22 09:41:04 CET 2024

Total time taken to generate the page: 0.04146 seconds