Comments (3)
Hey @LivC182, thanks for your interest in our code! Yes isomericSmiles is set to False
although this name can be somewhat misleading since it only affects stereoisomers (see rdkit doc). So for example this wouldn't preclude the differentiation between pentane and isopentane since these have different connectivities.
However, your point is valid in the more specific case of stereochemistry. This is an area for future work. I can only speculate on the authors' reason for excluding stereochemistry:
(a) for simplification since @@
and @
encode stereochemistry relatively (different canonicalization can flip this sign), and therefore does not encode the absolute (R/S) stereochemistry of the molecule. Therefore its not good enough to merely encode single and double @
individually)
(b) because many scoring functions (such as Morgan fingerprints which are used for the rediscovery and similarity benchmarks) also ignore stereochemistry by default
from guacamol.
HI @JoshuaMeyers,
I completely agree with all your points. I just wanted to make sure if isomers (except structural ones) are taken into account. If you have carbohydrate entries in your dataset not taking isomers into account might be problematic.
galactose_smiles = ''C([C@@H]1[C@@H]([C@@H]([C@H]([C@@H](O1)O)O)O)O)O''
glucose_smiles="C([C@@H]1[C@H]([C@@H]([C@H]([C@H](O1)O)O)O)O)O'"
gal = Chem.MolFromSmiles(galactose_smiles)
glc = Chem.MolFromSmiles(glucose_smiles)
print(Chem.MolToSmiles(gal, isomericSmiles=False))
'OCC1OC(O)C(O)C(O)C1O'
print(Chem.MolToSmiles(glc, isomericSmiles=False))
'OCC1OC(O)C(O)C(O)C1O'
from guacamol.
Yes this is true. In the GuacaMol v1 SMILES training dataset, there are in fact 45 occurrences of the SMILES you give as an example (as substrings of larger molecules). Thanks for highlighting this case.
from guacamol.
Related Issues (15)
- The hash value of the file is inconsistent with the hash value given in the code HOT 1
- TypeError: generate() missing 1 required positional argument: 'self'
- Error while assessing distribution learning benchmarks - FCD metric HOT 2
- Unknown segfault
- Support for input files of generated molecules
- how to test my Distribution-Learning Benchmarks?
- Something wrong in ranolazine_mpo() ?
- ChemNet file name has changed in FCD version 1.2
- Is it possible to evaluate a distribution without the models
- Latest scypi does not Support histogram HOT 1
- How the FCD value changes with the sample size of the reference molecule set and the padding length
- scipy imread HOT 1
- Now available on conda-forge HOT 2
- Support fcd_torch
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from guacamol.